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ome Further Investigations of Ionospheric Cross-Modulation 


By I. J. SHAW 
Late of the Cavendish Laboratory, Cambridge 


Communicated by W. L. Bragg; MS. received 5th Fune 1950 


BSTRACT. "The work in this paper forms the continuation and development of the 

experiments described in a paper by Ratcliffe and Shaw published in 1948. Measurements 

/ of the phase and amplitude of the modulation transferred to an unmodulated wave by a 

modulated disturbing wave in the ionosphere are used to provide information about that 
part of the ionosphere in which the cross-modulation occurs. 

‘ By measuring almost simultaneously the height of reflection of the wanted wave and the 
characteristics of the transferred modulation, a value of 1-4 x 10® per second for the collision 

¥ frequency of electrons with molecules at a height of 92 kilometres has been deduced. This 

measurement has been made possible by the use of a very sensitive automatic apparatus 
which is briefly described. 

y Special attention has been paid to the “ gyro-resonance effect ’ in cross-modulation first 

Ppointed out in a theoretical paper by Bailey in 1937. The experiments appear to show that 

$the predicted enhancement of the cross-modulation when the disturbing station radiates 

on the local gyro-magnetic frequency does not exist. 

A modification of the original theory is given which is of great use in the calculation of the 

ccross-modulation to be expected for a given pair of transmitters. ‘This modified theory is 

“used in the discussion of the gyro-resonance effect and of the observed decrease of the cross- 

{modulation at dawn. The agreement between the calculated and observed magnitudes of 

\cross-modulation is stressed in view of the fact that the calculations involve a knowledge of 

jah absorption of the wanted wave incident obliquely upon the ionosphere. 

/ 

| 


' See NOD Cre! ON 


N a previous paper, which will be referred to here as I, Ratcliffe and Shaw 
(1948) described experiments made to investigate the phenomenon of 
_ A ionospheric cross-modulation * and indicated how it was possible to deduce 
isome important facts about the ionosphere from the investigation. ‘This 
possibility has been discussed also by Huxley, Foster and Newton (1948) and 
Huxley and Ratcliffe (1949). 

It is the purpose of this paper to outline a modification of the apparatus f 
described in I, and to discuss some of the results of importance in the study of 


i a 


__ * The phenomenon has previously been called ‘ wave-interaction ’ ‘Luxembourg effect ’. 
Huxley and Ratcliffe (1949) have suggested the more precise description “ionospheric cross- 
‘modulation ’. 

t+ A full description of the experimental technique will be given elsewhere. 
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the ionosphere which have been obtained with it. The present paper fae a 
sequel to I but, where appropriate, the nomenclature has been altered to conform 
with that suggested by Huxley and Ratcliffe. | 

In I, it was shown how it is possible, by measuring the phase of the transferred 
modulation at a series of modulation frequencies, to determine the collision 
frequency in the ionosphere at a height which could be deduced from ti? 
experimental data. The deduction of the height required a knowledge of the 
height of reflection of the wanted wave, and as this could not be measured it 
was assumed to be 100 kilometres. Any error in this assumed height would give 
rise to a corresponding error in the estimated height for the observed value 
of the collision frequency. In §3 of the present paper a series of measurements 
of cross-modulation is described in which the height of reflection of the wanted 
wave was measured simultaneously. This leads to what is considered to be a 
reliable value of the collision frequency at a known height. 

In the experiments of I, the phase and amplitude of the transferred modulation 
were determined by a comparatively slow method of visual observation so that 
about 20 minutes were required to make a full ‘ sweep’ of modulation frequency 
from 60c/s. to 2,000c/s. During this time the ionospheric conditions were 
often found to have changed so that considerable inaccuracy was introduced. 
In §2 an automatic measuring equipment which considerably speeds up the rate 
of observation is briefly described. With this apparatus it is possible to obtain 
in four minutes a photographic record of the phase and amplitude of the 
transferred modulation over the required range of modulation frequencies. 
Errors due to alterations in the ionosphere during the period.of an experiment 
are thereby considerably reduced. ‘The automatic equipment is also more 
sensitive than visual observation so that much smaller depths of modulation 
can be measured. 

In I, the original theory of Bailey and Martyn which has so far proved adequate 
to explain the observations was re-stated in a form more closely related to 
Appleton’s magneto-ionic theory and this formulation, with some improvements 
in nomenclature, was extended by Huxley and Ratcliffe. If it is required to use 
the theory given in these papers to calculate the cross-modulation expected with 
a given pair of transmitting stations, it is necessary first to assume a distribution 
of electron density with height. It is shown in §5, that with certain simplifying 
assumptions it is possible to deduce a general expression which relates the 
amplitude of cross-modulation to the total power absorbed from the disturbing 
wave, whatever the distribution of the electrons may be. This expression is 
used in a discussion of the following points: 

(i) It has sometimes been suggested that markedly increased cross-modulation 
should be observed when the disturbing wave has a frequency equal to the 
‘ gyro-frequency’ of the magneto-ionic theory. Experiments described in §3 
have not shown this increased cross-modulation and the theory of §5 does not. 
lead us to expect it. | 

(i) The observed decrease of cross-modulation at dawn is discussed in terms 
of the theory. 

(iii) The theory in its simple form is found adequate for the calculation of 
the absolute magnitude of the cross-modulation observed with different pairs of 
transmitting stations. The implications of this result for the theory of the 
absorption of radio waves at oblique incidence are discussed in § 6. 


Further Investigations of Ionospheric Cross-Modulation 3 


§2. EXPERIMENTAL TECHNIQUE 


The general form of the experiments has been similar to that described in I 
and is illustrated diagrammatically in Figure 1. A ‘ wanted wave’ is transmitted 
from a station W and received, via the ionosphere, at a receiver R. A ‘ disturbing 
station’ D radiates a modulated wave which is received at R by the direct path 
DR over the ground. At some region B in the ionosphere modulation is. 
transferred from the disturbing wave to the wanted wave. The experiment 
consists in comparing, at R, the phase of the modulation received via the direct 
path DR with that of the transferred modulation received on the wanted wave, 
jover the range of modulation frequencies 60c/s. to 2,000c/s. Measurements. 

are also made of the magnitude of the transferred modulation. 


Figure 1. Schematic diagram of experimental arrangement of transmitters and receiver. 


I 


Wanted 
Station 
Receiver 


Phase 
Sensitive 
id Rectifier 


Disturbing Phase 
Station Sensitive 
Receiver Rectifier 


: Figure 2. Simplified schematic diagram of automatic recording apparatus. 


With the valuable help of Mr. S. W. Falloon equipment has been developed 
| with which the amplitude of the transferred modulation and its phase relative 
}to that of the direct wave from the disturbing station can be rapidly recorded 
| photographically. As this apparatus will be described in detail elsewhere, 
it is sufficient here to describe it in outline only. For this purpose reference 
t should be made to the block schematic diagram of Figure 2. 

The basic element of the apparatus is the phase-sensitive rectifier circuit. 
The cross-modulation signal T cos (wt+¢) and the reference signal M cos wt 
/are mixed in the rectifier, giving an output, when M is large compared with 7, 
' proportional to Tcos¢. From a second phase-sensitive rectifier, driven with the 
same cross-modulation signal and with the reference signal shifted 90° in phase, 
A-2 


4 I. F. Shaw | 
an output proportional to Tsing is obtained. The two output voltages | 
proportional to Tcosf and Tsing are applied to the x and Y plates of a : 
cathode-ray tube. The spot then moves from its rest position to a point defined 
by T and ¢, and this movement is recorded photographically as a vector. 
Figure 3 shows the type of record obtained. 


Taam recetttatteae OCCT ceceet ecco MAAN 


i . Atypical rd of the phase and amplitude of the transferred modulation at frequencies 
oe ae 60, 120, 300, 500, 700, 1,000, 1,250, 1,500 and 2,000 c/s. 


With this apparatus it is possible to obtain about twelve readings of the phase 
and amplitude of the transferred modulation in 30 seconds with an accuracy far 
above that obtained by the visual method used before. A complete variation of | 
modulation frequency in steps from 60 c/s. to 2,000c/s. can be made in four 
minutes. It is possible, on account of the great discrimination of the phase- 
sensitive rectifier against noise, to measure transferred modulation of the order 
of 0-1% in the presence of noise and adjacent channel interference. 


§3, THE COLLISION, FREQUENCY, Adv vA KEN OWEN SER iG ria 


In I it was shown that the phase difference ¢ between the modulation received 
via the path DR and that received via the path DBR (Figure 1) is related to the 
modulation frequency by the equation 


anf 2nfd 
= =) 
¢@ =tan Gt rie 


where d is the difference of the lengths of the two paths, f the modulation 
frequency, v the collision frequency of the electrons in the region of cross- 
modulation, and G the energy loss factor. 

From observations of ¢ at different modulation frequencies it is possible to 
determine both vy and d. If the trajectory of the wanted wave can be determined, 
it is then possible from geometrical considerations to determine the point on the 
trajectory at which cross-modulation occurs. 

The actual path of the wanted wave through the ionosphere is not known, 
but the equivalent triangulated path can be determined by observing the delay 
between pulses from the transmitter received via the ionosphere and those 
received directly along the ground. It is important to notice that the assumption 
is made here, and also in the phase-frequency equation quoted above, that the — 
velocity of propagation of the waves in the ionosphere is equal to the velocity in . 
free space, and that the wanted wave follows a linear path from the transmitter _ 
to the receiver with a sharp change of direction at the virtual level of reflection. — 
‘The consequences of these assumptions will be examined in the Appendix. 

The transmitters used in these experiments were normal B.B.C. stations — 
situated at Stagshaw and Ottringham (see map of Figure 7). Stagshaw radiated 
the wanted signal of 1,050 kc/s. and Ottringham the disturbing wave on a frequency 
of 167 kc/s.* 

The distance from Cambridge where the receiver was situated, and the 
power of the Stagshaw station were such that a ground wave which was, in general, 
small compared with the wave reflected from region E at night, reached the receiver. 


* I am indebted to the B.B.C. for allowing these transmitters to be used for the experiments. 
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t was therefore possible to determine the equivalent path of the wanted wave. 
t the same time, the intensity of the ground wave was not sufficiently high to: 
pset the phase measurements (see Appendix of paper I). 

Experiments were so arranged that a short pulse transmission from the 
anted station preceded and followed each determination of v and d from the 
odulation run. In this way the equivalent height of reflection of the wanted 
ave was determined at intervals of about four minutes throughout the night. 
Figures 4(a) and 4(4) show the values of Gv/2m and the heights of the region 
f cross-modulation throughout the nights of 28th April 1949 and 5th June 1949. 
A striking difference in the curves is immediately obvious, a difference which 
ill be discussed later. On three occasions on which measurements of collision 


Gy/ 20 


Height (km) 
Oo So 


Gy/2m 


Height (km) 
ee 


G.M.T. 


Figure 4. Values of Gv/2z7 and heights of regions of cross-modulation on (a) a quiet night, and 
(6) a disturbed night. 


» frequencies and the associated heights of the region of cross-modulation have 
» been made throughout the night, the following average results were obtained : 


) 28th April 1949.  Gyv/27=307, v=1-48 x 108, height =91-7 km. 
5th June 1949. Gyv/27=552, v=2-66 x 108, height = 86-1 km. 
21st June 1949. Gv/27=612, v=2-96 x 108, height = 86:5 km. 


Of these results, the first is undoubtedly the most accurate since the measure- 
* ments on this occasion showed that conditions remained constant, the mean 
deviation in the value of the collision frequency being +15%. As far as the 
- height measurement is concerned, although the equivalent height can be measured 
) from the observations to an accuracy of within +1°% it does not necessarily 
represent the actual height of the region of cross-modulation. Considerations 
- elaborated in the Appendix indicate, however, that although the real height may 
be a little lower than the equivalent height, the difference is not greater than 
2-5 kilometres. 
| 
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The last two results are practically identical, and if we take their average and 
combine with the first result, the scale height of the atmosphere at this level 
may be deduced. This calculation gives H =8-5km. a result in good agreement — 
with those obtained by other methods. On account of the rather small difference 
in the heights at which the two measurements were made no great accuracy 1s 
claimed for this determination of scale height. 


§4. VARIATIONS IN THE MEASURED VALUE OF yp 


On the majority of occasions on which measurements were made throughout 
the night, graphs of Gv/27 and heights of the region of cross-modulation had the 
form shown in Figure 4(b). Only very seldom were conditions as stable as those 
shown in Figure 4(a), when measurements over an extended period of time 
were constant or showed only a smooth variation. 

It might at first sight be considered that variations in the value of Gy/2a 
mean merely that the level of cross-modulation had changed, due, perhaps, to” 
some irregularity or inhomogeneity in the ionosphere. If this were the case, the 
values of Gy/27 should be related to the measured heights by the usual scale 
height formula. No such relation is found and indeed, on many occasions, an 
increase of Gv/27 is associated with an increase of height and vice versa. 

A change in the value of Gv/27 and in the height of the region of cross- 
modulation could be brought about by the displacement of one point of the 
phase-frequency curve by a large amount from its normal position. Such a 
change of phase could easily be produced by selective fading although it is unlikely 
that this effect would persist for 30 seconds (i.e. the time of measurement on each 
modulation frequency). In any case, it was found that on the majority of 
occasions the phase-frequency curve was of the normal type, with all the points 
lying nicely on the curve. 

To help in relating this type of disturbance to other disturbances in the 
ionosphere a ‘ disturbance figure’ for measurements of collision frequency has 
been defined. ‘The disturbance figure k, is obtained by multiplying the maximum 
spread of values of Gv/27 by the number of times the curve relating Gyv/27 and 
time changes direction, and dividing by the total number of points in the interval 
considered. ‘This somewhat arbitrary process leads to values of the disturbance 
figure which fit in well with intuitive ideas of the disturbance on different 
occasions. 

Attempts have been made to correlate k, with: (1) the magnetic K indices, 
(2) the disturbance of propagation conditions of long waves (data were kindly 
supplied by Dr. K. Weekes of the Cavendish Laboratory), (3) the incidence of 
sporadic E ionization, (4) meteoric activity. 

As far as (1), (2) and (3) are concerned, there appears to be no correlation. 
With meteoric activity, it is not easy to draw any definite conclusions on account 
of the rather small amount of meteoric data available for the nights on which 
experiments were carried out. It may, however, be significant that on the very 
few occasions on which conditions were not disturbed, no meteoric showers were 
forecast, whereas on all other occasions showers of varying intensities were 
expected to occur. 

. In the absence of further information, it is not safe to draw definite conclusions. 
It is hoped that more definite results, one way or the other, will be obtained when 
the analysis of certain records of meteoric activity has been completed. 
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§5. THEORETICAL CALCULATION OF THE ABSOLUTE MAGNITUDE 


OF CROSS-MODULATION 
The theory of ionospheric cross-modulation proposed originally by Bailey 


jy and Martyn (1934) and adapted recently (Ratcliffe and Shaw 1948, Huxley 

et al. 1948, Huxley and Ratcliffe 1949) to conform with the nomenclature of 
» Appleton’s magneto-ionic theory, gives an expression for the modulation 
| impressed by a disturbing wave of known field strength upon a wanted wave 
} passing simultaneously through a medium of uniform electron density and 


collision frequency. In attempting to apply this theory to the ionosphere we 


* must first make assumptions about the variation of electron density and collision 


frequency with height. We must remember also that the field strength of the 
disturbing wave is itself reduced in its passage through the ionosphere in a way 
which depends upon the assumptions made. We shall then be in a position to 


* calculate the modulation transferred at each level of the ionosphere. If the 
* calculation is done in this way it becomes apparent that the reduction in the 
+ field strength of the disturbing wave in any narrow layer in the ionosphere 1s 


closely related to the amount of modulation transferred in that layer and that, 


* in a straightforward case the amount of transferred modulation can be related 


to the reduction in amplitude of the disturbing wave without a detailed knowledge 


' of the distribution of electron density. ‘This important fact was pointed out 


to me by Mr. J. A. Ratcliffe and I am indebted to him for indicating the following 


u simple way of making the calculation. 


Let T, be the coefficient of modulation transferred from a disturbing wave 
to a wanted wave at a modulation frequency approaching zero. Suppose F is 
the power flux per unit area of tne disturbing wave at the region of cross- 


» modulation. The energy left per second in an elementary slab of unit area and 


thickness 6x is —(0¥/0x)5x and this energy is taken up by the electrons in the 
slab in the process of absorption. If the electron density is N and if P is the power 


+ expended on each electron in the slab by the disturbing wave, then 


" oF loF 
PNédx = — A OS 1.6, P= — Nox’ 

When the disturbing wave has been acting long enough to establish 
equilibrium conditions a stage is reached at which, on the average, the energy 
gained by an electron per second from the disturbing wave is equal to the energy 
lost to molecules in collisions. We then have P=vGAQ where » is the collision 
frequency of electrons with molecules. G is the energy loss factor which defines 
the average proportion of its excess energy lost by an electron in a collision. 
AQ is the difference between the energies of the electrons and the surrounding 
gas molecules. 

The electron energy Q is given by 


where m is the mass of the electron, v is its velocity, / the mean free path and 
v the collision frequency. Then AQ=ml?vAv and the change Av in the collision 
frequency of the electrons is given by 


Ca LOO as 
*= ‘ml vGmlv  3GkR0 
since O=4 mi? = $k, 


where k is Boltzmann’s constant and 6 is the temperature in degrees Kelvin, 


8 I. ¥. Shaw 


From a consideration of the exponential law describing the absorption ofa 


radio wave in the ionosphere, it can be seen that the element dT) of cross- 


modulation impressed upon a vertically travelling wanted wave by virtue of the 


change in collision frequency is given by dT, =Axw5x where ky is the coefficient 
of absorption of the wanted wave. 

If itisassumed that pw? >v2*, then «yw may be written in the formky =a(pw)Nv 
where «(py) is a function of the angular frequency of the wanted wave and the 
natural gyro-frequency, and depends upon the conditions of propagation. 

Then over the region 6x 
NP 


3GR0°* 


dT) =4( Pw) NAvdx =a(y) 


lee Ug, 
Pw) scp ox 


It should be noticed that, in this equation, neither of the quantities N and v 


= =a 


appear explicitly, and that the amount d7, of modulation transferred is related _ 


directly to the change (@7/dx)8x of the power flux of the disturbing wave as it 
passes through the element of ionosphere concerned. It is therefore possible to 
integrate and to obtain an expression for the total amount of transferred 
modulation in the passage through a layer of finite thickness in terms of the 
total decrease of power flux of the disturbing wave, and without any further 
assumptions about the variations of N and v through the layer. We then obtain 


a(pw) (OP 


3GRO J, Ox 


T)=|dT)= — 


_ a( pw) 
= 3686 (Fi=P 3) = eee (1) 


This result holds for a wave travelling upwards through the ionosphere. Ifa 
downcoming reflected wave is added, it can be shown that our result still holds 
if Y,; represents the power flux in the wave as it begins its upward journey, and 
FY, that as it returns to the earth. 

Equation (1) can be used to calculate the modulation transferred when both 
waves travel vertically and traverse the same region. But we know that the 
heights of reflection of different frequencies are different, and we must consider 
in more detail which parts of the region the waves traverse together. Let us call 
the wanted wave W and the disturbing wave D, and let their vertical paths be 
represented diagrammatically in Figure 5. We shall assume that D has the lower 
frequency and is reflected from the lower level. We further assume that for each 
wave the absorption is important only above the level where v=(pw-+-P,) or 
v=(Pp+Ppz,)t (the shaded regions in the figure). It is clear in this case that the 
whole of the power absorbed from D will be effective in producing cross- 
modulation, and equation (1) applies. Absorption of W in regions where D is not 
absorbed (e.g. above the level of reflection of D) will not produce cross- 


modulation, and is equivalent merely to a decrease in the power of the wanted 
station. 


*The assumption Pw’ >? is appropriate to the ordinary magneto-ionic component in 
quasi-transverse Propagation. In the quasi-longitudinal case, the corresponding assumption is 
(Dy tpyeS> ve. 

+ This point is discussed by Farmer and Ratcliffe (1935). 


} 


| 


| 
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If the frequency of D is greater than that of W, D is reflected from the higher 
level and the lettering in Figure 5 is reversed. The term (7, —,) in equation (1) 
does not now refer to the total power lost from D in its passage through the 
ionosphere. If equation (1) is to be used in this case, we are forced again to make 
assumptions about the electron distribution to decide what proportion of the total 
power lost by D is effective in producing cross-modulation. Alternatively, by 
comparing the observed and calculated values of cross-modulation in cases such 


# as these we may obtain an estimate of the ionization density gradient in the region. 


Considerations of this nature lead to the deduction of ionization gradients 
discussed in § 9. 

If D and W are incident obliquely upon the ionosphere the situation is more 
complicated. ‘The oblique incidence case is represented diagrammatically in 
Figure 6. jase 


VDP, PSS 


+ Figure 5. The relative positions of regions of appreciable absorption of wanted and disturbing 


waves for normal incidence. 


| Figure 6. The relative positions of regions of appreciable absorption of wanted and disturbing 


waves for oblique incidence. 


The value of p;, to be used for each wave depends upon the angle between the 
direction of propagation of the wave and the earth’s magnetic field. For the 
wanted wave, we shall call the appropriate value of p, for the upgoing portion of 
the trajectory p,, and for the downcoming portion py. For the disturbing wave 
we shall call the appropriate value of p,, py’ for the wave producing cross- 


modulation in region A and p,’ for that producing cross-modulation in region B. 


Consider first region A. If the frequencies and directions of propagation of 


the waves are such that hy >hp,, and (Pwt+Pa)>(Pp+Pa’) we have the case where 
the region of absorption of W completely encloses the region of absorption of D. 
We shall assume, as a reasonable approximation, that the trajectory of W is 
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straight over the region of cross-modulation, so that equation dTy=Akyw9dx is | 
replaced by dT) =Akw5x seciw and equation (1) becomes 


a ci 
Lo= Nah Sec ly (FP, —FPa)a- 


Similarly for region B we have 


a p , 
fhe = ee sec ty( Fy a Pr). | 
In practice it is generally found that since the disturbing station is somewhat 
to one side of the plane of propagation of the wanted wave, that the ‘ angles of fire’ 
of the disturbing station to the regions A and B are not greatly different. Hence © 
(P,—Pz2), will be nearly equal to (7,—Y,), and the final equation may be 
written 


a +o : 
La Sg oe (FP, — Pe) secty. 


It might be argued that in the case of oblique incidence the disturbing wave 
after reflection does not pass again through the region of cross-modulation as is 
the case with normal incidence. Hence for oblique incidence, the term (7, — PY.) 
should not refer to the total loss of power of the disturbing wave but only to that 
absorbed on the upward trajectory. This is, in fact, true, but we have reason to 
believe that the maximum loss of power per unit distance occurs before the 
disturbing wave reaches its level of reflection (see § 9) and that the power remaining 
for absorption on the downward path is a relatively small fraction of the incident 
power. It must be remembered also that although the particular trajectory 
which we have considered will not pass again through the region of cross- 
modulation, a wave which leaves the disturbing transmitter at a slightly greater 
angle of elevation will, after reflection, pass through the region of cross- 
modulation. It appears therefore, that as a good working approximation, we may 
take (7,—P,) in the oblique incidence equation to refer to the total power loss 
on the double journey through the absorbing region. 

There may, of course, be conditions for which hy >Ay or for which the lower 
edge of the absorption region for the disturbing station is below the corresponding 
level for the wanted station. A knowledge of the angles of incidence of the two 
waves on the ionosphere and of the angles between the directions of propagation 
of the various waves and the earth’s magnetic field is sufficient to determine, for 
give values of py and pp, whether the conditions illustrated in Figure 6 apply. 
In the majority of cases which have been investigated experimentally the 
conditions illustrated do apply. 

The function (py) is not a complicated one and has, in fact, been discussed 
by Booker (1935) although not explicitly. Booker has shown that under 
appropriate conditions, the propagation of a radio wave in the ionosphere may 
adequately be described by the use of two approximations to the general magneto- 
ionic formula. ‘These are the quasi-longitudinal (Q.L.) and quasi-transverse 
(Q.T.) approximations, and conditions are given for which each approximation 
holds. Applying Booker’s criterion, it is found that in England for the frequencies 
in which we are interested (i.e. up to about 1 Mc/s.) the propagation may generally 


en a 


—= 
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be described by the Q.L. approximation except for cases where the direction 


| of propagation is almost exactly perpendicular to the earth’s magnetic field, and 


provided that the frequency is not too near the critical frequency. For example, 
for a 200 kc/s. wave, the Q.L. approximation holds until the angle between the 
direction of propagation and the earth’s magnetic field is 89°; for 1-4 Mc/s., the 


{ Q.L. approximation holds up to 75°. Booker shows also that for wave frequencies 


Jess than the local gyro-magnetic frequency (approximately 1-3 Mc/s.) the 
ordinary component only is reflected, and we need consider, therefore, only the 


} ordinary component of the downcoming wanted wave in our work. 


Applying Booker’s approximations, the function «( py) becomes 


| 27e” 1 
" 1 dibs See 
( ) Q a( Py) €)mc (pw +pr) at p2 
27e” 1 


Q) Qt. (Py) = 


eymc py? +v?- 


If, in addition, the effective radio frequency is large compared with the collision 
_ frequency (as is nearly always the case) further simplification may be effected : 


277e? 1 
(1 id Des -)= — _— 
( ) Q a( Py) «mc (pwt pr)? 
Q) avr. 2s ee 
¥ epmec pw - 


j $6. COMPARISON OF CALCULATED AND OBSERVED MAGNITUDES 


In view of the simplicity of equation (1) and its extension, equation (2), for 
oblique incidence, it is no difficult matter to make a comparison of the observed 
and calculated values of cross-modulation for different pairs of stations. ‘This 
comparison is tabulated. ‘The following assumptions have been made in the 
calculation : 

(i) The disturbing wave is completely absorbed in the region of cross- 
modulation. If the reflection coefficient is about 4, about 90°% of the incident 
energy is absorbed so that this assumption will not lead to any great error. 

(ii) The wave radiated from the disturbing transmitter is assumed to have a 
cosine angular distribution of intensity in the vertical plane. ‘This is probably 
very nearly correct in the case of the long wave stations but may not be correct 
for the medium wave disturbing stations. 

(iii) The region of absorption of the disturbing wave (as discussed in §5 
above) is contained entirely within the region of absorption of the wanted wave. 
‘This is generally, but not invariably, the case in our experiments. 

(iv) The temperature of the ionosphere is taken as 300° K. 

(v) The energy loss factor G is taken as 1-3 x 10~? according to a recent 
determination by Huxley and Zaazou (1949). 

In making the calculation it is first decided whether the propagation of the 
wanted wave through the region of cross-modulation can be described by use 
of the Q.L. or Q.T. approximations. If the former is the case, the required 


values of p; for the upgoing and downcoming parts of the trajectory are 


calculated from the gyro-magnetic frequency and the angles between the earth’s 
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magnetic field and the appropriate parts of the trajectory. The two values of | 
a(Pw) are then found. 


Disturbing Frequency Power Wanted Frequency Magnitude of Ty ( WA 
Station (ke/s.) (kw.) Station (ke/s.) Calc. Obs. 
Ottringham 167 120 Westerglen 767 SOF) 6 
Lisnagarvey 1050 1:6 1 
Droitwich 1050 1-4 3 
Ottringham 1122 0-7 0:5 
Stagshaw 1050 1:8 2 
Droitwich 200 180  Lisnagarvey 1050 4-3 4 
Droitwich 1050 1-2 iho) 
Washford 804 2-6 2 
Decca 70 ALD) <0-1 
Rubgy 68 80 Lisnagarvey 1050 they) <0-1 
Moorside Edge 668 100 Lisnagarvey 1050 3-0 2 
Ottringham * 167 120 Droitwich 200 13 10-20 


* The observations on this pair of stations were made in Copenhagen by Professors Rybner and’ 
Nielsen, and I am indebted to them for permission to quote this result. vA 
The map of Figure 7 shows the approximate positions of the stations used; 
each station is marked with the initial letter of its name, and the Decca Navigation 
Station mentioned is situated at Warwick. The receiving station was at 
Cambridge in all cases. The Post Office transmitter at Rugby is constructed 
for single sideband telephony. !t was used as a disturbing station in two different 


Kilometres 


Figure 7. Map giving positions of transmitters used in the experiments. 


ways : in the first method, an audio-frequency tone injected into the speech 
circuit was keyed on and off, thus producing an intermittent sideband simulating 
keyed c.w.: in the second method a 1,000c/s. audio tone was itself modulated 
at 50c/s., so that the sideband emitted was effectively a single carrier modulated 
at 50c/s. All the measurements, with the exception of those using Decca and 
Rugby transmitters were made with the automatic apparatus. 

‘The agreement between the calculated and observed values is seen to be quite 
good except in the two cases where very low frequency stations are concerned. 
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In the case of Droitwich 200kc/s. disturbing Decca 70kc/s. the disagreement 
may be due to the fact that the level of reflection of the 200 kc/s. wave is above 
that of the wanted wave (i.e. assumption (iii) above does not hold). For Rugby 
68 kc/s. disturbing Lisnagarvey, 1,050 kc/s. it is difficult to see why the calculated 
and observed values do not agree. 

The agreement in the majority of cases is of importance in view of the fact 
that the absorption of the wanted wave at oblique incidence is involved. The 
phenomenon of cross-modulation offers a powerful method of investigating the 
absorption of radio waves incident obliquely upon the ionosphere, uncomplicated 
by considerations of the loss of energy by scattering and similar processes. It 
appears from these results that the true absorption of a wave incident obliquely 
upon the ionosphere at night can be calculated by the straightforward application 
of the magneto-ionic formula when due attention is paid to the increase in the 
length of the path of the wave in the region of absorption, and to the change in the 
value of p;, brought about by the obliquity of the trajectory. 


: §7. VARIATION OF CROSS-MODULATION AT DAWN 


Since the absorption of radio waves in the ionosphere begins to increase at 
dawn it might be expected that the depth of modulation should also increase. 
‘Previous experiments by Ratcliffe and Shaw (1948) using Lisnagarvey 1,050 ke/s. 
as wanted station and Droitwich 200 kc/s. as disturbing station showed, however, 
that the depth of transferred modulation decreased at sunrise. This result is 
difficult to understand, and in order to gain more information on the subject 
experiments have been carried out on different pairs of stations on frequencies 
separated as widely as possible. Combinations of stations for which results have 
been obtained are: 


(1) Droitwich 200 kc/s. disturbing Lisnagarvey 1,050 kc/s. 

(2) Ottringham 167kc/s. disturbing Westerglen 767 kc/s. 

(3) Ottringham 1,340 kc/s. disturbing Westerglen 767 kc/s. 

(4) Moorside Edge 668 kc/s. disturbing Lisnagarvey 1,050 kc/s. 


All the results show a decrease in cross-modulation as sunrise approaches. 
The decrease begins about 80 minutes before sunrise at ground level, a time 
which corresponds to sunrise at a level of about 100 km. at this latitude in summer 
when the experiments were made. By the time of ground sunrise there is 
generally only an extremely small amount of transferred modulation remaining. 
It should be noticed that Huxley (1950) has recently published experimental 
evidence which appears to show, on the contrary, an increase in the depth of 
cross-modulation at dawn. In the light of our own experimental results we shall 
assume that there is a decrease of cross-modulation at dawn and shall consider 
possible theoretical reasons for it. 

In order that cross-modulation should occur, it is necessary that both the 
wanted and disturbing waves should be absorbed in the same region of the 
ionosphere. It might be considered then, that the decrease at dawn is due to the 
fact that the regions of absorption of the two waves no longer overlap. 

We have seen that the absorption of a radio wave of angular frequency p is 
most important between the level where (p+ );) =v and the level of reflection. 
In the cases of Droitwich 200kc/s. disturbing Lisnagarvey 1,050kc/s. and 
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Ottringham 167 kc/s. disturbing Westerglen 767 kc/s., the regions of absorption 
of the disturbing waves are contained entirely within the regions of absorption of 
the wanted waves and no alteration of the ionization density, or of its gradient, | 
at dawn can alter the disposition of the absorbing regions. It follows, therefore, _ 
that at dawn, as the absorption of the disturbing wave increases, the term 
(P,—P.) in equation (1) increases, and the cross-modulation should increase 
also. ‘The decrease at dawn cannot, therefore, be explained in these cases on the 
hypothesis that the regions of absorption of the two waves drift apart as dawn 
approaches. 

Were there appreciable absorption of the disturbing station below the region 
where (Pp+p,)=v, the effect would be to reduce considerably the effective 
power of the disturbing station penetrating to levels higher than this, and hence 
to reduce the cross-modulation. But absorption below the region where 
(Pp+Ppr) =v is independent of frequency (the coefficient of absorption propor- 
tional to N/v) and if such a region is responsible for the increase of absorption of 
the disturbing wave, then it must be equally effective for the wanted wave, and i 
therefore there must be an increase in the cross-modulation. 

By no mechanism of change of reflection height or absorbing levels does it 
seem possible to explain the dawn decrease of cross-modulation in all cases. 

It will be remembered that in the derivation of equation (1) the assumption 
was made that (Pw+ pz)? >v?. If this restriction is removed and if it is assumed 
that v is approximately constant over the region of cross-modulation the final 
equation becomes 


admes sl oe (Pwepn eas gee 
o egme 3GRO {(pw t+ pr) +r? (Pi= Fa). eae (3) 


From this equation it is seen that any variation of the collision frequency caused 
by the absorption of the disturbing wave at the level where v=(py +p,) cannot 
cause cross-modulation. 

The region in the ionosphere where cross-modulation is most important 
is that where the absorption per unit distance of power from the disturbing wave 
is a maximum. ‘The loss of power per unit distance at a given level depends 
upon (1) the coefficient of absorption at that level, and (2) the amount of power 
left in the wave when it reaches the level concerned. As dawn approaches, there 
is a general increase in the ionization density at all levels of the ionosphere and a 
change from a steep gradient of ionization in the lower levels of region E to one 
which is rather less abrupt. Both of these effects will cause the level of maximum 
absorption per unit distance to become lower (i) because the coefficient of 
absorption at a given level is increased, and (ii) because the overall attenuation is 
increased so that the amount of power remaining in the wave when it reaches a 
given level is decreased. 

If the change in the ionization density and gradient is sufficiently great to 
bring the level of maximum power loss of the disturbing wave down near to the 
level at which v=(py+ py), then, as shown by equation (3), the cross-modulation 
will decrease to a very small value. 

If this is the true explanation of the dawn decrease of cross-modulation, the 
values of Gv/27 should increase as the amplitude of the cross-modulation decreases. 
There are indications that this increase does take place, though it has not been 
established definitely in our experiments. The reasons for this are (i) a large 
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value of Gv/27 gives a phase-frequency curve that is nearly a straight line and the 
consequent analysis is difficult and inaccurate, and (ii) the associated sudden 


decrease in amplitude renders it all the more difficult to make the measurements 
accurately. 


§8. CROSS-MODULATION ON THE GYRO-FREQUENCY 


Bailey (1937, 1938) predicted that if the disturbing wave had a frequency 
equal to the local gyro-magnetic frequency there should be strongly enhanced 
cross-modulation brought about by the great increase in ionospheric absorption. 
Cutolo (1947) has obtained some experimental evidence in support of this idea, 
for he claims to have observed depths of transferred modulation of the order 
of 3°, with disturbing stations of 5 or even 1 kilowatt power. This result would 
correspond to depth of modulation of more than 100% with powers of 
100-200 kilowatts such as we use. 

Experiments have been carried out to look for an increased cross-modulation 
when the disturbing station radiated on the gyro-frequency. Westerglen 767 kc/s. 
was used as the wanted station, and Ottringham, 120kw., on 1,122, 1,312, and 
1,474kc/s. as the disturbing station. Over about 24 hours’ observation on 
different nights, no significant difference was found in the depth of transferred 
modulation. It was considered that the frequencies might have been too widely 
separated and that a sharp resonance effect could possibly have been overlooked. 

For this reason a more complicated experiment was planned and carried out 
with the cooperation of the B.B.C. The frequency of the disturbing station at 
Ottringham was varied from 1,325kc/s. to 1,355kc/s. in steps of 5kc/s. The 
radiated power was of the order of 100kw. On each radio frequency, tone 
modulation at 60c/s. was radiated for two minutes and a measurement of the 
depth of modulation made at two-second intervals. ‘Thus each measurement 
was the mean of about 60 readings. The whole run from 15kc/s. on one side 
of the calculated gyro-frequency to 15 kc/s. on the other side was repeated twice 
each night, and the whole experiment was repeated on four nights. 

When the depth of modulation for any particular run from 1,325 to 1,355 kc/s. 
was plotted as a function of frequency, it was found that there were quite large 


variations in the depth of modulation from frequency to frequency. But points 


of high cross-modulation did not repeat from one experiment to the next, and 
these variations must therefore be attributed to changes in the ionosphere. When 
the means of all observations were plotted as a function of the frequency of the 
disturbing station there was found to be practically no variation of the depth of 
cross-modulation with frequency. ‘The general level of cross-modulation in 
these experiments was about 3%. 

Engineers of the B.B.C., using an entirely different technique of measurement 
and situated at a different receiving point, also carried out measurements during 
these tests. Their observations agree completely with those obtained at 
Cambridge in showing no sign of a maximum at the gyro-frequency.* 

We now consider the matter from the theoretical standpoint. In the equation 


a 
Ty= Ho) (P,P) 


* I am indebted to Mr. H. L. Kirke, Head of the Research Department of the B.B.C., for 
arranging for these observations to be made and for permission to quote the results. 
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it is the term (P, — 2) which determines the behaviour of Ty as the frquency of the 
disturbing station is altered. For a disturbing wave of low frequency, say 
200 ke/s., the reflection coefficient of region E at night is of the order of 0:2. This 
means that about 96%, of the total incident power remains in the ionosphere. 

Now consider a disturbing wave of frequency equal to the local gyro-magnetic 
frequency. The coefficient of absorption per unit path in the ionosphere may be 
many hundreds of times greater than that of the low-frequency wave, but it is the 
term (P,—,) which is important, and this can never be greater than Y,. The 
cross-modulation with the disturbing station on the gyro-frequency can thus be 
only about 4°%, greater than that for the low-frequency disturbing station. ; 

We conclude, therefore, that there is not likely to be any great increase in the 
cross-modulation when the disturbing station radiates on the local gyro-frequency. 
This theoretical conclusion is in agreement with our experimental results. 

In his discussion of the phenomenon of gyro-resonant cross-modulation, 
Bailey made detailed calculations of a special hypothetical case where the disturbing 
wave was on the gyro-frequency and the wanted wave on a frequency of about 
50kc/s. and showed that the transferred modulation would be unusually great 
(20% for 100kw. power radiated by the disturbing station). It is of interest to 
examine Bailey’s result in detail, and to relate it to our conclusions. 

Examination shows that the large value of the cross-modulation is brought 
about, in the main, by the large absorption of the wanted wave under the chosen 
conditions, and that the increase of the collision frequency of the electrons with 
the change of frequency of the disturbing station is relatively unimportant. 
Translating Bailey’s nomenclature to ours, the change in (7,—#7,) is unim- 
portant, but the large value of the cross-modulation is brought about by a greatly 
increased value of «(pw). 


$9. SOME FURTHER POINTS 
(i) Gradients of Ionization Density 


From a consideration of the difference between the calculated and observed 
value of cross-modulation in the cases where the two regions of absorption do not 
completely overlap it is possible to make an estimate of the gradient of the ioni- 
zation in the lower part of region E at night. For instance, in the case of Droitwich 
200 ke/s. disturbing a transmission from the Decca Navigation station on 70 ke/s.. 
the low value of cross-modulation can be explained on the assumption of ar 
ionization density which increases by a factor of 1-7 per kilometre height. 

Ottringham 167kc/s. disturbing Westerglen 767 kc/s. gives 6°/ cross-modu-— 
lation, whereas Ottringham 1,112 kc/s, disturbing the same station gives only 2%. 
‘The geometrical arrangement and the frequencies of the second pair of stations 
are such that the level of reflection of the 767 kc/s. wave is below that of the 1,122 ke/s. 
disturbing wave. Again this gives sufficient data for the calculation of an 
ionization gradient and the factor of increase is here found to be 1-6 per kilometre. 

The calculations are approximate and the good agreement between the two 
cases must be regarded as fortuitous, but the order magnitude may be assumed to. 
be correct. For purpose of comparison, it may be pointed out that this gradient 
is considerably steeper than the maximum gradient of a Chapman region with a 
scale height of 10 kilometres. This figure supports the impressions gained from 
other experiments that, at night, the lower part of region E is sharply bounded. 
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(ii) Position of Region of Cross-Modulation 

Reference once again to equation (1) leads to the conclusion that the level at 
hich cross-modulation occurs is controlled by the absorption of energy from the 
fdisturbing wave alone. It is to be expected that the level of maximum cross- 
tmodulation will be that level at which the power extracted from the disturbing 
ave per unit distance isa maximum. At any level the power absorbed from the 
ave will depend upon the coefficient of absorption and upon the power entering 
tthe level. If the gradient of ionization density in region E at night is as great as we 
thave calculated, then the coefficient of absorption must increase rapidly with height. 
jlt might reasonably be expected, therefore, that the region of maximum power loss 
from the disturbing wave will be somewhere near, though not necessarily at, its 
level of reflection. It is further to be expected that, provided the geometrical con- 
‘ditions are the same, the heights of the regions of cross-modulation for different 
»wanted stations acted upon by the same disturbing station will be identical. 

The first point was checked in experiments in which Ottringham 167 kc/s. 
‘disturbed a wave from Stagshaw 1,050kc/s. In these experiments all the necessary 
heights are known, and using our calculated value of the ionization gradient, 
“it can be shown that the level of cross-modulation is about 2 kilometres below the 
Slevel of reflection of the wanted wave. ‘This is in reasonable agreement with 


The second point was verified by experiments carried out with two wanted 
c ransmitters at Westerglen, radiating on 767 kc/s. and 1,149 ke/s.; the disturbing 
station was Ottringham in both cases. ‘The average values of Gv/27 measured on 
the two wanted stations were 671 and 668, showing that the level at which cross- 


(iii) Cross- Modulation with High-power Disturbing Station 


According to the full theory of cross-modulation (Ratcliffe and Shaw 1948) the 
depth of the cross-modulation should be directly proportional to the power of the 
Pdisturbing station, and to its depth of modulation, provided that the increase of 
ithe energy of the electrons due to absorption is not comparable with their mean 
wthermalenergy. In other words, proportionality to power and depth of modulation 
liof the disturbing station holds so long as P/Gv< Q,. 

; Experiments in which a wave from Westerglen 767 kc/s. was disturbed by 
Behan 167 ke/s. radiating 173 kw., 346 kw. and 520 kw., showed that even at 
‘the highest power used, the depth of the transferred modulation was still 
) oroportional to the power of the disturbing station. From this we may deduce 
that the increase in electron temperature due to absorption is certainly not greater 
» chan 10°% of the mean thermal value, even when 520 kw. of power are radiated from 
| the disturbing station. 

) Other experiments, in which the depth of modulation of the disturbing station 
was varied gave, as expected, a completely linear relation between the transferred 
' nodulation and the depth of modulation of the disturbing station, 


i (iv) The Energy Loss Factor 

| Ratcliffe and Shaw in I used the value G=2-6 x 10-*. In this paper the value 
»ased, 1-3 x 10-3, is the result of a more recent determination of the quantity by 
+ Huxley and Zaazou (1949). Some of the experiments which have been carried 
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out are capable of being used to give some idea of the magnitude of G in the 
ionosphere. ‘This isimportant for the following reason : Laboratory experiments 
indicate that G is a function of Z/p, where Z is the accelerating field and p the pres- 
sure of the gas. Although the actual value of Z/p which obtains in the ionosphere 
for the disturbing powers used in these experiments is approximately equal to that 
used in the laboratory determination, the actual values of Z and p separately are far 
larger in the laboratory determination than in the ionosphere. It is important 
to know whether the assumption that G is a function of Z/p holds over the rather 
large range of values of Z and p involved. 

If cross-modulation takes place below the level where v = (py + pz), equation (3) 
shows that there will be a change of phase of 180° compared with the more normal 
conditions. If we take the greatest value of Gy/27 measured, for which the normal 
phase relationship applies, we know then that the value of v deduced from this 
maximum value of Gy/27 must be less than (Pw+p,). If we assume this upper 
limit for the value of v we can deduce a lower limit for the value of G. For the case 
of Droitwich disturbing Lisnagarvey we find G>0-9 x 10-3 and for Ottringham 
disturbing Westerglen, G>1-2 x 107°. ‘ 

The experiments with the high power disturbing station can also give a lower 
limit to G if the temperature of the ionosphere is assumed. ‘Taking the tempera- 
ature to be 300°K. we find G>0-6 x 10-°. 

None of our experiments gives an upper limit to G, but the good agreement 
between the calculated and observed magnitudes of the cross-modulation using 
G=1-3x10-2 1s reassuring, 

The above results show that it appears safe to assume that the value of G 
determined for air in the laboratory holds also for the ionosphere. 
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APPENDIX 


Since the experiments described in this paper appear to provide the oppor- 
tunity for measuring the collision frequency at a known height, it is important to 
discuss in more detail the accuracy to which that height is known. For this purpose 
consider Figure 8 in which B represents the region of cross-modulation and in 
which the actual and triangulated trajectories of the wanted and disturbing waves 
areshown. ‘To simplify the problem we shall neglect the earth’s magnetic field. 

As is well known, it can be shown (Breit and Tuve 1926, Appleton 1928) that 
for any distribution of electron density, the time for a group (pulse) to traverse a 
* certain horizontal component of the actual path is equal to the time required to 
{ cover the same horizontal distance by travelling along the triangulated path with 
+ the free-space velocity. As a special case, the total time from the transmitter to 

the receiver is the same for the actual path traversed with the group velocity and 
) the triangulated path traversed with the free-space velocity.* 


w? Figure 8. Diagram to illustrate the difference in the heights of the true and deduced regions of 
cross-modulation. 


. |! Now in the cross-modulation experiment mentioned in §3 we determine the 


J 


a time t=t/(DLB) (fp)+?t/(BCR) (fw) where ¢’(DLB) (fp) and t’(BCR) (fy) 
‘(| represent the times for groups, of frequencies fp and fy, to traverse the paths DLB 
» and BCR respectively. But by the theorem just quoted 


{(DLB) (fp)=— and (BCR) (fy) = 2 


+ so that the time ¢ corresponds to the passage of the wave over the distance 


© (DH+GR) with the free space velocity, 


i From the experiment, the triangulated path WER is known and on the simplest 


x theory the apparent region of cross- -modulation is located at B’ by a geometrical 


ut construction which makes 

| DB’ +B’R=ct=DH+GR. 
tg * Breit and Tuve and Appleton prove the theorem only for the total path but the extension 
© mentioned above is implicit in their work. 
B-2 
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In the case illustrated the actual region B of cross-modulation is below the 
calculated region B’. 

Before we can decide how far B is below B’ weneed some knowledge of the actual 
trajectory of the wave, which depends upon the distribution of electron density. 
It will be sufficient for our purpose to assume that the electron gradient is a linear 
function of the height above the lower edge AYC of the ionosphere. Appleton 
has shown that with this assumption, the curved parts of the trajectories are 
parabolae and that 

COS" 1g 


EX = XY y= 


a 


where i, is the angle of incidence, and « is a constant which determines the rate of 
increase of electron density through the formula p?=1—«y, y being the height 
above the reference level AYC. 

The calculation of the difference y, of the heights B and B’ is made difficult by 
the fact that B is not necessarily vertically below B’. In order to get some idea of 
the uncertainty we may notice that if the cross-modulation takes place near the 
level AYC, the apparent and true heights will be practically coincident, whereas if 
cross-modulation takes place at the top of the trajectory of the wanted wave, the 
maximum uncertainty which can exist will be equal to the maximum difference in 
the heights of the true and triangulated paths. We shall therefore take this 
difference (y,) as giving an approximate maximum value of the uncertainty. 

We attempt to estimate y, as follows. For waves incident vertically upon the 
ionosphere the equivalent height of reflection is at a distance 2/« above AYC, and 
for waves incident at an angle 7p, it is at a distance (2 cos*p)/a above AYC. ‘The 
difference in equivalent heights measured at vertical and oblique incidence is 
therefore (sin?z))/x. Now in the experiments which were carried out with 
Stagshaw, the equivalent height of reflection for the particular angle of incidence, 
60°, was found to be 96 kilometres. ‘The equivalent height at vertical incidence 
was not actually measured at the time of the experiments, but we may assume, 
probably without serious error, that it was about 110 kilometres. Hence we have 
14 =(2 sin? 60)/« whence « = 4 approximately and using this value of « we then 
obtain y, = cos? z/a =10 x }=2-5 kilometres. 

The uncertainty in the height measurement is seen to be not greater than 
2:5 kilometres. 
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ABSTRACT. Three problems of meteorological interest concerning the formation of radio 
ducts are treated. The first concerns the rate of growth of ducts, with special reference to 
ducts over the sea; it is shown that ducts form quickly at first and then more slowly as the 
maximum duct width is approached. The second problem concerns the formation of ducts 
under conditions of nocturnal cooling; it is shown that there is a limit to the height which the 
duct may attain. The third problem concerns the width of a radio duct necessary to enclose 
the track width of a mode associated with a propagated signal of wavelength A assuming a 
power law of modified refractive index with height; it is shown that the necessary width is 
») proportional to A?/3, and the factor of proportionality, which depends upon the index in the 
assumed power law, is calculated. 

The three problems have a common meteorological basis, for each is concerned with 
turbulent motion near the earth’s surface. The law of turbulent diffusion assumed here is 
that the coefficient of eddy diffusion is proportional to a power of the height, following the 
theory developed by O. G. Sutton and others. 


NEN OD Cet OlN) 


HE following symbols are used: x, 2, ¢, horizontal, vertical and time co- 
ordinates; u, v, horizontal and vertical velocities; K, eddy diffusivity in the 
atmosphere, assumed given by K=az'”, a and m being parameters; V,a 
meteorological variable, either 6 the potential temperature in the atmosphere, or 
H gq the water vapour content of the atmosphere; «, the critical value of the gradient 
of V necessary for duct formation, equal to 5-4° F/100 itxorl0-*se/cm. tor V =0- 
yand —0-54gm. water vapour/kg. air/100 ft. for V=q; p, density of the air; 
©, specific heat of air at constant pressure; 7, temperature in the earth; K,, 
+ conductivity in the upper layers of the earth; p,, density in the upper layers of the 
earth; c,, specific heat in the upper layers of the earth; R, rate of loss of heat by 
» radiation at the earth’s surface; A, wavelength of a propagated signal; jy, ratio of 
' modified refractive index of the atmosphere to that at the surface. 
A radio duct is formed in the atmosphere when the variation with height of 
' the modified refractive index obeys certain conditions; in particular there is a 
surface based duct if the modified refractive index decreases from the surface 
_ upwards to a minimum and then increases again, and the top of the duct is at the 
same height as the minimum of modified refractive index. ‘The variation with 
height of modified refractive index may be expressed in terms of gradients of 
temperature and humidity in the atmosphere; if there is a layer above the earth’s 
surface in which the gradient of potential temperature is positive and exceeds 
5-4°r/110 ft. or 10-3 °c/cm., then this layer will form a radio duct, and the top of 
the duct is the height at which the gradient falls below this critical value. This 
criterion assumes that the humidity remains constant in the vertical, and the 
duct is termed a ‘dry’ duct. Similarly if the gradient of segue humidity is 
: negative and is less than —0-54gm. water vapour/kg. air/100 ft. a ‘ wet’ duct is 
formed, it being assumed that the potential temperature remains constant with 


height. 
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The temperature and humidity gradients near the earth’s surface are controlled 
by turbulent exchange of heat and water vapour. ‘The three problems treated 
below are aspects of the same fundamental meteorological problem, the effect of 
a boundary upon the air mass above it, although there is dissimilarity in the methods 
of attack. The theory of turbulence developed by Sutton (1949) and others in 
connection with the diffusion of smoke near the earth’s surface is used and this 
implies taking the coefficient of eddy diffusivity to be a power law of the height 
above the surface, 1.e. 

Keagee 


where a and m are assumed constant in any one meteorological situation. 

The three problems arise in practical forecasting for radio purposes and the 
object is to give numerical estimates of important quantities arising from the 
problems. 


§2. THE RATE OF GROWTH OF RADIO DUCTS OVER THE SEA 


When relatively warm dry air flows over the sea, or any cool moist surface, 
there is turbulent exchange of heat and moisture and a radio duct is formed. 
No estimate seems to have been made of the time interval for the duct to reach a 
given height, although the problem is of importance to the meteorologist fore- 
casting for radio purposes. 

The problem considered here is associated with the modification of an air mass 
in juxtaposition with a surface of different physical properties. The conditions 
assumed are: (a) At the lower boundary, z =0, V remains constant; the sea surface 
is a boundary satisfying this condition approximately both for V =6 and for 
V=q. (6) Above the turbulent layer the boundary condition is the same; V 
remains constant. ‘This condition is satisfied in the case of a well mixed homo- 
geneous air mass passing over, and being modified in the lower layers by, a surface 
of different physical properties. 


(1) The Differential Equation 
Assuming that turbulence is horizontally isotropic and that the diffusivity 
follows a power law, then the governing differential equation is 


Ay nd) eve ava 
az = — +w— 
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where x is measured in the direction of horizontal velocity u. Except in special 
cases this equation is difficult to deal with; it is noted, however, that if w=0, the 
term u(dV/dx) is such as to increase the rate of transfer of V, compared with that 
when u=0, by continuously replacing the air at any point by less modified air. 
‘Thus in the case considered the solution of 


0 OV OV 
1—m a 
ae (: =: pice eee pe (1) 


will give an upper limit to the time required to form a duct of given height. The 
actual time required may well be much less than that calculated from the above 
equation, but it cannot be greater. 

If the values of V at the boundaries are: 


V=V, at 2=0) V=V, as 2>0, 
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hen the appropriate solution of i is 
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(11) The Maximum Duct Height 
If « is the critical value of the gradient of V at the top of the duct then 
eras tal 
D{m/(m + 1)} 
is simply an equation in the duct height z and the time ¢. Noting that e-Ypm/+) 


has a maximum value at %=m/(m+ 1), then there is an upper limit to the height 
: of the duct z given by 
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eee wo 
max = Tin] m+ Iya 
This is a somewhat different maximum height from that obtained by considering 


the steady state, but the two maxima only differ numerically by about 10% in 
extreme cases. Figure 1 shows the maximum duct width for various values of m. 


2 (m+ Lee ie”. p= 
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Time in seconds for duct formation 


Ory Duct Height in feet 


a 5 PAR no a : Profile - index 

Figure 2. ‘Time in seconds for forma- 
tion of maximum duct for values of 
m, assuming azl—m=3x10* at 
z=100 metres referring to 1°F. 
difference between boundaries. For 
T °F. difference, multiply time by 
Ti+m, For a wet duct, difference 
x gm. water vap/kg. air between 
boundaries, multiply time by 
(10x)m+1, 


. Figure 1. Maximum duct height for 
values of m. For 1° F. difference 
between upper and lower bound- 
aries. For difference T° F. multiply 
d by T. For wet duct multiply d 
| by 10 for each 1 gm. water vap/kg. 
air difference between upper and 
lower boundaries. 
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The time taken to reach the maximum is given by the equations ? | 
m+. 
joes ee (5) 
m+1  (m+1)at 


Some determination of the values of a and m is necessary in order to give some 
idea of the magnitude of the upper limit to the time taken to reach the maximum 
duct width. Sir G. I. Taylor found values of K of the order 3 x 10? in c.c.s. units 
for moderate winds over the sea (Brunt 1939, Chap. XII). With this value of | 
K, Figure 2 shows the dependence on m of the time taken to reach the maximum 
duct width. | 

(iii) Ducts less than Maximum | 


From Figure 1 it can be determined whether a duct of given height will or | 
will not be formed. For duct heights less than the maximum the time ¢ taken to — 
reach any specified height z) is determined from 


az_l {m|/(m+ 1)} = oS a Be 

(m+1)(V,—Vo)’ ~ (m+1)Pat” 

As an illustration take m=4 and let the temperature difference between the upper 
and lower boundaries be 9°F. Then from Figure 1 the maximum duct width 
is about 50 ft. and from Figure 2 it will be formed in less than about $ hour. 
If the humidity difference between the boundaries is 2gm. water vapour per 
kg. air, the maximum duct height is about 80 ft. and the upper limit to the time 
taken for the duct to grow to this height is about 1} hours. For a 9° Fr. difference 
between the boundaries the upper limit for the time for the duct to rise to a 
height of 25 ft. is calculated as follows. 


Setting V,—Vy=5° c (=9" F)),. w= 10-2 cjem., 25750 cm = 2500, 
750 x 10-3 x T(4) 
™ So 


EP yfgrl (m+) _ 


we find Eye = = 02267, 


whence # = 1-4 and 

G Bee a <: 7503/2 

~— (m+1Pab (3/2)? x 30x 1-4 
The half maximum duct width is thus reached very quickly. 


t = 220 sec. 


(iv) Conclusion 
It may be concluded in general that the duct forms rapidly at first and approaches 
its maximum height more slowly. In the case of a warm dry airstream passing 
over the sea, where the air speed is approximated from synoptic observations, the 
duct probably forms sufficiently quickly for the time lag to be neglected in 
comparison with the errors in speed estimation, since the duct must extend laterally 
for a great distance to be effective in a radio sense. 


§3. THE FORMATION OF RADIO DUCTS UNDER CONDITIONS OF 
NOCTURNAL COOLING 


At night heat is radiated from the earth’s surface and the heat loss is greatest 
when the sky is clear; the rate at which heat is lost is sensibly constant (Brunt 1939) 
and is derived both from the top layers of the earth and from the lower layers of 
the atmosphere. An inversion of temperature forms in the atmosphere and if 
radiation continues for a sufficient length of time there is a layer in the atmosphere 
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which the gradient of potential temperature exceeds the critical value for duct 
ormation and a radio duct is formed. It is assumed that the humidity remains. 
onstant in the vertical; a remark about this assumption will be made later. 
Conduction of heat in the earth’s surface is assumed to follow the classical law. 
n theatmosphere the diffusion coefficient is again assumed to follow the power law. 


(1) Mathematical Equations 


The differential equations governing the system are: 


0 ks 00 alo 
= (# 5) =F 220), f= Macnee (6) 

otod Deas db 
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here z is measured from the earth’s surface. 
The boundary conditions are: 

Geet eee mee tl) eel A ae ie ie a cas (8) 
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here 6 is the change in potential temperature of the atmosphere from the value 
tt=0, T is the change in temperature in the earth’s upper layers from the value 
: t ¢=0, and the other symbols are as defined above. ‘These equations have been 
wziven by Jaeger (1945). 

_ It should be pointed out that at some time shortly after dusk the temperature 
jn the first hundred metres of the air becomes almost constant with changing 
reight, and so a natural time zero is provided (Johnson and Heywood 1938). 
Since the gradients in the atmosphere with which we are concerned in this problem 
are large this time zero can be taken as a time at which the potential temperature 
is almost independent of height. 


(ii) The Solution of the Equations 
The solution of equations (6) to (9) is effected by the operational method, in 


which the transform F(z, p) of F(z, t) is defined by: 


| F(z,p)=p { "exp (—pt)F(z, t) dt. 
0 
We find (Knighting 1950) 
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where K,, is the modified Bessel function. 
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(iii) Application of the Solution to the Formation of Radio Ducts 


For a particular value « of 00/0z equation (11) is simply an equation between 
z and p, or alternatively between z and ¢, and we have 


m i-m EU ciggepie lad 
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For the purpose of calculating duct heights « is 10~° “c/em. To interpret (12). 
split the right-hand side into two factors, the first of these being | 


g UK mews ); | 


and apply Duhamel’s theorem. The transform of this factor is always positive ) 
and less than 24/"+D-1[P{] /(m+1)}. Since m is known to be less than unity, 
an upper limit to the value of z given by equation (12) is found by putting 


EMD Kn gE) SO DA (m+ 1}. 


Making this approximation and carrying out some simplifications the upper limit 
of the value of z, say 2, is given by: 


R(m+ 1m tp AD (1 (m+ 1)} 
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apc,(m+ 1)2"+ OP + 1/(m+1)} 


where B = ~ al p.c,K Ye {m|(m+ 1} ee Peer) 


m=1/7. 

here values of the constants near the earth’s surface are, in C.G.S. 
units, R=2:1x 10°, p;=1-6, c,=0-2, K,=4-7 x 10-8, ¢,=0-24, p=1-27 x 10-3 
The value given for R is applicable to Southern England and the surface values 
for fine dry loam. Values of m and a vary with the meteorological conditions. 


Taking m=1/7, and a=1-12 so that the diffusivity at 100 metres is 3 x 108, equation 
(13) becomes 


ga, _ 21 x (8/780 (7/8) _ 
16p>8(1 + 2-5 x 10-3p-38) ° 


Expanding in negative powers of p and interpreting term by term, the series 
in ¢ is 


Zo 


{3/8 2-55 10-3 ¢6/8 
PUesS Ao ss Pee 
This series is rapidly convergent for t = 3,600 sec. and with this value of t only the 
first term need be retained; it is found that z29~11m. or 35ft. Hence the 


maximum height the duct can attain after radiative cooling has been operative 
for one hour is about 35 ft. 


i= Nis. 
The method of the preceding section is useful when f is sufficiently small 


and the resulting series is rapidly convergent. Cowling and White (1941) 
analysed some experiments carried out by Johnson and Heywood (1938) and gave 
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alues for the diffusivity at various times through the night. In the mean the 
alues may be represented by the law az!-™=25z!-"3_ Using these values for 
and m equation (13) reduces to 


15 15 
2/3 — 
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O interpret this equation a contour integration is carried out and it is found that 
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» Taking a value of t=3,600sec. and performing the integration it is found that 
“%)~ 10m. or 30 ft. Hence again we find that the maximum height the duct 
‘an attain after radiative cooling has been operative for one hour is about 30 ft, 


. (iv) The Maximum Possible Duct Height 


i An important point to note is that the term in 1/f in equation (16) is independent 

fof the approximation made for €/"+)K,,,,,(€), for it arises from integration 
»around a small circle at the origin where the approximation is exact. Hence 
: che maximum possible duct height is given in the case treated in sub-section (iii) 
Soy 223=15/8 or z~40m. or 120 ft. From the definition of 8 the maximum 
“possible duct height is given by z2--”=R/aapc,. ‘The maximum is large for small 
wvalues of a as might be expected since small values of a indicate that diffusion is 
wacting slowly. A time calculation shows that for small values of a the maximum 
pis approached after a very long time and the physical conditions would have 
wchanged, making equations (6) to (9) inapplicable. 


(v) Calculation of the Duct Height from the Ground Temperatures 
We can also use the fall in temperature at the ground to calculate the height 
‘ofthe duct. Consider again the case treated in sub-section (iii). From equation 
(10) the fall in temperature at the surface is given by 
nae —R 
9° eK BP 


—R 1 1 
~ Bp ce, Ky¥? {a - seep ee (17) 
“The height of the duct is given, from equation (11), by 
ie, pe a (18) 


Picea )(Peetp)e ~~ ts 
d emstituting from (18) into (17) and interpreting the equation so formed 
Rt4 2 CE SG) 

Bete Lays) EGE 
“where AT denotes the fall in temperature at the ground in the time interval (0, ¢). 
‘’With the values of the constants used above equation (19) reduces to 

Pile omeeesN Te a eG (20) 

| This formula refers to summer conditions in Southern England, since the para- 
‘meters used correspond to those observed in summer. At the end of a June 
bnight t~ 7 x 3,600sec., AT ~ 10°c. which gives z ~ 30 ft. Equation (20) does 
‘not hold everywhere, the constants varying both with place and time of year, but 
lit should be possible to obtain a similar formula for many parts of the world. 
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It is important to note that the duct height attained after radiative cooling 
has been operative for seven hours is about the same as the maximum duct height | 
attainable after one hour. We may infer that ducts grow rapidly at first and then | 
more slowly approach their final height. 


(vi) The Effect of Water Vapour 


The calculations of the dry duct heights have been made on the assumption 
that the vertical gradient of water vapour remains constant. Measurements of 
temperature and water vapour have been made at Rye at various heights in the 
atmosphere. The analysed results, which have not yet been published, show 
that during the night there is generally a positive gradient of water vapour. The 
effect of this inversion of water vapour is to inhibit duct formation and hence the. 
duct heights given in this paper will generally be overestimated. A quantitative 
attack on the problem has not yet been made. 


§4. THE DUCT WIDTHS ASSOCIATED WFEH A POWER LAW 
PROFILE OF MODIFIED REFRACTIVE INDEX 


Assuming that the diffusion coefficient in the atmosphere follows the power 


law an approximation to the variation of modified refractive index with height has. 
been given (Booker and Walkinshaw 1946) as 


ees jn 
peal 2c(2 D -) ee (21) 


where y 1s the ratio of the modified refractive index at height z to that at the 
surface, m is the index 1n the power law, and C, D are constants. Since p is very 
approximately unity 


d 
hea G as. => 00, 


dz 


It is known (Booker and Walkinshaw 1946) that du/dz—4K,/5 as so, where 
K, is the curvature of the earth, and hence C= —4K,/5. From equation (21) 
uw has a minimum at z=D and hence a surface-based duct exists to the height D, 
The minimum value of uw is given by 

peal =20 DES iin) eee (22) 
Although the value of given by equation (21) depends essentially upon a single 


parameter, a modified refractive index profile of this form can often be fitted to 
experimental data. 


(i) The Eigenvalue Problem 


Following the electromagnetic theory of mode propagation as presented by 
Booker and Walkinshaw the eigenvalue problem to be solved is 


| * dae=TiaN, | eee (23) 
0 
where mis a positive integer, w,, is the track width of the mth mode, A the propagated 
wavelength, and q is defined by 
Q? p= C088 ek eee (24) 


y being the angle of propagation with the horizon at the transmitter. w 
further related to the variables being a solution of the equation 


q=0. REA FL 
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quation (23) is regarded by the authors as an eigenvalue equation in y and leads 

o solutions of the problems that they formulate. Equation (23) can be regarded 
s an eigenvalue equation in w,, and so treated leads to a simple result concerning 
he duct width necessary to enclose the track widths of the first 2 modes. 


(it) The Duct Width enclosing n Track Widths 


From equation (23) if w,, is identified with D then the problem of determining 
he duct width enclosing the first 7 track widths is formulated. For the trapping 
f electromagnetic energy in a duct in the atmosphere the angle of propagation is 
mall, being certainly less than one degree. From equations (21) and (24) 


g=y—2C (2-De =| es ey Dae (26) 
Now from equation (25) g=0 at s=D, and hence 


#=2C(D- Dist D™ =). me (27) 
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te Figure 4. Duct width for wavelength 
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‘The eigenvalue equation, equation (23), is now 
; D D ym 1/2 
I, | 2c(p- — —2+Dim =) | iE Sie hey: (28) 


Setting z= Dt in this equation 
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The integral in equation (29) is purely a function of m, and writing 
il | pm 1/2 
= = —— al (A? yk IE 30 
1=I(m) Re 14t “) dt, (30) 
1(y — 1)) 2/3 
then D -{ CaCO, a re (1) 


The values of /(m) have been plotted against m in Figure 3. Equation (31) gives 
in general the duct width required to enclose the first 7 track widths. In practice 
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i 


the first mode is often the most important, and setting »=1 in equation (31)) 
and inserting the value of C( = —4K,/5) 


X) 2/3 
D=12'5 {it fi. --. = (32) 


where A is measured in centimetres. Figure 4 gives D in terms of m for A=1 cm. | 
It is easily seen that for m=4 the duct width necessary to enclose the track Wee | 


of the first mode, and hence to suppress leakage, is about 50 ft. for a wavelength 
of 3cm. and about 120 ft. for a wavelength of 10cm. Equation (32) thus provides. 
an estimate of the duct width necessary to enclose the first track width of a 
propagated signal of wavelength A, with the assumed relation between modified | 
refractive index and height. | 
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ABSTRACT. The design and operation of Rosenblum-type spark counters having 
plane, concave and convex cathodes is described, and discussed with reference to the 
counting characteristics, useful life and background counting rate. The absolute efficiency 
of the counter in air has been determined for «-particles, knock-on protons, and f-particles, 
giving values of about 90%, 5% and 0:0002%, respectively. A marked dependence of 
efficiency on specific ionization and air pressure is evident. The counter is directional in 
all its forms. The effect of employing gases other than air is also investigated and a 
qualitative description of the counting characteristics is given. Alpha-particle ranges may 


be determined using either a single counter or two counters used in anticoincidence as is. 
described in the text. 


§1. INTRODUCTION 
HE original paper describing an arrangement whereby ionizing particles 
| or quanta, passing through the sensitive region of an electrode assembly, 
are registered by sparks passing across the electrode gap, is that due to 
Greinacher (1934 a). ‘The name ‘ Spark Counter’ was proposed, and subsequent 
papers by the same author described in general terms the properties of various 
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ypes of electrodes and various geometrical arrangements (Greinacher 1934b, 

935 a,b, c, 1938 a,b). Later more detailed investigations were reported by other 

orkers (Stuber 1939, Jeanguenin 1944, Frey 1946, Rytz 1949). In all these 
ounters the electrode separation was small and one electrode, usually the anode, 
ad a small radius of curvature. The counters described by Greinacher and 
tuber were sensitive to gamma rays and ultra-violet photons as well as to alpha- 

d beta-particles, but they had many important disadvantages: 

(i) No ‘ plateau’ region in the counting characteristics appeared for any of 
he radiations mentioned, unless very particular electrode systems were used 
Frey 1946). 

(ii) The sensitive volume increased with applied voltage, and was different 
for the different ionizing radiations. Similarly, the threshold voltage was different 
for radiations of different ionizing power. 

(111) Operation was SLES sensitive to the surface conditions of te metal 
lectrodes. 

(iv) A high ‘ background’ counting rate was usually observed and the 
esolving time was long (about 0-03 sec.). 

_ It was shown by Chang and Rosenblum (1945) that the principal disadvantages 
ere removed if a flat metal plate acted as the cathode and the anode was a fine 
ire stretched in front of the plate at a distance of about 15mm. The counter 
perated in air at atmospheric pressure. A quench circuit was used to enable 

‘the counter to recover after the passage of a spark. Applied voltages were of the 
rder of 3 ky. and the counter was sensitive to alpha-particles only. Rosenblum 

‘found the background to be small (about one count in three minutes for a wire 
cm. long). -When the properties of this type of spark counter were still 
npublished, an examination of the counter was undertaken by the present author 

with a view to its development, and a brief summary of part of this work has 

Jappeared elsewhere (Connor 1949). It is the purpose of this communication to 

“present a fuller account of these researches. Since the appearance of the 

‘preliminary account, a report of similar but independent work on some of these 

“topics has been published by Payne (1949). 


§2. APPARATUS 


In the course of the present work, it has been shown that the original wire- plate 
Ycounter of Rosenblum is a particular case of a more general type of counter in 
‘which the anode is a fine wire, and the cathode may be convex, planar or concave 
‘with respect to the anode. Throughout most of the work the anode—cathode 
“separation was about 1mm., and under these conditions the counting threshold 
‘voltage was of the order of 3kv. ‘The counters were operated in air at atmospheric 
» pressure in all but a few experiments which are mentioned separately below. In 
order to facilitate changes in the counter geometry, two electrode carriages were 
constructed, one of which is shown in Figure 1. The essential features are two 
»ebonite blocks on which the anode and cathode can be mounted. ‘The lower 
- block is fixed to the vertical brass plate and carries the cathode clamped in position, 
swhile the upper block, which carries the anode, can be made to move vertically 
jinto any desired position with respect to the lower block by the action of an 
,accurate screw. Side play is avoided by the careful milling of the slot in which 
,the upper block moves. The screw is provided with a divided head and acts 
against the controlling action of two springs so that backlash is avoided. ‘The 
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anode, a fine wire of diameter 0-05 mm., is stretched horizontally across the ends : 
of two glass pillars or other insulators, which protrude from the lower surface of | 
the upper ebonite block, the wire being pulled taut and secured at each side of the : 
block by friction between two nuts on a screw sunk into the side of the block. | 
Tungsten is chosen as the anode material on account of the tension to which the 
anode wire is subjected. The cathode, formed either of wires of larger diameter 
than the anode, plates or hollow cylinders (representing convex, planar or concave : 
geometries with respect to the anode), was secured by clamping to the other | 
ebonite block. In the first two cases, the anode and cathode were aligned by, 
screwing the electrodes together and adjusting them until the tungsten anode : 
wire lay evenly along the cathode in a symmetrical fashion. By means of the screw, 


SIDE VIEW FRONT VIEW 


Figure 1. 


the gap could then be set at any desired value, the parallelism between the 
electrodes being maintained by the carriage. In the case of the hollow cylinder 
arrangement, the anode wire passed centrally down the cylinder and a slit or 
window was provided to permit entry of «-particles from a direction normal to 
the axis of symmetry. It is clearly of importance to have the electrode separation 
uniform along the entire length of the anode wire, for the counter will tend to 
spark preferentially at points of minimum spacing. In order to avoid end and 
edge effects, a vertical slit was placed between the source of «-particles and the 
electrodes, thus restricting the particles to the central region of the counter wire. 
In the majority of these investigations, the source of «-particles was the active 
deposit of thorium collected, by recoil from a dry emanating source of radio- 
thorium, on the surface of a small 8B.a. screw head which had been filed flat and 
polished. 

The cathode was connected directly to the negative output of a stabilized 
power supply, having an output voltage variable up to 4,000 volts, and the anode 
was connected to earth via a resistance capacity potentiometer. This served two 
purposes: firstly, it acted as a quench circuit to extinguish the discharge, and 
secondly, it provided a known fraction of the voltage surge for application to a 
scaler. ‘Ihe potentiometer circuit was connected to the anode with the shortest 
possible length of lead. ‘This precaution was necessary, since, if the connection 
was more than a few inches long, a damped oscillation appeared on the pulse after 
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he leading edge and in certain arrangements lasted for several microseconds and 
as of such amplitude as to produce ‘ double counts’ in a sensitive scaler. Visual 
xamination on a calibrated cathode-ray oscillograph of pulses, obtained under 
ood conditions, revealed an exceedingly sharp leading edge, the rise time being 
stimated as less than 10-7 sec. ‘The recovery time was of course largely controlled 
by the time constant of the quench circuit. A time constant of 50 microseconds 
was found to be suitable and the smallest value consistent with stable operation 
was shown to be 20 microseconds for the counters here described. 

Measurement of the total pulse height indicated that on discharge all but a 
few per cent of the total anode—cathode potential difference is developed across 
the quench circuit. Clearly, no amplifier is needed when pulses of such magnitude 
are available and, as Figure 2 shows, the quench circuit used supplied about 1/51 
of the total surge directly to the scaler. 
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$3, THE COUNTING CHARACTERISTICS 


From the geometrical details given above, it will be seen that for counters 
- having wire and plate cathodes, «particles can approach the counter gap from two 
} principal directions, both normal to the anode, the one parallel (‘ End-on’ 
, direction) and the other perpendicular (‘ Broadside’ direction) to the plane 
- containing the axes of the electrodes. For counters with hollow cylindrical 
, cathodes, in which the anode passes centrally down the axis of the cylinder, only 
one principal direction can be specified, corresponding to the ‘ End-on’ direction 
of the other arrangements, in which the «-particles pass through a slit in the 
cathode into the sensitive volume in the direction of the normal to the axis of 
symmetry. Figure 3 shows the Broadside characteristic obtained with an anode 
consisting of a tungsten wire 0:05 mm. in diameter and a cathode of copper wire 
0-9 mm. in diameter, the electrode spacing being 0:75mm. Figure 4 shows an 
End-on characteristic for the same counter. In these figures the counting rate 
has been corrected for source decay and normalized to a constant source—counter 
geometry. The radioactive source was placed 7:9cm. from the anode wire in 
each case, the atmospheric pressure being 76-7cm.Hg. In these experiments 
only the ThC’ «-particles could reach the counter. At this stage it is appropriate 
to mention that at no time was it found possible to record, with counters of the 
PROC. PHYS. SOC. LXIV, I—B c 
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types here used, any appreciable effect due to the B- or y-rays from sources up 


to 10 mc. in activity. Furthermore, «-particles having ranges too short to reach 


the counter gap produce no effect, although they are producing intense ionization _ 


only a few millimetres away from the gap. The Broadside curve is of a type not 
found with any other counter. After a rapid rise from the counting threshold, 
the curve passes through a sharp maximum marked B in Figure 3, subsequently 
falling to a minimum value. The counting rate therafter increases to a flat 
maximum (point A on Figure 3) as the applied voltage is further increased. 


About A, the curve exhibits a ‘usable’ region of some 300 to 400 volts in which | 
the dependence of counting rate on voltage is small. In general, it is found that _ 


for a given anode the voltage at which the flat maximum occurs is greater, the larger _ 


the cathode diameter and the larger the electrode spacing. Also the flat maximum 
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counting rate increases with spacing. It should be mentioned that the abscissae 
of all counting curves shown in this paper represent the voltage as read on the 
output meter of the power supply. For reasons which will appear later, this. 
voltage is not that appearing across the electrode gap. 

The End-on characteristic is quite different from the other. After a more 


gradual rise, the counting-rate curve passes through a small peak and thereafter 


enters a long plateau region which can extend for 1,000 volts or more. Figure 5 
shows the End-on curve for a Rosenblum wire-plate arrangement of spacing 1 mm. 
and Figure 6 the corresponding curve for a hollow cathode counter, the cathode 
being of brass and of internal diameter 2:07 mm. The same anode arrangement 
was used throughout. With wire cathodes, the plateau is not strictly flat but 
falls by 1 or 2°%, over the entire plateau region. In this way, Figure 4 is seen 
to resemble the curves for the other two geometries but with features less marked. 
Wire cathode counters have been investigated for cathode diameters ranging 
from 4-7 to 0:-4mm., using several different materials and with spacings ranging 
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rom 0:5 to 1:25mm. The graphs of Figures 3 to 6 inclusive refer to counters 
which have been suitably conditioned. This procedure consists in cleaning 
the cathode with metal polish, then washing with chloroform and alcohol, and 
nally, allowing the counter to operate for some 5,000 to 10,000 counts, which 
recaution, according to Stuber and others, removes the small irregularities 
rom the electrodes. Should this precaution be omitted, the counting rate is 
ound to be a few per cent too low. Counters with tungsten anodes of 0-1 mm. 
diameter exhibit characteristics similar to those using 0-5mm. diameter wires, 
but higher voltages are required. No stable characteristic could be obtained when 

he polarity of the counter was reversed, that is, with a fine wire cathode and an 
sanode having any of the geometrical forms mentioned above. 


§4. THE CHARACTERISTICS OF THE COUNTER IN THE 
QUIESCENT STATE 

If a microammeter is incorporated in series with the quench circuit, some 
information can be obtained as to the current passing between the electrodes in 
the quiescent, that is, non-sparking state. It is observed that, when a counter 
shaving any one of the three types of cathode described above operates in air at 
jatmospheric pressure in the absence of a source, a corona discharge occurs round 
‘the anode wire as the voltage applied i is raised above that corresponding to the 
‘counting threshold. With anode wires approximately 2 cm. in length, the corona 
lis first observed in a darkened room when the current registered by the meter is 
of the order of 5x 10-*amps. ‘This current rapidly increases with applied voltage 
until, with a spacing of the order of 1 mm. and 4,000 volts applied from the power 
‘supply, the quiescent current reaches a value between 100 and 200 pQ depending 
‘on the cathode geometry. To the eye the glow is not at all diffuse, appearing as a 
fine line of light somewhat thicker and more intense in the direction of the cathode, 
and extremely faint if not actually invisible on the side of the anode away from the 
‘cathode. ‘This visible corona extends about 0-01 mm. into the gap from the 
surface of the anode in the direction of the shortest path between anode and 
‘cathode, and the luminous boundary appears to the eye to be sharply defined. 
The glow possesses a sharp ‘ offset’ or extinguishing voltage and a higher but 
‘somewhat variable ‘ onset’ voltage, this being a common effect noticed by other 
observers in a positive point to plane corona. The ‘ offset ’ voltage has been shown 
to coincide with the voltage of the peak B of the Broadside, and also with that of 
the small peak of the End-on characteristic. It may be noted here that, because 
of the existence of the quiescent (corona) current, the voltage actually appearing 
across the electrode gap during any counting experiment 1s less than that read on 
the meter of the power supply as already briefly mentioned. 

It should also be mentioned that for counters having wire cathodes of copper, 
nickel silver and aluminium, of the same diameter and the same electrode spacing, 
the counting rate at the flat maximum of the characteristic was the same to within 
the experimental error, irradiation being under standard Broadside conditions 
throughout. 

When the quiescent discharge was examined on a cathode-ray oscillograph, 
“corona pulses’ were observed corresponding to those photographed by ‘T'richel 
(1939). The corona pulses were of particularly large amplitude at the voltage 
corresponding to the peak B of the Broadside characteristic. ‘The corona 
fluctuations at this and adjacent voltages were sufficiently large as to be audible. 
c-2 
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$5, ‘THE ERFECT (OF SPECIFIC IONIZATION 

The probability of recording the passage of an a-particle through the 
potentially sensitive region of the counter (or, alternatively, the extent of the truly 
sensitive region) has been shown to be a function of the specific ionization of 4 
a-particle as was to be expected from the almost complete insensitivity of the 
counter to 6-radiation. In Figure 7 the normalized counting rate-range cur | 
are given for the Broadside and End-on positions, and for 15°¢%] 
760mm. Hg. The term ‘ Range’ is here taken as being the distance of the source ; 
from the anode wire, thus it has not precisely the same significance in the two | 
cases, but it is interesting to note that both curves lie within Henderson's (1921) | 
ionization.curve. In each case the voltage corresponding to the flat maximum of | 
the Broadside counting curve was used. 
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Figure 7. 


In view of this effect of specific ionization, it is to be expected that the counting 
characteristic will show some variation if the residual range of the incident 
a-particles is altered. Using the Broadside position, it was found that for residual 
ranges greater than about 6 mm. (in standard air), using cylindrical wire cathodes, 
the sharp peak of the counting curve no longer occurred at the corona offset 
voltage, the counting threshold being some 20 volts or more above the latter. 
Under such conditions the rise in the counting rate from the threshold is 
particularly steep, and the fall, in passing through the sharp peak, is corres- 
pondingly rapid. ‘The normalized height of this peak decreases as the residual 
range is increased, and the flat maximum shifts to regions of higher voltage. The 
shift of the threshold is small, but the shift of the flat maximum can be as much as 
300 volts for a change in residual range from 0:5 to 25cm. 


$or THE USER UL LIP Ore RE +COWINTEER 
Five Broadside counting characteristics were taken with one counter in which 
more than four million sparks passed between the taking of the first and last 
curves, and no significant change in the characteristics was found. The length 
of time for which a counter of this type can be used is determined principally by 
the extent to which dust, fibres, etc. collect on the electrodes. If the counter is 
operated in a bell jar using filtered air, it does not require cleaning so frequently 
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s does the same counter used in open air, but if anode and cathode are cleaned 
ach time prior to use, the same electrodes can be used over and over again, until 
nally it is observed that the counter is discharging locally. Usually local 
ischarge arises from an irregularity in the electrodes, some pitting or erosion, 
ind both electrodes should be replaced when this occurs. 


§7. BACKGROUND 


The smallness of the background counting rate is one of the attractive features 
»f counters of the Rosenblum type. Under ‘ clean’ conditions, viz. freedom 
rom «-emitting radioactive contamination and dust, the background has been 
letermined on occasion as less than two counts per hour. With the counter 
incovered, the background is generally less than ten counts per hour, but naturally 
t varies with the dust content of the air. In a room with a wooden floor, vigorous 
rushing of the floor has been known to result in a temporary rise in rate to the 
tartling figure of five counts per minute when the counter was uncovered ! 


§8. COUNTING IN AIR UNDER REDUCED PRESSURE 


Counting characteristics for both Broadside and End-on positions have been 
ybtained for air pressures lower than 1 atmosphere. The general shape of the 
‘urve remains unaltered, but the counting threshold moves to lower and lower 
roltages as the pressure is reduced. ‘The entire curve also contracts so that the 
roltage range between threshold and continuous discharge becomes smaller. 
3roadside counting curves were obtained for a counter having a spacing of 
75 mm. with a cathode 0-9 mm. in diameter and a series of values was obtained 
or the normalized flat maximum counting rate as a function of pressure. The 
yosition of the source was adjusted so that the greatest normalized counting rate 
vas obtained in each case. The graph exhibiting these values is equivalent to an 
fficiency—pressure graph and is shown in Figure 8. The rate of variation of 
ficiency with pressure is seen to be small in the region of 76cm. Hg but not 
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Figure 8. 


egligible, but at half an atmosphere pressure, the rate of variation is relatively 
urge. It is therefore both convenient and expedient to operate the counter in 
ir at atmospheric pressure. 

This result has obvious correlations with that concerning the effect of residual 
ange and provides further confirmation for the view that the probability of spark 
xrmation depends on the amount of ionization produced by the individual 
-particles in a specified region of the counter. 
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§9. THE ABSOLUTE EFFICIENCY OF THE COMUNGER p | 

Since the efficiency of the counter was found to vary with the air pressure and — 
the residual range of the «-particle, it became important to measure the absolute _ 
efficiency of the counter under standard conditions, viz. in air at atmospheric 
pressure and with the source-counter distance such as to give the maximum 
normalized counting rate. This was done for the Broadside position, using a 
collimating slit to define the length of counter employed. ‘The number of | 
a-particles per second passing this slit and subsequently passing between the | 
counter electrodes was determined in a separate experiment with a thin windowed 
proportional counter of known aperture. ‘The ratio of actual counting rate to the | 
number of «-particles passing through the spark counter gap per second is taken | 
as a measure of the absolute efficiency (Table 1). 


Table 1. Spark Counter 


Cathode . Measured 
Spacing (mm.) 


diameter (mm.) Efficiency 
0-93 0-71 0-98-+.0-04 
0-93 0-71 1:04-+0-04 
0-93 0-78 0:98 40-06 
0-70 0-80 0:90 40-04 
0-70 0-61 0:96 +0-04 
0-70 0-65 0-84-+0-03 
0-70 0-74 0-90-+0-05 
0-70 0:87 0-83 40-06 


Unit efficiency can only be obtained if the whole width of the gap is completely 
effective, but the presence of the corona discharge around the anode wire indicates 
the occurrence of space charge and a consequent reduction in the field in the 
immediate vicinity of the anode. Workers on positive corona in molecular nitrogen 
(Rubens and Henderson 1940) have shown that the voltage drop through the 
visible discharge to the anode surface is small, being close to the first ionization 
potential of that gas. Thus the actual field in the electrode gap is entirely different 
from that calculated on the assumption of no space-charge distortion, and we 
might well believe that «-particles passing through the gap in the immediate 
vicinity of the anode will not cause spark breakdown. 

On the other hand, it appears from these investigations that an efficiency of the 
order of 90°, is commonly attained. 


§10. DIRECTIONAL EFFECTS 
It has been observed that the counter in all its forms displays marked directional 
effects. In the Broadside position, if the direction of the incident «-particle 
beam is altered to be other than that of the normal to the plane of symmetry, the 
counting rate is observed to fall. Maximum sensitivity occurs when the direction 
of incidence is normal to the plane of symmetry. The counter has a very small 
sensitivity for «-particles whose direction is that of the axes of the electrodes. Ina 
similar way in the End-on position, the counting rate falls as the «-particle beam 
is rotaced in the plane containing the axes of the two electrodes, to make increasing 

angles with the normal to the axes which lies in this plane. 
These effects may be described as follows. In the Broadside case, when an 
a-particle passes through the counter gap normal to the plane of symmetry, the 
electrons produced in the sensitive region can all be expected to converge on the 
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ode. The lateral diffusion of the electrons in the direction along the anode wire 
ill be small, thus a dense local concentration of electrons will arrive at the anode 

d the probability of initiating the spark mechanism will be large. When an 
-particle traverses the gap obliquely, the electrons are no longer concentrated 
y the field to arrive at one locality on the anode surface. The avalanche is 
iffuse and may arrive at the anode over a considerable length of the electrode. 

he more oblique the incidence, the more dispersed along the anode will be the 
avalanche, and under such circumstances the probability of sparking will fall 
sorrespondingly. A similar qualitative description obviously applies in the 
<nd-on case. 


§11. QUALITATIVE DESCRIPTION OF THE COUNTING 
CHARACTERISTICS 


In the pre-visible corona-onset region of voltage, electron avalanches will give 
rise to photons, and these in turn may liberate electrons from the cathode. These 
electrons will produce further avalanches and the likelihood of a second avalanche 
following the first will increase with the voltage. In this way, with large and 
frequent avalanches occurring according to the theory postulated by Meek 
Meek 1940, Loeb and Meek 1940), streamer formation near the anode can be 
expected to occur. In a highly diverging field such as is employed in these 
Sounters, the streamers are unlikely to bridge the gap and so cause spark breakdown, 
Dut a corona will be set up when the critical voltage is reached. At voltages 
somewhat below this critical voltage, it is likely that the dense highly localized 
nization of an «-particle will produce an avalanche large enough to cause a 
spark. ‘The counting threshold is defined as the voltage at which this process 
commences. At this stage the field in the gap will be comparatively free from 
space-charge distortion, and if the field is maintained, the very intense burst of 
»hotons produced by the spark will have a large probability of liberating sufficient 
electrons from the gas molecules and the cathode to produce further avalanches 
of sparking magnitude. That this ‘ after effect’ can persist for relatively long 
ntervals is shown by the necessity of using very long (0:05 sec.) quench times 
vith parallel plate spark counters (Keuffel 1949, Pidd and Madansky 1949). 
[he rapid rise to the high counting rate (peak B) of the Broadside counting 
characteristic is attributed to such effects. 

Conditions alter greatly once the anode corona occurs. At the voltage 
-orresponding to corona onset, the quiescent space charge in the gap becomes 
normously greater than that existing in the pre-corona region of voltage. 
Consequently that part of the columnar ionization of the «-particle which is 
ormed in the gap now finds itself surrounded by relatively large space charges. 
Vieek’s criterion for streamer formation is that the radial field of any localized 
pace charge must be of the same order as the applied field at the location of the 
pace charge. If we consider the dense, highly localized ionization of an «-particle 
n the light of this criterion, we see that the ionization of the «-particle is likely 
o cause such distortion of the quiescent corona space charge as to have a likelihood 
f leading to spark breakdown. ‘The quiescent space charge will swamp any 
mall ionizing event in the electrode gap. Thus we see from the counting 
haracteristic that once the corona onset is reached, the counting rate falls rapidly, 
ndicating that the spurious electron-producing events following the passage of a 
park are no longer of sufficient magnitude as to lead to spark breakdown. It is 
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thought unlikely that any spurious electron avalanche will be able to cause | 
sparking once the trough of the counting characteristic has been passed. These 
views are in accord with the observed dependence of counting rate on specific 
ionization and gas pressure. 

With increasing applied voltage, the a priori probability that a given electron 
avalanche will reach sparking magnitude increases, but at the same time the 
corona space charges in the gap will increase, thereby tending to reduce the 
sparking probability. Finally, the stage is reached with increasing voltage that _ 
the a priori probability of an avalanche leading to a spark can increase only slightly, | 
but corona-space charges will continue to increase in effectiveness, and so it is not | 


surprising to find the counting rate falling somewhat after reaching a maximum, 1 

Similar reasoning can be applied to the End-on case, with the difference that 
here the ionization of the «-particle provides a conducting path between anode — 
and cathode, and we would expect the voltage dependence to be slight as is indeed — 
the case. 

Counters having hollow cathodes yield important information on these points. 
The general shape of characteristic for a hollow cathode counter of about 2mm, 
inside diameter is seen to be similar to that of the curve of Figure 4, but when the 
diameter is increased beyond 3mm., the glow discharge (corona) has its onset 
potential about 250 volts below the spark threshold for «-particles. In two cases 
examined the diameters were 3-25 and 3-4mm. and the respective corona onset 
potentials and spark thresholds were 2,280 and 2,500 volts, and 2,500 and 
2,850 volts, respectively. On amplifying the quench-circuit output about 
one hundredfold, the electron avalanches produced by «-particles at voltages 
below spark threshold could be made to activate a scaler. Under these conditions 
the graph of counting rate against voltage gave a curve similar to that of Figure 4, . 
the peak occurring at corona onset. Taken under normal working conditions the 
spark-counting characteristic showed a slow rise from threshold to a short plateau 
terminated by the onset of continuous discharge. It is significant that the 
plateau of the ‘non-sparking’ curve was found to link up with that of the 
‘ sparking’ curve. 

When the non-sparking pulses were examined on a cathode-ray oscillograph, 
a wide range of pulse size was apparent, the maximum pulse size increasing with 
voltage until the spark threshold was reached. Above this voltage the number of 
non-sparking pulses decreased (the spark rate increasing), until on the plateau 
of the spark curve, no non-spark pulses were seen. 


$125 THE SYMME DRICALNCOUN DER 

The properties of counters having anode and cathode of exactly similar form 
have been examined for the case in which the electrodes are parallel wires. Such 
counters are in general not very satisfactory, in that the counting characteristics 
showed no plateau for quench-circuit time constants ranging from 2 to 400 psec. 
It was not thought worth while to increase the time constant further, so work in 
this field was discontinued after some of the properties of these counters had been 
investigated. It was shown that the voltage interval between the a-counting 
threshold and the onset of continuous discharge was small, and that this interval 
increased with the electrode spacing. It depended on the diameter of the 
electrodes, being a maximum for electrodes of diameter about 0:45 mm. in 
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ratory air. Electrodes of diameters in the range 0-05-1-0 mm. were used in 
investigation. 

Some dependence on the electrode material was also noticed. Electrodes of 

per, brass, iron and platinum were used, the latter material giving the largest 
age interval for a-counting. Considerable dependence on surface conditions 
observed, so for all these reasons such counters were not thought to be of any 
t service. Under no conditions was an anode corona observed. 


§13. THE USE OF GASES OTHER THAN AIR 


The behaviour of the ordinary (asymmetric) type of counter has been examined 
ng several gases and gaseous mixtures. Only the main features of the results 
| be summarized here. 

No counting plateau was obtained with either argon or helium at atmospheric 
ssure. With argon—air mixtures a plateau appears for a small admixture of air, 
| lengthens as the proportion of air is increased, until the typical air 
racteristic is obtained when the mixture is rich in air. In an argon—air 
‘ture containing 1:5°% air the dependence of counting efficiency on specific 
ization is less than in air; in this mixture a fair approximation to a true 
nber—distance curve may be obtained. Water vapour was found to cause a 
erioration in the argon—air counting characteristic, and a drying agent was 
refore usually incorporated in the chamber. The addition of air to helium 
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tuced no improvement in the characteristic. Oxygen is a reasonably suitable 
ater gas: a rather large plateau slope of about 0-2% per volt was recorded 
1 this gas. Nitrogen is unsuitable; in nitrogen a beaded corona forms on the 
le, not a sheath corona as is observed in air, thus the sensitivity is variable 
g the wire. Nitrogen—argon mixtures are unsuitable for similar reasons, 
nitrogen—-oxygen mixtures are generally satisfactory. For small percentages 
xygen (~8%,), the Broadside curve, after the large peak (Figure 9), shows a 
flat plateau region. The counter is somewhat unstable in the region of the 
s peak, but this instability disappears after the corona is established. Carbon 
ide at atmospheric pressure was found to be unsuitable as a counter gas; 
lateau was obtained and no corona was observed. Carbon dioxide—oxygen 
ures were no better, but with carbon dioxide—nitrogen mixtures a poor 
sau of 0-24% per volt was plotted. 
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§14. THE COUNTING OF PARTICLES OTHER THAN a-PARTICLES 


As has been previously mentioned, the counter in air does not regis 
B- or y-rays to any appreciable extent. A direct measurement of the efficien 
of the counter in the End-on position for B-particles was made using a strom 
source of thorium active deposit, and the efficiency was found to be of the orde 
of 2x 10-*. Nevertheless, a counting characteristic for B-particles was ee 
with a plate cathode using the End-on position, the curve having a plateau a 
shape similar to that obtained for «-particles. | 

Gamma-rays produce no observable effect, but the counter was found t 
count ‘ knock-on’ protons in air, the End-on position being used again, with} 
plate cathode. 

The protons were produced by bombarding Styrafoil with the a-particle 
from polonium, and by comparison with the counting rate in a thin windowe' 
Geiger counter, the efficiency of the spark counter for these protons was shown t 
be about 5%. 

Payne (1949) has shown that B-particles and knock-on protons are counted i 
carbon dioxide at atmospheric pressure. This has been verified, but an extensiv 
examination has failed to show even the short plateau mentioned in that paper. 


§15. THE DIFFERENTIAL COUNTER 


In this arrangement, two wire cathode counters are placed one behind th 
other so that «-particles approaching from the Broadside position can pas 
through the gap of each in turn. In the arrangement studied, the separatio 
between the axial planes of the two counters was about 2mm. Separate powe 
supplies were used for the two individual counters and their quench-cirew 
output pulses were fed to a simple anticoincidence unit. In this way, by recordin 
only those «-particles which cause the first counter to spark alone, it was hope 
that a differential range curve could be obtained. On approaching an «-particl 
source to such a counter pair, it was expected that the anticoincidence countin 
rate would reach a maximum and subsequently decrease as the range of th 
particles was traversed. 

Clearly, for successful operation of the arrangement, the individual counter 
must exhibit stable and reproducible counting, and there must be no couplin 
between them. ‘The rear counter should only spark when «-particles actuall 
pass through its gap. It was to eliminate coupling that the separate power suppli¢ 
were used, and the smallest quench capacity consistent with stable operation wé 
employed. The general result of the latter precaution was that all sparks wet 
extremely faint; under these conditions the protons produced in one count 
gap did not produce spark breakdown in the other. Thus it was possible to ha\ 
the front counter sparking rapidly due to a strong source of a-particles who: 
range was too short to reach the rear counter, and the rear counter would ni 
discharge. Figure 10 gives details of the quench circuits employed as well as « 
the anticoincidence circuit, which is essentially a modification of the cathode 
coupled amplifier due to Schmitt (Schmitt 1938 a, b, 1941). From a study of th 
circuit, it will be apparent that, when the front counter sparks alone, the outpt 
to the scaler consists of negative pulses while coincidences and pulses from th 
rear counter alone give positive pulses to the scaler. With the scaler set to regist: 
only negative pulses, anticoincidences of the type required will thus be recorde 
A differential range curve obtained with this arrangement is given in Figure 1 
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at maximum was found to exist in the counting characteristics of both front 
rear counters; and the voltages appropriate to these maxima were used 
»ughout the operation of the differential arrangement. From graphs like that 
n in Figure 11, an ‘extreme’ and an ‘ extrapolated’ range can be defined 
$5). For ThC’ «-particles, the ‘extrapolated’ range so obtained was 
9+ 0-02 cm. (at 15°c. and 760mm. Hg pressure). This is somewhat shorter 
the extrapolated range of 8-56 +0-03 cm. obtained with a single counter in 
Broadside position—and the latter figure is a few tenths of a millimetre less 
fh the standard value. This result indicates that in the double counter a larger 
‘zing event is needed than in the single counter to produce the same probability 
)parking, and this again is probably due to the fact that in the double counter 
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I field is considerably distorted due to the close proximity of the two systems. 
practice therefore, any double counter arrangement must be calibrated with 
tae standard «-particle range, and even with a single counter this is necessary 
the range value obtained varies slightly with the geometry of the system and 
§), quiescent corona current used. 

‘With each arrangement, however, range differences as between «-particle 
eups agree with standard values. Using ThC’ and ThC «-particles, for 
‘imple, the following range differences at 15°c. and 760 mm. Hg were obtained : 
szle counter: extreme 3:90+0-015cm., extrapolated 3-84+0-015 cm. ; 
Gerential counter: extreme 3-86+0-02cm., extrapolated 3-83 +0-02 cm. 
ese values are to be compared with 3-88+0-01cm. and 3-838+0-002cm., 
pectively, deduced from Henderson’s ionization curve (Henderson 1921). 
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ABSTRACT. An instrument is described which was designed to measure the elastic an 
viscous properties of liquids, primarily of solutions of linear high polymers whose elast 
properties had already been observed qualitatively. The liquid is contained in a narro’ 
annulus between two long vertical cylinders. The outer cylinder can be driven in stead 
rotation at a known speed, or in harmonic angular oscillations of known amplitude an 
frequency. In steady rotation, the torque on the stationary inner cylinder is measured : 
different speeds, and interpreted to give a viscosity coefficient as a function of shear stres 
In oscillatory motion, the amplitude of oscillation of the inner cylinder, when constrained b 
a torsion wire, is measured and plotted against the frequency. The shape of the resultin 
graph can be interpreted so as to distinguish and measure different types of elastic behaviou 
The annular gap, the restoring constant of the torsion wire and the amplitude of oscillatio 
of the outer cylinder can be varied in discrete steps. 'The angular velocity of the out 
cylinder can be varied continuously from 0-3 to 3,000 rev/min. in steady rotation, and tt 
frequency can be varied continuously from 0-25 to 25 c/s. in steady oscillation. Ti 
performance of the instrument in both steady-state and oscillatory experiments is assess¢ 
by comparing observed and predicted results for Newtonian liquids of known viscosit 
An example is given of the behaviour of a typical polymer solution in both kinds of exper 
ment, together with a suggested quantitative interpretation of the results. 


§1. INTRODUCTION 
HE observed rheological behaviour of certain solutions of linear hig 
polymers in simple solvents cannot be completely explained in terms of 
single variable viscosity coefficient depending on the rate of strain. Fi 
example, the recoil of air bubbles in a mixture of poly(methyl methacrylate) ar 
cyclohexanone*, contained in a bottle which has been suddenly turned through < 
angle and then brought to rest, is interpreted as direct evidence of liquid elasticit 


* Made, for example, by dissolving 3 gm. ‘ Perspex ’ in 100 ml. solvent. 
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‘he instrument described in the present paper was designed for the purpose of 
asuring both the elastic and viscous properties of such elastico-viscous liquids. 
1€ properties in a simple shearing motion can be measured directly by this 
trument; the properties of the liquid in more complicated types of motion can 
=n be deduced by a process of generalization (Oldroyd 1950), in which further 
deriments of a different type may well be required to resolve ambiguities. 

In order to observe the liquid in (as nearly as possible) a simple shearing 
»tion, a coaxial-cylinder elastoviscometer was envisaged, with an annular gap 
hall in relation to the diameter and with cylinders sufficiently long to make end- 
rrections relatively small. These requirements called for a high degree of 
ecision in the construction and alignment of the cylinders. The simplest type 
experiment in such an instrument, involving a steady relative rotation of the 
?linders at various speeds, can readily be interpreted to give the viscosity of the 
‘aterial as a function of the shear stress. Stability of the flow at high speeds 
mands that the outer cylinder should rotate (Taylor 1922) at known speeds; 
%e torsional couple on the stationary inner cylinder is then to be measured. For 
purpose of detecting and measuring any wall effects which may be shown by 
lymer solutions (Toms 1949), it must be possible to change the width of the 


‘cillation provides a type of motion which is especially amenable to mathematical 
ealysis. If it is possible to make observations on the same liquid at various 
quencies, different types of elastic behaviour can be distinguished and measured. 
vhen the outer cylinder oscillates in a known way, it is not very easy to arrange 
J at the inner cylinder remains perfectly stationary while the forces on it are 
corded. It seemed preferable, therefore, at the expense of some facility in the 
‘athematical analysis of results, to observe the small amplitude of oscillation of 
fe inner cylinder, when constrained by a torsion wire of known elastic constant. 
| he phase difference in the motion of the inner and outer cylinders could also 
i: measured, if required. 

+ Asa result of the above considerations, it was decided to construct an instru- 
sent of the following type. A vertical outer cylinder, surrounded by a water 
-cket for temperature centrol, can either be rotated steadily at any speed within a 
ide range or put into forced harmonic angular oscillations at a wide range of 
own frequencies, with one of several possible known amplitudes. Any one of 
‘ree inner cylinders, with slightly different diameters, can be accurately aligned 
axially with the bore of the outer cylinder. It is possible to measure the torsional 
uple on the inner cylinder in a steady state by means of a torsion wire, and to 
‘serve optically the amplitude of oscillation of the inner cylinder in a state of 
jniform oscillation. The apparatus described in the present paper differs in 
sveral important features from that of Goldberg and Sandvik (1947). 


| {2 DESCRIP TLON 
_ The parts of the elastoviscometer are separately and rigidly fixed on a concrete 


ed. Three major components may be distinguished : the coaxial-cylinder 


ssembly, the driving mechanism, and the measuring devices. 
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(i) The Coaxial-Cylinder Assembly 7 
The general appearance of this unit and the framework in which it is mounte 
represented in Figure 1; important details of the construction are shown 
Figures 2 and 3. For convenience, we shall refer to the inner and outer cylinde 
as the ‘ bob’ and the ‘ pot’ respectively. ; 
The coaxial-cylinder assembly is held in horizontal bridges, or cross-piec 
extending between two vertical pillars 120 cm. in length made from steel rod 
2-5cm. diameter. These pillars are set 25-4cm. apart and normal to the fl 
surface of a heavy steel base-plate. The cross-pieces are of steel plate, 1-Siey 
thick, and they have various widths between 10 and 15cm,; holes are drilled ne 
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their ends so that they slide upon the pillars, being held apart and in position b 
steel sleeves of the correct lengths. ‘The top end of each pillar is threaded so tha 
the whole framework comprising pillars, cross-pieces and sleeves can be locke: 
together rigidly by tightening two hexagonal nuts (Figure 1). This method 


construction makes it possible to dismantle and re-assemble the apparatus whe: 
necessary. | 


The pot is essentially a (vertical) tube closed at the bottom end by a plug, 
downward extension of which forms the driving shaft (Figure 2). The eae | 
axes of this shaft and the bore of the pot body are coincident. The pot is mounte 
in two ball bearings, one at the top and one on the driving shaft. These bearings 
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ch are held in cross-pieces, are coaxial with the bore of the pot body; the pot 
tes freely in them, and the setting is such that its longitudinal axis is as nearly 
ical as possible. 

The body of the pot was made from a solid billet of stainless steel, and the whole 
e inside surface was finally honed to a mirror finish. The liquid-tight vessel, 
Ocm. in depth, obtained by pressing a stainless steel plug into one end of this 
has the form of aright cylinder. Its internal diameter is given in the Table; 
r dimensions can be deduced from Figure 2, which is drawn to scale. 


Dimensions of the Cylinders and Annuli. 


Cylinder Annulus 
Functional Length Width 
Type Designation diameter * (cm.) (cm.) 
(cm.) 
Outer Pot 4-3424 — — 
Inner Bob I 4-2423 25-40 0-0500 
t Inner Bob II 4-1402 25-40 0-1011 
Inner Bob III 3-9383 25-40 0-2021 


The dimensions to which the instrument was made were specified in inches, and the variation 
iameter along the length of a cylinder is less than 0-0001 inch. 
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Figure 3. Details of bearings supporting inner cylinder. 


-ralumin, honed smooth, but the bottom end has been shaped to a cone of large 
ical angle (approximately 145°) terminating in a cylindrical spigot. There is a 
iametrically tapered hole in the centre of the flat end of the bob. The 
‘gitudinal axes of this hole, of the cylindrical body of the bob, and of the 


got, are coincident. Three bobs have so far been made to this pattern: they 


vA bob has the form shown in Figure 2. Basically it is a solid right cylinder of 
| 
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ditfer in diameter (see Table) but have the same length measured from the 
surface to the end of the cylindrical body (25-40 cm.), or to the end of the spi 
(26-80 cm.). ei 

It has been arranged that when a bob is mounted inside the pot the verti 
axes of the two coincide; this is accomplished by means of special bearings wh! 
locate the top and bottom ends of the bob. Details of these bearings are showr 
Figure 3. 

A miniature thrust ball bearing serves to locate the bottom end of the k 
centrally with respect to the bore of the pot. ‘This bearing 1s held in a steel slee 
which fits the spigot of the bob; it bears on a 5 mm. diameter steel ball locate 
a conical depression at the centre of the floor of the pot. A removable screw p! 
forms the central part of the pot floor (Figure 2): replacement of this plug is 
only alteration of the pot that need be occasioned by a change in the form of ' 
bottom bearing. (Similarly, different devices could readily be attached to 1 
spigot of the bob.) 

In order to locate the top end of the bob, a hollow steel spindle with a symmet 
cally tapered end (Figure 2) is inserted into the tapered hole provided ; this spin: 
is locked in place by means of a threaded rod which passes through it and 
screwed into the body of the bob (cf. Edsall et al. 1944). The longitudinal a 
of the spindle then coincides with that of the bob. ‘The bob is lowered into the q 
through a heavy brass ring of rectangular section, held in a cross-piece; the f 
end of the spindle projects upwards through this ring when the bob rests on 1 
bottom bearing (Figures 1 and 2). It is arranged* that the inside cylindri 
surface of the brass ring is coaxial with the bore of the pot. It now remains 
locate the spindle centrally with respect to the brass ring: this is done by mez 
of an air-lubricated bearing which is self-centring and almost ideally frictionle 

The air bearing (Figure 3) has two parts, one a cylindrical boss, made 
Duralumin, which slides on to the spindle attached to the bob (Figure 2), and t 
other, an annular air chamber having twelve radial holes, each of 0-08 cm. diamet 
drilled symmetrically around the inner periphery. The air chamber fits the br: 
ring, and the position of the boss on the vertical spindle is adjusted so that t 
chamber surrounds it like a collar (Figures 3 and 1), with an annular gap 0-00 
cm. in width between them. Air supplied under pressure to the chamber isst 
from the holes and impinges on the boss, keeping the spindle in a central positi 
relative to the ring. ‘The escaping air mingles with the atmosphere above a 
below the gap; the down draught is prevented from entering the pot by a deflect 
skirt at the bottom of the boss. There is symmetry about the vertical axis a 
no resultant torque on the boss. 

Experience has shown that the air bearing will operate satisfactorily when 
excess pressure within the chamber is maintained at 0-5 atmosphere by a sm 
rotary compressor. Under these conditions air lubrication is maintained ever 
a resultant force of 50 gm. wt. is applied normally to the spindle at a point above t 


* The brass ring is held in the cross-piece by means of adjustable locking-screws (Figure 
and 3). A brass disc, 1:2 cm. in thickness, which fits the annular space between the spindle and 
ring, is placed on the free end of the spindle and pushed down until it rests inside the ring. " 
position of the ring in the cross-piece is then correctly fixed by testing that the vertical axis of 
spindle remains coincident with that of the bore uf the pot while the locking screws are tighter 
Deviations of the spindle from the central position are detected by means of two specially shai 
brass plugs which fit the annular space between it and the wall of the pot. The position of the t 
is adjusted until the tightness of fit of the two plugs is the same (as judged by manipulating tk 
with the finger-tips). There is reason for thinking that by this method the centre of the brass 1 
can be located within 0-0003 cm. of the vertical axis of the bore of the pct. The setting of the 1 
is checked at intervals. 
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ss. The excellent lubrication provided by the air film increases the relative 
portance of even the slightest imperfection in the miniature ball bearing which 
*“pports the weight of the bob, but effects due to this cause can be considerably 
duced by lessening the load on the ball bearing (see § 3). 
|) Liquid is introduced into the annular gap between the bob and the pot by 
rcing it upwards through an axial hole in the pot driving shaft (Figure 2). This 
jle bifurcates at the top so as to discharge liquid on either side of the steel ball in 
%e centre of the pot floor, and it can be closed by a tap near the bottom end of the 
jaft. In the end of the shaft, below the tap, there is a tapered socket which 
ceives astandard glass cone attached toa pneumatically operated filling apparatus. 
e level of the liquid in the annular gap can be deduced from observations of its 
el in a vertical sight-tube attached to the filling apparatus. These arrangements 
jake it possible to fill the annular gap with liquid to any desired level without 
Atrapping air bubbles and without altering the relative positions of the bob and 


The temperature of the coaxial cylinder assembly is controlled by surrounding 
e pot with a lagged water-jacket (Figure 1). Water kept at constant temperature 
+ 0-03° at 25°c.) in a storage tank is pumped continuously through this jacket 
= approximately 2-5 litres/min. Leakage of water where the pot driving shaft 
asses through the floor of the jacket is prevented by a Neoprene gland. The 
+mperature of the ambient air is maintained at 20+1°c. During an experiment, 
1e top of the pot can be enclosed by a removable cover (not shown in Figure 1) 
hich serves to improve temperature control within the annulus and to reduce 
Vaporation of the contained liquid. ‘The cover is designed so that it does not 
aterfere with the movement of the pot or of the bob spindle. 


(ui) The Driving Mechanism 


The pot is moved about its vertical axis by connecting its driving shaft (Figure 1) 

‘ith a mechanism actuated by a 0-5 horse-power electric motor (240 v. 50 ¢/s. 
aduction type). "Twokinds of motion are possible: the pot can be rotated steadily 
t any speed between 0-3 and 3000 rev/min. ; alternatively, a simple harmonic 
ngular motion of the pot can be maintained at any frequency between 0-25 and 
5c/s. The amplitude of oscillation can be varied in steps: at present, five 
mplitudes, in the ratios 1 :2 :4 :8 : 16, are available, the smallest corresponding 
> an amplitude of the pot equal to 0-005 radian. 
- The driving mechanism has been designed so that the electric motor runs 
ontinuously at its maximum speed (1,420 rev/min.). Essentially, it causes a vertical 
haft (the ‘ primary shaft’), carrying a flywheel, to rotate steadily at any speed 
etween 15 and 3,000 rev/min. ‘This is accomplished by interposing an H-type 
ear, a clutch, one of three trains of spur gears, and a pair of bevel gears (in that 
rder) between the motor and the primary shaft. ‘The three speed ranges 
rovided by the trains of spur gears overlap, and so, by employing the H-type 
ear, the speed of rotation of the primary shaft is made continuously variable 
etween the limits specified. When speeds less than 15 rev/min. are required, an 
uxiliary 50:1 reduction gear-box is temporarily engaged: this is connected to 
ne primary shaft by means of a Vee-belt and serves to drive another vertical shaft 
he ‘ secondary shaft ’) at any desired speed between 0-3 and 60 rev/min. 

The pot is made to rotate continuously by fixing a pulley on the pot driving 
naft and connecting it, by means of a Vee-belt, with a pulley on either the primary 
r secondary shaft. 

PROC. PHYS. SOC. LXIV, I--B D 
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The pot is given a uniform angular oscillatory motion in the following wa 
The pulley on the top end of the primary shaft has five threaded holes in its A 
horizontal surface at different calculated distances from its centre. Above t ! 
pulley (and just clear of it) is mounted a rectangular Duralumin plate, or carriag 
having two pairs of self-lubricating bushes fixed to its corners. ‘The bushes slic 
on two parallel steel rods, 25 cm. in length and 12 cm. apart, which are horizont 
and are supported at their ends in a frame. ‘The direction of movement of th 
carriage along the rods is perpendicular to the line joining the primary shaft and 4 
pot driving shaft. A slot, fitted with adjustable hardened steel edges, has been cu 
in the carriage at right angles to the direction of movement. A stout steel pin { 
inserted through the slot and screwed into one of the holes in the pulley beneat 
this pin carries two ball races, one above the other, which run on opposite sides « 
the slot. When this connection between the slot and the pulley has been made 
rotation of the primary shaft causes the carriage to move along its supporting roc 
with simple harmonic motion. The amplitude of this motion is constant an 
depends upon the eccentricity of the connecting-pin relative to the axis of th 
primary shaft, but the frequency, which corresponds with the speed of th 
primary shaft, can be varied between 0-25 and 25 c/s.* 

The moving carriage is made to drive the pot by having on it a centrall 
mounted vertical pin which engages in aslot in the free end ofa horizontal Duralumi 
beam fixed to the pot driving shaft. Since the effective length of the bear 
(50-8 cm., equal to the distance between the axes of the primary shaft and the pc 
driving shaft) is large compared with the linear throw of the carriage, the angula 
oscillatory motion given to the pot is almost exactly simple harmonic. The pe 
moves with the same frequency as the carriage and its angular displacement in 
half-cycle is readily calculable. 


(iii) The Measuring Devices 

In both the kinds of experiment that we make at present with the elastovisco 
meter it is only necessary to know the speed with which the pot is moving when it 
motion is steady or uniformly oscillatory. It has therefore been possible to measur’ 
both speeds and frequencies by a relatively simple stroboscopic method. For thi 
purpose repeating patterns of special design have been marked on the pulley whicl 
fits on the pot driving shaft, and on the rim of the flywheel mounted on the primar 
shaft. ‘The speed of rotation of either wheel is then estimated by means of : 
commercial ‘ strobolamp’. ‘The frequency of the oscillatory motion is deduce 
from the speed of the flywheel. The strobolamp is calibrated before eack 
experiment by means of a standard vibrating reed. 

At present it is arranged that one of two other observations can be made 
depending on the type of motion executed by the pot. When this is made t 
oscillate, angular displacement of the bob is constrained by means of a torsion wire 
and the amplitude of the uniform oscillatory motion assumed by the bob unde 
these conditions is then observed as a function of the frequency. When the po 
rotates steadily, the torque on the bob is measured. "These determinations ar 
made directly by means of simple devices, which will now be described. 

A removable boss (Figures 1 and 2), on which are mounted a screw chuck, : 
pulley of radius 2:54 cm. and a small plane mirror, is fixed on the free end of the bol 
spindle above the air bearing. ‘The chuck serves to grip the lower end of a torsiot 


* Higher frequencies (up to 50 c/s.) are attainable when the amplitude of oscillation is small, , 
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e which extends vertically between it and another similar chuck attached to an 
ustable torsion head, rigidly held in a cross-piece which can be clamped to the 
ical sleeves on which it slides (Figure 1). Straight torsion wires up to 30 cm. 
length can be accommodated in this assembly, and wires of three different 
gths and three different diameters have so far been provided. Their ends are 
anently fixed in cylindrical bushes which fit in the collets of the chucks. The 
es were made to our dimensional specification from a beryllium—copper alloy, 
their torsional constants are, approximately, in geometrical progression 
ween 0)-0037 x 10® and 7-25 x 10° dyne cm/radian. The constants of the wires 
d the moments of inertia of the bobs were determined by a dynamical method. 
If, when the pot rotates and its speed is increased, the torque on the bob 
omes too large, the restoring force due to a torsion wire can be supplemented 
ha suitable couple applied to the pulley on the bob spindle. This is done by 
aching two threads to the periphery of the pulley: these pass over auxiliary 
lleys (not shown in Figure 1) which are mounted in brackets fixed on opposite 
es of the cross-piece containing the air bearing. By loading the threads 
ually a known axial couple can be applied to the bob spindle without intro- 
cing a resultant force tending to displace its vertical axis. 
Angular movements of the bob about its vertical axis are observed with the aid 
an optical lever 100cm. in length. A beam of light is reflected from the 
rtical mirror on the boss attached to the bob spindle, and is focused on a 
rizontal curved scale engraved on a flexible strip of translucent bakelite. The 
le is approximately 100 cm. in length and is set so that one division corresponds 
an angular movement of the bob equal to 0:000635 radian. When necessary, 
gular movements of the pot can also be observed by reflecting a beam of light 
pm a plane mirror fixed to the side of the pot driving shaft (Figure 1). 
4 §3. EXPERIMENTAL METHODS 
) It has been found that the bob will respond to very small couples when its 
itting on the bottom bearing is finely adjusted: this is done by raising the bob 
trough a small distance, thus reducing the load on the miniature ball bearing, so 
‘at friction, and the effects due to minute irregularities, are virtually eliminated. 
he bob is lifted by retracting the chuck attached to the torsion head, and its 
‘rtical displacement is observed by means of a simple clock gauge actuated by the 
yse of the chuck on the bob spindle. ‘he vertical displacement is varied 
cording to the type of experiment: it amounts to 0-0025 cm. when the pot is to be 
sade to oscillate, but is less than 0:00025 cm. when rotation of the pot is intended. 
he latter displacement is sufficient to ensure that the weight of the bob ts entirely 
ipported by the torsion wire. Experience has shown that the detrimental effect 
* these adjustments on the centring capacity of the bottom bearing is negligible. 
' When the setting of the bob has been completed the top chuck is locked to the 
»rsion head. 

In the operation of filling the annulus with liquid, sufficient time is allowed for 
1e liquid, and the bob, to attain thermal equilibrium with the pot. ‘The level of 
1e meniscus in the sight-tube of the filling apparatus is then observed with the 
r bearing in operation, the tap in the pot driving shaft is closed, and the filling 
yparatus removed. 

The annulus is filled to different levels according to the type of experiment 
hich it is intended to make. It is completely filled (so that the liquid surface is 
D-2 
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level with the top of the bob) for experiments in which the pot is oscillated ; b 
in experiments involving continuous rotation of the pot, often at high speeds, le 
liquidisused. In these latter experiments, moreover, measurements of the torq 
on the bob, over a suitable range of pot speeds, are always made with two levels 
liquid in the annulus differing by 10-20cm. ‘This is done so that end effects m 
be eliminated (see § 4). 

Although no provision has yet been made for determining the temperature 
the liquid in the annulus while the elastoviscometer is working, measuremen 
made with the pot at rest have shown that at 25°c. it differs by less than 0-2° 
from the known (mean) temperature of the water passing through the jack 
Moreover, there is reason for thinking that the liquid temperature remai 
sensibly constant throughout an experiment in which the pot is oscillated, b 
when the pot rotates at constant speed it sometimes happens that heat is produc 
so rapidly by shearing the liquid that the torque on the bob is observed to chan; 
with time. It then becomes necessary to operate the apparatus discontinuous! 
stopping after each measurement has been made in order that the liquid 
the annulus may cool. It is not possible, however, to accelerate the pot fro 
rest to a steady speed, and to make the necessary observations of speed ar 
torque, in less than 10 seconds. There is therefore a limiting rate of shez 
depending on the viscosity of the liquid, above which reliable observations cann 


be made. 
§4. PERFORMANCE 


(1) Experiments involving Steady Rotation 


The observations made when the pot is rotating steadily are interpreted to gi 
the viscosity of a liquid as a function of shear stress (or of rate of strain). Ws 
effects, if present, can be detected and measured. 

In order to eliminate end effects, two experiments, in which the annular gz 
is filled to different heights (h’, h” measured on an arbitrary scale), are performed « 
each liquid. ‘The couples G’, G” exerted on the stationary bob in a steady stat 
are recorded as functions of the angular velocity Q of the pot, and the relatie 
between torque and angular velocity in each experiment is then defined t 
drawing a smooth curve through a plot of the observed values. For a Newtoni: 
liquid the viscosity 7 is given by the formula 


_ (PF =7yP)[E"Q)—G'(Q)] 1 
n= tart’ 0, 2 (1) 
where 7, is the radius of the bob, and 7, is the internal radius of the pot; C.G.s. unt 
are used in all calculations, and © is measured in radians/sec., so that is obtaine 
in poises. In this case the rate of shear y and the shear stress 7 at r= /(14 19) a 
given by 


— 2ryr,Q 
y= ior? (2) 
iGO) GQ) 
— Qnr_rh" —h’) i PEN Ger 5c (3) 
and, of course, n=T/y. 1 ~—* See (4) 


When the annular gap is narrow and the liquid is non-Newtonian, the sari 
formulae (1)—(3) give the relations between viscosity, rate of shear and shear a 
parametrically, with only a second-order error. In this case + has the sari 
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Sysical significance as before, and y is a first-order approximation to the rate of 
sar at r= /(r, 72) (reckoning the annulus/diameter ratio as a first-order small 
antity), so long as wall effects are absent. 

Wall effects are detected by plotting » against 7 (or y), as given by equations 
(3), for different annular gaps. Any appreciable wall effect would make the 
tves for different annular gaps distinct, but if wall effects are absent the 
es should coincide with only a second-order error. 

Liquids having viscosities between 0-02 and 20 poises approximately, can 
adily be studied. It is possible to make rotational experiments with more viscous 
juids, but two practical difficulties are encountered. In the first place, much 
ne is needed to adjust and determine the level of liquid in the annulus. Secondly, 
lays due to heating of the liquid by shearing are both frequent and protracted, 
en when Bob III is used and the pot rotates at relatively small speeds. The 
‘oblem presented by very mobile liquids (7 <0-02) is naturally different, since the 
yntrolling factors are the sensitivity of the instrument and the stability of laminar 
ow in the annular gap. The conditions for stability of laminar flow between 
axial cylinders when the outer one is rotated have been investigated by Taylor 
~936), and our own observations in experiments with water are broadly consistent 
ith his. ‘The sensitivity of the elastoviscometer in its present form is sufficiently 
dicated by the fact that reliable observations can be made when 7 is as small as 
‘dynes/cm*. Because of these limitations, measurements on liquids with 7 <0-02 
ive to be made within a relatively narrow range of fairly large pot speeds. Never- 
seless, in an experiment with distilled water at 25°c., the value 7 =0-0087 was 
dtained, and when allowance is made for the error of a single determination (see 
Slow), this estimate can be considered to be in good agreement with a standard 
alue, 7 =0-00894 (Bingham and Jackson 1917). 

| The accuracy of a mean value of 7 determined with the elastoviscometer has 
ben assessed by making nine independent experiments with aspecimen of dimethyl 
athalate. The viscosity of this Newtonian liquid, measured at 25-00 + 0-02°c. 
y means of a standard U-tube viscometer (B.5.5. No. 2 pattern, calibrated by the 
‘ational Physical Laboratory), was found to be 0-130 + 0-000, poise; its density at 
1e same temperature was 1:184gm/ml. Three experiments were performed with 
ich bob, and in each experiment measurements were made with three different 
vels of liquid in the annulus. ‘The same torsion wire was used throughout, and 
| was always small enough for the torque on the bob to be measured by means of 
1e wire alone. ‘The temperature obtaining in all the tests was 25:0 +0-2°c. ‘The 
reparations for every independent experiment included a complete readjustment 
f the setting of the bob within the pot. 

The mean values of 7 determined with Bobs I, II, and III were 0-127, 0-128 and 
127,, respectively, indicating that differences due to the replacement of one bob 
y another are without significance. Attaching the same weight to each of the 
ne experiments, we find the mean viscosity 7 to be 0-128 (standard deviation 
002). This standard deviation is not larger than might be expected in view of 
ie number and nature of the possible independent sources of error. 

Individual errors are difficult to assess, even when they are systematic and their 
use is known ; but in those cases where the determination is practicable, no 
ngle error has been observed which is large enough to be a likely source 
_ ambiguity in the kind of work that will ordinarily be done with the 
istoviscometer. 
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(ii) Experiments involving Oscillatory Motion 


Results obtained when the pot is oscillating uniformly are more difficult to— 
interpret. The elastic properties exhibited by any particular liquid are measured — 
by comparing experimental results with theoretical predictions for a series of 
idealized liquids characterized by different elastic constants. Agreement 
between theory and experiment is to be considered satisfactory only if the deduc- 
tions from experiments with one bob and one torsion wire can be used to predict 
correctly the behaviour of the same liquid with different bobs and wires. . 

Oscillatory experiments are carried out with the annular gap completely filled 
with liquid. The amplitude of oscillation of the bob, with a torsion wire in | 
position, is measured at frequencies in the range 0-25 to 25 c/s. The ratio 9 of the 
bob and pot amplitudes is then plotted against the frequency 1, several curves 
corresponding to different bobs and wires being obtained in the case of each liquid. 
A check is made that the resulting curves are (very nearly) independent of the pot 
amplitude when this is small; for practical convenience, most of the work is then 
carried out with a pot amplitude of 0-02 radian. 

Predicted curves have so far been obtained by considering an elastico-viscous 
liquid whose behaviour in shear at small amplitudes is characterized by three 
constants, a coefficient of viscosity y) and two relaxation times 4,, A,. A linear 
relation between the shear stress 7, the rate of shear y, and their rates of change 
7 and y is assumed : 


THAT =Hyly tAsy)-— = a (5) 


It is seen that the relaxation times have the physical significance that if the motion is 
stopped, the shear stress decays as exp (—t/A,), and if the stress is removed the rate 
of strain decays as exp(—t/A,). The relation (5) represents the behaviour of a 
type of material which includes as particular cases the Maxwellian elastico- 
viscous liquid (A, =0) and the Newtonian viscous liquid (A, =0, A, =0). 

In making theoretical predictions, the amplitude of oscillation is supposed 
small, so that the ratio of the amplitudes of the bob and the pot is independent of the 
absolute value of the amplitude, and 7, can be identified with the viscosity observed 
in steady rotation at very low rates of shear. The full length of the annulus is 
assumed filled with liquid, and end effects are ignored. The predicted (9, n) 
curves for different bobs and torsion wires can then be obtained for a series of 
idealized elastico-viscous liquids. ‘The mathematical theory involved will be 
outlined in another paper by one of us (J. G. O.). 

The agreement between theory and experiment in the case of a Newtonian liquid 
(Liquid Paraffin B.P. of density 0-877 gm/ml. and viscosity 1-76 poises at 25° c.) is 
illustrated in Figure 4. Typical experimental points obtained at 25°c. using 
Bob II and a torsion wire of restoring constant 2-79 x 10® dyne cm/radian are 
shown in relation to the theoretical curve. The frequency my of free oscillation of 
the wire and bob is indicated. Inset in the figure are shown some of the individual 
experimental points near the peak, on an increased frequency scale, as an indication 
of consistency. ‘The agreement appears to justify the neglect of end effects in the 
theory of small oscillations, as well as indicating the accuracy which can be expected 
of the apparatus in this kind of experiment. 

It has been found possible to carry out oscillatory experiments with liquids 
having viscosities between 0:05 and 50 poises. In the tests made with liquid 
paraffin it was found that results obtained with the smallest pot amplitude 
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(0-005 radian) were unreliable, apparently for purely mechanical reasons. In 


practice, therefore, an agreement between observations at amplitudes 0-01 and 
0-02 radian is looked for. 


§5. BEHAVIOUR OF A SOLUTION OF A LINEAR POLYMER 
As an example, we have studied a mixture of poly(methyl methacrylate) and 
pyridine at 25°c. This liquid, of density 0-98 gm/ml., contained 30:5 gm. of 


polymer per litre; the intrinsic viscosity [7] of the polymer *, measured in pyridine 
at 25°c., was 5-50. 
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Figure 4. The relation between # and 7 in the case of Liquid Paraffin, B.P. 
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Figure 5. The relation between 4 and y in the case of a mixture of poly (methyl methacrylate) 
and pyridine. 


Figure 5 shows that 7 is not independent of y in the range 0-1<y< 10*sece*; 
Nevertheless, it is reasonable to suppose that at very small rates of shear (y<0-1) 
the viscosity of the mixture attains a constant (limiting) value 7, and this is 
estimated to be 7-9 poises. 

Results obtained in oscillatory experiments are given in Figures6. | Lhe 
curves refer to experiments with three bobs (I, HI, I11) whose moments of inertia 
are 2,677, 2,468 and 2,098 gm./cm?, and two torsion wires (numbers 5A and 8) 
with restoring constants 0-171 x 10° and 2°79 x 10® dyne cm/radian, all carried 
out with a pot amplitude of 0-01 radian. In each case % has a maximum value 


* [y]= ioe (nsp/c), where the polymer concentration ¢c is expressed in gm. per 100 ml. of 


solution. 
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greater than unity, and this occurs at a value of m greater than mq, the appropriate 
frequency of free oscillation of the wire and bob. The contrast with Newtonian 
behaviour is very marked: all (9,») curves for Newtonian liquids, whether 
predicted or observed, show a maximum not exceeding unity at a frequency 
almost exactly equal to mg. 
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Figure 6. The observed relations between $ and w in the case of a mixture of poly (methyl 
methacrylate) and pyridine. 
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Figure7. The predicted relations between $ and n in the case of an idealized elastice-visoous quid. 


The experimental curves in Figure 6 differ only shghtly from those predicted 
for an idealized material of density 0-98: gm ml. characterized by the constants 
0=7-9 (poises), A, = 0-065 (sec.), A, =0-015 (sec.), which are shown in Figure 7. 
It is therefore concluded that, to a first approximation, the mixture of poly (methyl 


* The values of 4, and A, which give the best agreement between theory and experiment in the 
Case of a single bob and torsion Wire can be estimated to within about 5%; the values quoted are 
those giving the best fit m the case of Bob II and Tossion Wire 8. 
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The Maxwellian prototype is inadequate to explain the behaviour of this polymer 
solution. 

Results obtained in both rotational and oscillatory experiments with other 
polymer solutions will be described elsewhere by two of us (D.J.S. and B.A.T.). 
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ABSTRACT. Data are presented on the evaporation rate of § titanium over the tem 
perature range 1,650—1,810° k. 
These are represented by 
: 2-42 = 
log E712 — Ts +9-79, 


where £ is in gm. cm~? sec"! and T is in degrees Kelvin. 
The latent heat of vaporization is deduced both from the measured values of evaporation 
rate, together with known thermodynamic functions of the solid and gas, and also from the 


temperature coefficient of evaporation rate. 
An estimate is made of the variation of total radiation with temperature, and of the 


vapour pressure at the melting point. 
Experimental details are given. 


§1. INTRODUCTION 
T= vapour pressure of titanium is of interest both technologically, since the 
metal is attaining engineering importance, and theoretically, since, for 
example, its temperature coefficient gives the latent heat, and hence the 
binding energy of the lattice. 

It was, therefore, decided, as part of the titanium research programme of the 
Royal Aircraft Establishment, to measure the rate of evaporation and, incidentally, 
the total radiation of the pure metal in vacuo. Some preliminary results and 
measurements have already been published (Carpenter and Reavell 1949, Blocher 
and Campbell 1949) and the present paper gives the final results of the Farnborough 
work with an account of the experimental technique. 
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$2) “THEORY. OP CURE METHOD 
The method of determining the vapour pressure p was, in essence, that of 
Langmuir*, in which the rate of evaporation E of a wire at known temperature T © 
is determined in vacuo. It follows from kinetic theory that 


Fal MEN 2 OOS eee (1) 
a 2n7RT 


where M is the molecular weight of the evaporating particle, R the gas constant 
and « the condensation coefficient. 

Using the Clausius-Clapeyron equation and assuming that the vapour is large 
in volume compared with the solid and behaves as a perfect gas, it can be shown 
that 


d mm 4s a ee 
ay eT qT)" } R” (2) 


where AH, is the heat of sublimation at the temperature T (Guggenheim 1949, 
Speiser and Johnston 1950). Hence, if In« does not change significantly over 
the temperature range of the measurements, it is possible to deduce AH, from the 
variation of E with T. AH, can also be deduced from individual values of p 
provided that the difference of entropy (AS), between solid and gas is known, for 


AG; = SH DNS 5, 
and AG,,, the standard Gibbs free energy of vaporization, is given by 
AGp=—RT In (p/po), 


where fy is the pressure of the standard atmosphere. Hence equation (3) can be 
solved for AH. 

If, during evaporation experiments, the watts supplied to the wire are measured, 
the relation between total radiation and temperature may be estimated. In 
principle, a correction should be made for the cooling of the wire due to the 
evaporation itself, but in practice this is negligible. 


§3. APPARATUS AND EXPERIMENTAL TECHNIQUE 

The experimental procedure consisted in the determination of the weight of 
titanium evaporated from filaments of known surface area which were heated 
electrically in a high vacuum for a known time, the brightness temperatures being 
measured optically. 

The filaments were made from titanium wire, of about 0-6 mm. diameter, 
which had been prepared either from the chloride by the magnesium-reduction 
process or by thermal dissociation of the iodide. 

The magnesium-reduction titanium was prepared in this establishment 
(Bickerdike and Sutcliffe 1949) the analysis (°%, by weight) being 


Titanium 99-68 Silicon 0-07 Iron 0-05 Aluminium 0-02 
Manganese 0-002 Chromium 0)-0035 Cobalt 0-001 Nickel trace 
Oxygen 0-2 Nitrogen 0-09 


The ‘iodide’ titanium, which was probably of higher purity than the 
‘magnesium’ titanium, was prepared in the Philips’ Laboratories at Eindhoven, 
and was kindly supplied by Dr. Fast, to whom we are much indebted. 


* The Knudsen method was inapplicable as a non-reactive crucible could not be found. 
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The wire was dipped in 30°, hydrofluoric acid for 20 seconds to clean the 
surface, and each new filament was heated in vacuo to about 1,800° kK. for a few 
minutes in order to remove volatile impurities, 

The vacuum system was made of Pyrex glass and pressures lower than 10- mm, 
Hg were produced by an oil-diffusion pump. Simple mechanical arrangements 
were used to make dismantling and reassembly easy (Figure 1), the middle of the 
three cone and socket joints being waxed and water-cooled. A 32cm, length of 
wire was bent into a hairpin-shaped filament D the ends being spot-welded on to 
lugs of titanium sheet. In order to avoid sagging, the filament was mounted 


Figuie 1. Evaporation apparatus. 


A. Tungsten rod. G. Optical window. 

B. Terminal block. H. Etched hole. 

C. Nickel screen. J. Pyrex pegs. | 

D. Titanium filament. K. Slotted driving plate. 

E. Pyrex collecting sleeve. L. Flattened supporting rod. 
F. Outer containing vessel. M. Cone and socket seal. 


vertically by clamping the titanium lugs on to a copper terminal-block B which was 
attached by a grub-screw to a short tungsten rod (A). This rod was sealed into the 
lower end of a Pyrex tube, the other end of which was integral with the cone 
portion of the upper of the three ground joints shown in Figure 1. 

Thin platinum washers were used to separate the titanium lugs from the copper, 
since it was found that hot titanium reacts with copper to form a silvery alloy, 
whereas no reaction takes place between titanium and platinum at any rate up to 
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1,000°c. The terminal block was connected to a battery supply by flexible 
copper leads which passed out through copper-to-Pyrex vacuum seals. The 


power was kept constant during an evaporation run by means of a rheostat and 


wattmeter, the duration of the experiment varying from 8 to 1 15 minutes according 
to circumstances. ‘The evaporated metal was collected on a cylindrical Pyrex 
sleeve E, made thin to aid rapid attainment of thermal equilibrium of the filament. 

As the ends of the titanium filament were cooler than the rest, an artifice of the 
‘guard-ring’ type was used to ensure that only evaporated metal appropriate to a 
sharply defined isothermal portion was caught on the sleeve and subsequently 


weighed. ‘To this end the upper few centimetres of the sleeve were lined with a — 


cylindrical nickel screen C which was in place during the evaporation but was 
removed for the weighings. Hence the effective part of the filament was that 
below the plane passing through the bottom of the sleeve. The validity of this. 
artifice depends upon the ‘catch’ per unit area of sleeve (which is of course 
abnormally low opposite the cool portions of the filament) having become uniform 
when the level of the bottom of the screen is reached. ‘To check this the screen C 
was made in three separate portions just touching and it was confirmed by weighing 
that the ‘catch’ per unit area on the lower two was the same. 

To prevent the development of local hot spots, the mass of titanium evaporated 
during a run was restricted to about 7 mg. and was determined by weighing the 
collecting sleeve E, before and after a run, ona balance having asensitivity of 0-1 mg. 
Suitable precautions were taken to avoid errors due to absorption of atmospheric 
water vapour. ‘The evaporating area of the filament was calculated from its initial 
diameter and the masses evaporated in successive runs. Owing to gas pick-up, 
discussed in § 4, it was not possible to use the weight loss or resistance change as a 
measure of the amount evaporated. 

A rectangular hole H (1-6x0-8cm.) was etched through the side of the 
collecting sleeve. The sleeve of diameter 4cm. fitted snugly inside the outer 
containing vessel F and was supported on a greased cone which formed part of the 
lowest ground joint M, the socket being sealed on to the bottom of the outer 
envelope. ‘Thus the inner sleeve could be rotated, from outside, about its own 
vertical axis. ‘The inset to Figure 1 shows the detailed mechanical arrangement of 
this rotating vacuum seal. ‘The flattened steel supporting rod L and the slotted 
driving plate K allowed sufficient movement in two dimensions to obviate fracture 
of the pegs J at the bottom of the Pyrex sleeve. Thus the etched hole could be 
brought opposite a side-arm carried by the outer envelope. The end of this 
side-arm was sealed by a plane window G through which the filament could be 
viewed by an optical pyrometer. When temperature readings were not being taken, 
the inner sleeve was turned so that evaporation metal did not blacken the window. 
Owing to the loss of evaporated metal through the hole in the collecting surface, 
it was necessary to add a small correction (0:8°%) to the mass of titanium collected 
on the sleeve. 

Temperatures were measured with a disappearing filament type of optical. 
pyrometer with a red filter (A=0-64,.). A series of temperature readings of one 
particular part of the filament (the ‘ pyrometer observation station ”) were taken 
during each evaporation run, and the individual values were averaged to give 
the mean temperature of the observation station. Before and after each 
evaporation run, a temperature distribution experiment was performed. For 
this, the outer containing vessel and inner collecting sleeve were replaced by a 
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cylindrical envelope. ‘The filament was heated to a temperature below that at 
which evaporation starts (a brightness temperature of the order of 1,100°c.) and 
the temperature was observed at 0:5 cm. intervals over the effective length of the 
filament. Filaments which, owing to wire-drawing defects, showed variations * 
greater than +10°c. from the mean temperature of the ‘isothermal’ region 
were rejected. The difference between the temperature at the ‘ pyrometer 
observation station’ and the mean temperature of the filament during the 
distribution runs was scaled up, in the ratio of the absolute temperature of the 
evaporation run and that of the distribution run, to obtain the difference between 
the mean of the observed temperatures of the ‘ observation station’ and the 
required mean temperature of the filament during the evaporation run (Table 1). 
The pyrometer was checked against standard tungsten strip-filament lamps which 
had been calibrated at the National Physical Laboratory. Corrections for the 
absorption of the window were determined by calibrations with and without 
the window present. The absorption of the window was checked after each 
evaporation run to ensure that no blackening of the window, due to evaporated 
titanium, had occurred during the run. The observed brightness temperatures 
were converted to true temperatures by means of a linear extrapolation to 
1,810°k. of the known emissivity values for titanium (Bradshaw 1950) f. 

In the preliminary work (Carpenter and Reavell 1949) a constant emissivity 
of 0-5 was assumed, based on unpublished work of Went (Van Arkel 1939) and 
some preliminary results of Bradshaw. The use of Bradshaw’s final results (1950), 
obtained on material from the same manufacturing batch as that used in the 
present work, has slightly altered the temperature scale and, in particular, leads 
to a value of the temperature coefficient of evaporation higher than that previously 
reported. The results of Michels and Wilford (1949) have not been used because 
the anomalous results they give for the resistivity and its temperature coefficient 
suggests that their titanium was of doubtful purity. 


§4. EXPERIMENTAL RESULTS 


Titanium changes from close-packed hexagonal to body-centred cubic 
structure at 1,155 +20°k. (de Boer, Burgers and Fast 1936) and the evaporation 
rate below about 1,500°x. is negligible. Hence, the results given below refer 
to the £ phase. 

Eleven filaments (numbered F.1 to F.11) were used in the experiments, but 
as the nickel ‘ guard screen’ technique was not introduced until F.7 was reached, 
only the results of the last four filaments are recorded here. F.9 was unusable 
as it deformed and accidentally touched the sleeve E during its first run. All 
runs on F.7, F.8, F.10 and F.11 are included in Table 1 and are considered to be 
of approximately equal weight. The evaporation rate of F.10, prepared from 
‘iodide’ titanium, was not significantly different from that of ‘magnesium ’ 
titanium (Figure 2). Actually the results from the earlier filaments (in which the 
‘cold end’ effect was allowed for by computation) are in good agreement with 
the later runs but have nevertheless been discarded as being of inferior reliability. 


* This limit of -+-10°c. was fixed after consideration of the error introduced by using the 
arithmetical mean temperature instead of the effective evaporation temperature, the evaporation 
rate not being linear with 7. 

+ We are indebted to our colleague Mr. F. J. Bradshaw who kindly made his results available to 
us prior to publication. 
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To get rid of impurities more volatile than titanium, a filament on first heating 
was raised to a temperature as high as it was subsequently intended to use and 
about one milligram evaporated off. It was then taken down, reweighed, = 

for use in the first recorded evaporation run. 

In Table | the brightness temperature recorded is corrected for absorption 
of optical window and for errors in pyrometer scale. The true temperature given 
in column 5 is deduced from the brightness temperature using known values of | 
emissivity (Bradshaw 1950). Columns 6 and 7 give the actual masses (corrected _ 
for evaporated material lost through the etched hole) and times, and column ae 
the corresponding rate of evaporation. 


Table 1. Experimental Results 


() 2 S) (®) (5) (6) (7) (8) 
F7A2 1695, 104, +8 1808 8-5 458 4-832 
F8A2 1633, 95. 0 1729 3-7 690 1-071 
F8A3 1571 88 —1 1658 4-6 3240 0-370 
F8A4 1572 88 0 1660 10-2 6900 0-390 

FIOA2 = 1617 93. —5 1706 4-0 1063 1-052 
FIOA3 169%, 104. —5 1796 8-3 466 5-111 
FILA2 1711, 106, —14 1804 10-3 452 6-001 


(1) Run number; (2) brightness temperature (° K.); (3) emissivity correction (° K.); 
(4) distribution correction; (5) true temperature_(° K.); (6) m (mg.); (7) t (sec.); 
(8) Ex 10* (em. cm™ sec). 


Log (“7”) 


10*x1/T (7 in°kK) 


Figure 2. Evaporation rate of titanium. 


It was found that the mass evaporated and caught on the sleeve usually 
exceeded the mass lost by the filament by about 10°, . The additional matter 
was no doubt obtained from the gas phase, as the ok of the ‘ gettering’ 
property of titanium (Carpenter and Reavell 1948), and, as the filament was 
much hotter than the deposit on the collecting sleeve, the former is the more 
probable absorbent. This view was supported i the fact that at the end of its 
life a filament was considerably more brittle than when new, a fact consistent with 
gas pick-up (Fast 1938). Moreover, on several occasions the temperature 
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coefficient of resistance at room temperature was measured on a filament when 
new and after use, and found to have decreased from about 40 to about 30 x 10-4. 
Also the resistance itself increased by more than was consistent with the amount 
evaporated. Both these observations point to the entry of gas into the lattice. 

The probable error in E£ is estimated to be between 3 and 4%, assuming 
the effective area of the filament to be its geometrical area. The major part of 
this error arises from the uncertainty in the evaporated mass. The probable 
error in 7'is estimated to be about 5°, and arises from a combination of calibration 
and reading errors. 

The Battelle results (Blocher and Campbell 1949) and the Farnborough 
results give plots of log p against 1/7 which are approximately parallel, the Battelle 
line lying about 12% above the Farnborough one. This could be accounted for 
by a systematic difference of about 7° in the temperature scales, or by a difference 
in the ratio of true to apparent surface area. 

In Figure 2 the values of log ET"? are plotted against 1/T and a straight line 
has been drawn through the points using the method of least squares. The 
equation of the line so obtained is 


log ET!? = — a 1 O37) ae i re (5) 
and the correlation coefficient is — 0-992. 

The slope of this line leads to a value of AH,,.;=110-9kcal/mole and the 
application of Student’s t-test shows that the probability of the true value of the 
slope lying within +5 kcal. of this is one half. 

The main purpose of the present work was the determination of the vapour 
pressure, but it seems worth while to put on record the data obtained incidentally 
on the variation of total radiation with temperature. The filament voltage and 
current were noted, the filament areas calculated, and the corresponding 
temperatures and values of log JV/A plotted against log 7. ‘The result was a 
straight line. The intercept on the log /V/A axis gives the value of C in the 
equation 

pyatte Co Ce) Ge FA eee: (6) 
and the slope the value of n. 

In calculating A a crude attempt (based on temperature distribution) was made 
to allow for the presence of the cold ends of the filament, as it was not thought 
worth while to attach special potential leads and thus define a definite length of 
filament at a uniform temperature. The resulting uncertainty in C in equation (6) 
probably lies between 5 and 10%, but the value of is considered to be reliable. 
The final result is 

atts cies loo eed ee gan (7) 


and is based on a plot of 31 points. 


§5. DISCUSSION OF RESULTS 
(i) The Value of AHyog 
Equations (3) and (4) show that AH, can be deduced from individual values 
of the vapour pressure provided that the difference of entropy of solid and gas 
is known. 
Prof. L. Brewer has computed values of H, and of the function — (Gp — Hy9)/T 
for the solid and the gas from specific heat data compiled by Dr. K. K. Kelley 
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of the United States Bureau of Mines, Berkeley, California, and his results * 
are given in Table 2. , 
Using linear interpolation we have obtained the appropriate value of 
—(AGp—AHg9,)/T for each of our experimental temperatures. Hence we 
construct Table 3, computing AG, using Table 1, equation (1), assuming « =1 
and monatomic vapour, and equation (4), and thus obtaining AZo. a 
The assumption that «=1 is not inconsistent with the data; the point is 
discussed later. 
Assuming that there is no significant error in the values of the function 
—(Gp—AHyg)/T derived from Table 2, the effects on AH, of given errors in. 


Table 2. Thermodynamic Functions for Titanium 


298° kK. 1000° k. 2000° k. 
[—(G—HAlags)/T eas (cal. deg—* mole) 43-07 45 -84 48-72 
[—(G—Ho98)/T |sotia (cal. deg? mole?) 7-24 10:52 14-70 
[—(AG—AAie98)/T Jevap (cal. deg—* mole) 35-83 Sores 34-02 
[H — Hagslgas (kcal.mole~*) 0) 3:70 8-90 


Table 3. Values of AHygog calculated from Individual Experimental Points 


(1) (2) (3) (4) (5) (6) 
F8.A3 1658 5-1781 55-44 34-47 112-59 
F8.A4 1660 5-2015 oR! 34-47 142-52 

F10.A2 1706 5-6381 53-43 34-41 112.12 
F8.A2 1729 56487 54-07 34:37 113-50 
F10.A3 1796 4-3357 50-54 34-29 112-13 
F11.A2 1804 44063 50-17 34-28 112-00 
F7.A2 1808 43127 51-06 34-27 113-02 


AHeyo= 112-55. 
S:D= 20:55 


(1) Run number; (2) temperature (° K.); (3) log (ET?) (c.c.s.); (4) AG (kcal/mole) ; 
(5) —(AG—AAig95)/T (cal/mole); (6) AHog (kcal/mole). 


our experimentally determined values of E and T can be calculated. If we take 
(see § 4) the probable error of E as 3-5°/,, and that of Tas 5°K. the probable error 
of AH, should not exceed 0:4 kcal. This is consistent with the standard 
deviation of 0-55 kcal. given in Table 3. 

Using Brewer’s results for the heat content of the gas as a function of 
temperature together with those of Kelley for the solid (Kelley 1949) one obtains 


AHsg, —AHive, 33:85 kcal. ee (8) 


Hence, adding this to the value of AH,;3; obtained from equation (5) (viz. 
110-8 kcal.) we obtain AH,,=114-7 kcal. as opposed to 112-6kcal. which is the 
mean value obtained in Table 3. 


* We are much obliged to Professor Brewer for kindly making the figures available to us prior 
to publication. 
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This discrepancy of 2 kcal. between the value AHyog, deduced by the method 
of Table 3, and that obtained from the slope of the line in Figure 2 is well within 
the +5kcal. estimated as the probable error from Student’s t-test. The slope 
method is clearly of inferior accuracy when only a small number of points are 
available. 


(ii) The Condensation Coefficient 


The present experiments do not provide a sensitive test for the value of «. 
As, however, not much is known about this quantity for metals at high temperature 
it is worth while to draw what conclusions are possible. 

The information available can be used in three ways. We can examine 
Figure 2 for evidence of variation of « with temperature; we can examine 
‘Table 3 for evidence of the actual value of «, assuming it does not vary with T, 
and we can enquire what assumed constant value of « would be consistent with the 
difference between (AF 29s) temp. coeff. and (AZZ 508 )ontropy . 

In the first place, equation (2) can be integrated and the total change of In«w 
obtained over the experimental range T, to T,, provided that AH, is known as a 
function of T. Now, from thermodynamics, (d/dT) (AH) =AC,, where AC, is 
the difference of molar heat of solid and vapour, which we will assume independent 


of T. Then 
sane) ML 21), SG a). 


~T's T. 
d(Ina) = | 
| 7 ( oe) iis T; R tes ihe 


We will insert experimental values between T, =1,660 and T,=1,800°x. The 
first term with an estimate of its maximum error can be obtained from Figure 2 
and is 2-65+0-23; the second term is —2-57 and its maximum error (corre- 
sponding to a maximum uncertainty of about 2-5 kcal. in AA yg69) is +0-06. The 


third term is found to be negligible, so that the truth of the assumption that AC, 
1800 


is temperature independent is unimportant. Now physically, | xe d(In«) is 
almost certainly negative, and the most negative value it can reasonably have is 
(2:65 —0-23 —2:57 —0-06) = —0-21. Hence the possible decrease of « from 
1,660° to 1,800° k. cannot exceed about 20%. 

Secondly, if « is constant but less than unity, it is in principle possible to 
discover the fact by examining the variation of (AH go) enizopy 28 deduced from 
evaporation observations over a range of temperatures. It follows from 
equations (1), (3) and (4) that the true value of (AH o9¢) entropy exceeds by RT In « 
the value deduced on the assumption that «=1 from evaporation observations 
at temperature 7. The effect is small. For example, if «= 3, the difference in 
the values of (AA os)entropy Calculated from evaporation observations at 1,660° 
and 1,800°x. respectively, is 0:28 kcal/mole. 

Hence very accurate determinations of evaporation rates and entropy data 
would be necessary before a systematic trend of (AH a9) entropy With temperature 
could be attributed to a departure of « from unity. Examination of the last 
column of Table 3 shows that our accuracy was not sufficient to reveal the 
existence of such a trend. ‘ 

Thirdly, as mentioned in the preceding section, the value of (AH 59s) temp. coer. 
derived from the slope of the line in Figure 2 is about 2% higher than that of 


(AH oos)entropy> Which is calculated from the entropy data combined with individual 
E 
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values of the evaporation rate, assuming «=1. However (AHoog)temp. cooff. has a 
probable error of 5% compared with less than $% for (AHo9¢)entropy: eee: | 
a=} and 1/10, the calculated values of (AH29¢)entropy decrease by 2 and 8% 
respectively. Thus, assuming «=1/10 would make a discrepancy between 
(AF oo8)temp. coe. 224 (AH o298) entropy of nearly twice the probable error of the former. 

Summarizing, the experimental data are consistent with « being as low as 
4 but not as low as 1/10, while a temperature variation over the temperature range 


of more than 20% is not likely. 


(iii) The Vapour Pressure 


Using equation (1) to convert equation (5) to atmospheres on the assumption 
that « =1, we obtain 


log Patm= — = x 104-+- 7:30 -= ~~ Tiare (5 a) 


Linear extrapolation gives a value of 1:2x10-?mm. Hg for the vapour 
pressure at the melting point, 2,000 + 10°. (Van Arkel 1939). 


(iv) The Radiation Characteristics 


The total radiation results do not call for much discussion except to note 
that for metals the value of usually exceeds the black body value of 4 (see for 
example Blanc 1921) and the equation (7) implies a total emissivity for titantum 
at 1,735°x. of about 32° of that of a black body. 
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ABSTRACT. The description is given of an electrostatic focusing system by which 
electrons emitted from the tip of a hairpin filament are reunited at a distance of 2 to 3 cm. 
No intermediate potentials are required. An x-ray method is used for the experimental 
investigation of its properties. The electron spot is a reduced image of the tip of the 


_ filament and has a diameter of about 40 ». On the basis of this system, a fine focus x-ray 


tube is designed, with a peak specific target loading of about 11 kw/mm?. Some performance 
data are given. 


Sie UNGER OD Cale lOAN 

HE development of x-ray analysis and its application to problems demanding 
| increased intensity and resolution led early to the construction of tubes 
with rotating anodes (e.g. Miiller 1929, Astbury and Preston 1934, Miller 
and Clay 1939). It has been found possible in this way to increase the specific 
load W (the energy per unit area per second) applied to the anode by a factor 
5 to 10 as compared with tubes employing a stationary anode and a focal spot 

of the same dimensions. 

There is, however, another approach to the problem of increasing the specific 
load in the focus; an investigation of the heat dissipation in the stationary, water- 
cooled anode of an x-ray tube shows that as the radius R of the focal spot, assumed 
to be circular, is reduced to dimensions much smaller than the thickness of the 
anode, then W «1/R for a given target material (Goldsztaub 1947, Oosterkamp 
1948). Thus a focus of R=15 may safely be loaded 30 to 40 times as high as 
one of normal dimensions, with a corresponding increase in X-ray intensity. 
Other advantages of a fine x-ray source in connection with crystallography have 
been pointed out by Bernal (1945). The production of fine, highly loaded electron 
foci is a problem of electron optics; in view of the considerable advantage of such 
foci in x-ray tubes, it is surprising to find that comparatively few attempts have been 
made to apply electron optical ideas to this problem. 

Demountable fine focus tubes employing magnetic focusing have been 
described by Goldsztaub (1947) and Metropolitan-Vickers Electrical Co., Ltd. 
(1948). The tubes follow closely in construction and design the example of the 
cathode-ray tube or the electron microscope. ‘lhe crossover, or a first image, 
produced by the electron gun, is demagnified by a magnetic lens whose image 
plane coincides with the anticathode. ‘The results so obtained are good, but the 
need for a stabilized u.7. and lens current, and therefore complicated power 
unit, is a serious disadvantage. In addition, the object-image distance has to be 
comparatively large to obtain the required demagnification, which also makes the 
focusing system very sensitive to stray fields. 

This points to the use of electrostatic focusing systems which do not require 
a stabilized u.T. unit. One would like this system to have the following 


* This paper forms part of a Ph.D. thesis presented by W. E. S. to the University of London. 
E-2 
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properties: (a) concentration of electrons to a spot of less than 0-1 mm. diameter 
at the maximum permissible load, (b) no intermediate potentials which would 
complicate the H.T. unit, (c) simple external adjustment to optimum focus, 
(d) comparatively short filament-anode distance. 

The object of this paper is to give an account of a simple electrostatic focusing 
system, and to describe its application to a fine focus x-ray tube. 


§2. THE FOCUSING SYSTEM 
A preliminary experimental and theoretical study of the more conventional 
electrostatic lens systems (Myers 1939, Zworykin et al. 1945) led to the conclusion 
that the electrostatic focusing system of an x-ray tube presents a distinctly separate 
electron optical problem. ‘The system eventually developed and examined in 
detail is shown in Figure 1. It consists of the grid cylinder G, mounted coaxially 


Figure 1. Section through focusing system. 
Dimensions are given in cm. 


with the anode rod A. G surrounds the hairpin filament F, whose tip lies on the 
axis of symmetry of the system. C is a cathode shield, containing a central 
aperture of 0-25 cm. diameter. Anode and grid cylinder are enclosed by a coaxial 
electrode E, which also forms the outer wall of the tube and carries the x-ray 
window W. ‘The H.T. is applied between A and G; the filament is kept at a 
small positive potential V, with respect to the grid cylinder. The electrode E 
forms an integral part of the electron optical system and the focusing properties 
are affected fundamentally by the potential applied to it. The system is suitable 
for the intended purpose if E is at or near the potential of G. The experimental 
data given in the following section refer to the case Vy=Vg. The spot size 
observed on the anode surface depends critically on the grid aperture—anode 
distance d, and the position of the tip of the filament with respect to the plane of 
the aperture. The bias V,, of the order of 50v., contributes little to the actual 
focusing ; inaddition to controlling the beam current, it largely suppresses emission 
from the sides of the filament. Electron concentrations of half width 30 are 
obtained on the anode surface, at V, =50 kv. and beam currents up to 500 ja. 


§3. EXPERIMENTAL INVESTIGATION OF THE FOCUS 
The intensity of the x-rays produced at any point of a metal anode is proportional 
to the electron intensity, so that, by means of an x-ray pinhole camera, a magnified 
image of the electron spot can be observed on a fluorescent screen or recorded 
photographically. 
In fact, it was found advantageous to use a fine slit rather than a pinhole for 
imaging the focus. Obviously, for a faithful representation, the pinhole should be 
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small compared with the electron spot. It was found difficult to make suitable 
pinholes of diameter appreciably less than 30, but good slits could be made 
well below 10. The x-ray slit used in the experiments (Ehrenberg 1949) 
consists of two slightly etched silver steel rods of 4 mm. diameter, pressed together 
by set screws in a small brass frame and separated by a pair of spacers (74 nylon 
sheet). ‘The use of the slit technique has the further advantage that the intensity 
distribution across the slit image may be measured with a normal microphotometer. 
It should be noted, however, that the slit image will only represent the true 
distribution of intensity across the centre of the spot if the distribution is of the 
form 4 exp(—r’), where A and B are constants. As this condition is very 


nearly satisfied in the present case, the intensity curves given may be taken as 
essentially correct. 


Mahine Ds 2. 00. 150 


Rip) 
Figure 3. Intensity distribution curves of Figure 4. Variation of spot size, 2R, with grid 
three of the slit images shown in Figure 2. aperture—anode distance d for [(R)=0°5 Ip, 


(half width), and J(R)=0°1 Ip. 


Photographic records were usually made on Kodak R20 recording film, with 
a magnification of 10, and aslit—spot distance of 6cm. In order to obtain the spot 
dimensions without foreshortening, the slit was set in a plane normal to the surface 
of the anticathode. With a 35mm. camera up to 40 exposures could be taken 
in succession and developed under the same conditions. ‘The width of the slit 
was measured on a photograph taken with the slit placed close to the film. 

Figure 2 (Plate) shows a series of slit images, recorded to investigate the effect 
of changing the grid aperture—anode distance d. For each value of d (given in 
the central column of Figure 2) two slit images are shown, one exposed for 
5 minutes, the other for 25 minutes. In all exposures Vy, =50 ky., and Vy is 
adjusted to give a beam current of 300ja. ‘The plane of the filament is horizontal. 
In the case of x-rays the photographic density D 1s directly proportional to the 
X-ray intensity up to about D=1-5. The maximum density obtained in the 
5 minute exposure is less than unity and these photographs may therefore be 
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regarded as faithfully representing relative x-ray intensities. The same applies 
to the outer regions of the 25 minute exposure, which is used to investigate the 
“tail’ of the intensity distribution. The densities were measured with a Hilger 
Spectral Photometer. Figure 3 shows intensity distribution curves for some of 
the slit images. Their variation with d appears to be similar to that which would 
be produced by a lens system, imaging a point on the axis with positive spherical 
aberration. Thus curve a, for the largest value of d, corresponds to a section of 
the beam beyond the paraxial image plane. Curve 6 represents a section close 
to the paraxial image plane; a fine central peak, of diameter less than 10 p is 
formed, and the tails have become less intense. Curve c shows a section near the 
circle of least confusion. As expected, the peak of the intensity distribution 
curve is now wider, and a further widening is observed on decreasing d. 

The change in spot size with d is more clearly brought out in Figure4. The 
curves give spot diameters 2R measured at J/J,)=0-5, and J/J)=0-1, where J is 
the intensity at R microns from the centre, and J) the intensity at the centre. 
The minimum half width is about 30 p, with fairly large ‘tails’ of 125 « diameter 
at0-1J). At the circle of least confusion, these are reduced to 85, but the half 
width is now 45 p. 

Further series of slit photographs were taken to investigate the effect of 
asymmetry of the filament and the change in spot size and quality with the beam 
current and with the filament position relative to the aperture. It is found that 
the spot is slightly elliptic, the length of the axes differing by about 25%. The 
spot size is practically constant for beam currents between 100 ma and 500ya. 
The position of the filament, however, is critical. ‘The results show that the tip 
of the filament should always lie in the plane of the aperture, or slightly above; 
withdrawing it into the aperture leads to a considerable increase in the intensity 
thrown into the tails of the distribution curves. 


$4. THE NATURE OF THE ELECTRON CONCENTRATION 

In order to understand the formation and the nature of the electron concentration 
on the anode, the potential distribution produced by the electrodes shown in 
Figure 1 was determined by means of the relaxation method for the cases Vp=V 
(system I), and Vyg=V, (system II). The filament was idealized to a thin 
circular cylinder, with a flat base in the plane of the aperture. Space charge was 
neglected. The potential distributions form the basis for the determination of 
the electron trajectories through the systems. A more detailed account of this 
investigation will be given in a second publication; here it is sufficient to state 
some of the general results. 

(a) In system I, the concentration of electrons on the anode surface is due 
to the formation of a demagnified image of the emitting area of the filament. This 
is clearly brought out by Figure 5, which shows a set of trajectories leaving the 
axial point, and one off-axis point of the emitting surface at 42° to the axis and 
with different initial velocities. The rays are denoted by a group of three numbers ; 
the first represents the velocity of emission (x 10-7 cm/sec.), the second the angle 
of emission, and the third the initial radial height in millimetres. The most 
probable velocity of emission is 3-6 x 107 cm/sec. A demagnified image is formed 
in a plane situated a small distance beyond the anode surface. The magnification 
is 0-32. Experimental values of spot size and object-image distance agree closely 
with the theoretical results. The experimentally determined optimum position 
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of the anode does not coincide with the image plane found, due to the occurrence 
of a small group of high velocity electrons emitted at angles outside the range 
of the ray tracing. These rise to comparatively large radial heights, are over- 


focused by the field, and thus define a plane of least confusion slightly in front 
of the image plane. 


7 (mm) 


aed ! 
Emitting Anode _ Paraxial 
surface Image Plane 


Figure 5. Electron trajectories through system I. Distances along the axis of the system (z-axis) 
are measured in ‘ units’, where 1 unit = 0-790 mm.; radial distances are in mm. Rays are 
denoted by a group of three numbers: the first represents the velocity of emission 
(x 10-7 cm/sec.), the second the angle of emission, and the third the initial radial height 
in mm. 


rT (mm) 


-0:12 Anode 
Emitting Surface 


Figure 6. Trajectories through system II. Units and notation as in Figure 5. 


(b) Figure 6 shows a similar series of rays traced through the field of system II 
whose focusing action resembles closely that of the more conventional type of 
immersion lens. It can be seen that an image is formed with considerable 
longitudinal aberration close to the grid aperture. ‘The rays, having passed the 
image plane, are now bent back towards the axis, forming a disc on the anode 
surface whose diameter is strongly dependent on the velocity of emission. 
Obviously such a method of concentration of the electrons is unsuitable for use 
in a fine focus tube. 


72 W. Ehrenberg and W. E. Spear 


§5. THE DESIGN OF THE FINE FOCUS TUBE 
Figure 7(a) shows a complete section through the tube based on system I. 


The dimensions of the focusing system are given in Figure 1. ‘The main body, 


supported by the exhaust tube 'T’, carries a metal cone at its lower end which fits 
directly into the top of the oil diffusion pump P. To the right of T, the internal 
diameter of the outer tube is slightly increased to accommodate the glass cylinder 
G, sealed into the metal body with Apiezon wax. E is an ebonite ring, waxed on 
to the glass, which effectively suppresses any corona. ; 
The cone joint C carrying the anode A is sealed to the other end of the glass 
cylinder. The anode is soft-soldered into the cone after the assembly of the tube 
has been completed, by the use of a jig which standardizes its position. It carries 
a disc of the target metal, of thickness 3mm. brazed into itsend. Several anodes 
with different target metals can be prepared and centred for use with a particular 
tube. Contrary to conventional x-ray tube design, the anode is at H.T., for the 
metal outer wall and the external focusing adjustment should be at earth potential. 


(b) linch in Figure 7a) 


M 


(€) td) 


Figure 7. Section through fine focus x-ray tube. 


For this reason the anode is cooled by a slow stream of carbon tetrachloride. 
The circulating system consists of a small centrifugal pump and a reservoir; 
connections are made with flexible plastic tubing. 

The whole of the cathode assembly is detachable from the main body of the 
tube. A vacuum-tight joint between the two parts, which terminate in accurately 
machined flanges, is provided by slight compression of the rubber ring R. The 
nickel cathode cylinder is held centrally in the main tube by the brass barrel, 
shown in detail in Figure 7(6). It is turned to give a good sliding fit in the outer 
tube and contains a simple mechanism which facilitates lateral and horizontal 
adjustment of the tip of the filament in the grid aperture. The nickel tube N 
carries 2 small ceramic filament holder shown in Figure 7(d): the filament is 
made from 0-006 in. tungsten wire, bent into the shape of a hairpin which is 
spot-welded to the nickel supports M of the filament holder. The tip should be 
sharply bent, and the sides close to it almost parallel. The height of the filament 
tip above the supports M is about 3mm. ‘Tube N is locked by a set screw to 
the brass disc D, which has the shape shown in the section, Figure 7(c). After 
inserting a new holder with filament, it is only necessary to loosen screws H and 
readjust disc D by gently pushing it with a screwdriver through the holes provided, 
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until the new filament lies centrally within the grid aperture. The grid cylinder 
is then adjusted axially until the tip lies in the plane of the grid aperture. The 
filament leads (not shown) are clipped to the nickel supports of the filament 
holder and pass through the tube N to the terminals F, riveted to a mica disc. 

The cathode assembly is connected by rod J to the top of the bellows B. 
The cathode cylinder will therefore move towards the anode if the bellows are 
compressed by turning screw S. It is thus possible to alter the grid aperture— 
anode distance d whilst the tube is running, and it becomes a simple matter to 
refocus visually after the filament has been replaced. 

Figure 8 (Plate) shows a photograph of one of the tubes which forms part of 
asmall table unit. The supporting tube is held here by a bracket screwed to the 
steel frame of the table. This method of fixture allows adjustment of the tube 
height above the table top. The high-tension lead is led through the steel frame 
and table top in an ebonite cylinder, which also surrounds it above the table. 
A spring-loaded rod makes contact with the anode part of the tube. 

Figure 9 shows the simple circuit that has been found satisfactory. T, is a 
compact H.T. transformer which incorporates metal rectifiers for half-wave 
rectification. It has been designed for this tube and gives a maximum current 
of lma. at 50kv. The bias V, is produced by passing the beam current through 
the variable resistor R,. This arrangement provides a certain amount of negative 
feedback which considerably stabilizes the beam current. 


E 
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Figure 9. Diagrain of electric circuit. 


T,: H.T. transformer, 50 kv., 1 ma., incor- lie & variable resistor, 300 kQ. 4 
porating metal rectifiers. iby 2 variable resistor, 2,500 Q, 25 w. 

T.: filament transforme1, 5 amp., 5 v. R;: resistor, 300 kQ, 25 w. 

C,: condenser, 0-05 pr., 50 kv. working. M,: moving-iron voltmeter, 0-250 v. 

C,: condenser, 0:01 pF. M,: moving-coil microammeter, 0-500 pa. 

R,: potentiometer, 1,000 Q, 200 w. M;: moving-iron ammeter, 0-5 amp. 

R,: resistor, 100kQ. N: neon lamp, VS 70. 


The tube described has now been in use, on and off, for over one year, and only 
very slight anode pitting has been detected. It is run at 50kv. with a copper 
anode and with beam currents between 300 and 500ya. For a half width of 40, 
this leads to a mean specific load (load within circular area of 40 diameter, divided 
by area) of 8kw/mm2. The peak specific load (at the centre of the spot) is. 
11kw/mm?. Calculation shows that under these conditions the maximum 
temperature on the anode surface rises to approximately 450° c. For comparison 
it may be mentioned that typical values of the load in a commercial sealed-off tube 
and a rotating anode tube are 0-1 kw/mm’, and 0:5-1-0 kw/mm? respectively. 
These two figures are probably too high, because in arriving at them no allowance 
is made for the tails of the distribution. 
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The question arises whether the combination of fine focus and rotating anode | 
would not permit a further substantial increase in the specific load. Itis possible 
to arrive at an estimate of the improvement by making the following assumptions: 
(a) The load W is distributed uniformly over the focal spot. (b) The anode is 
regarded as a semi-infinite, homogeneous and isotropic solid of thermal con- 
ductivity k, density p and specific heat s, the boundary plane of which coincides 
withthe anode surface. The assumption is justified for fine foci, whose dimensions 
are much smaller than those of the anode. (c) The maximum temperature 
attained at a point P, moving with velocity v in the centre of the focal path, will be 
approximately the same as that produced by a stationary load W after a time | 
t=2R/v, i.e. the time it takes P to traverse the focus. | 

In view of assumption (6), the solution of the time-dependent equation of heat _ 
flow will be of the form | 


_ 2W(spy"2.d A dt 


pro eno 
ar= 4 kent exp ( — spr?/4kt) 


(Oosterkamp 1948); dT denotes the contribution to the temperature at any 
point P arising from an instantaneous energy inflow dE =2W dA dt, t seconds 
before the time under test andrcm. from P. For an element in the centre of the ~ 
uniformly loaded circular focal area dA =2nrdr. Assumption (c) makes the 
variable r independent of t, and the temperature T, in the centre of the focal path 
of the rotating anode, is found by integrating dT from r=0 to r=R and t=0 to 
t=2v/R. This leads to 
WR 1 

T, =~ <p {1 exp (—22) + om! — O(a), 
where « =(Ruvsp/8k)!”, and @(«) denotes the well-known error integral. 

To the same approximation, the temperature 7, in the centre of a stationary 
anode is obtained from the above expression by putting v=0, so that T, = WR/R. 
For the comparatively high value of v=785 cm/sec., corresponding to a disc of 
5cm. diameter making 3,000 rev/min., and a focal spot of radius 20 pu, «=0-43, 
if the material constants of copper are used. T/T, now becomes 1-2 for a given 
value of W. At the given maximum temperature, the loading of the rotating 
anode is therefore only 20°, higher than that of a stationary anode. The technical 
difficulties of increasing the speed of rotation appreciably above 50 rev/sec. are 
considerable, particularly as it is essential that the anode should run true to within 
a fraction of the focal diameter. 

The results show that the improvement in loading that can be obtained by 
rotating the anode of a fine focus tube with approximately circular spot is in no way 


comparable with that for focal spots of normal size, and seems hardly to warrant 
the additional complication and cost. 


§6. THE PERFORMANCE OF THE TUBE 
Figure 10 (Plate)* shows a shadowgraph of a platinum grid supporting two 
10 platinum wires at right angles. The grid was placed against the window of 
the tube and the shadowgraph taken at a distance giving a magnification of 20. 
The sharp shadow of the vertical wire is due to the effective horizontal source 
width of 5 produced by foreshortening. Although the vertical half width of 


* See end of issue. 
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the source is at least of the order of the size of the object, the shadow of the hori- 
zontal wire is still fairly well defined; its definition depends, of course, on the grid 
aperture—anode distance. Observation of the shadow image provides a simple 
method of adjusting the tube visually to optimum focus. Figure 10 also indicates 
the possibility of using the tube for micro-radiography. 

The successful application of the tube to problems in x-ray crystallography 
involves essentially the design of a special set of cameras, for focal spot and camera 
really form one single unit (Bernal 1928). To use the highly loaded fine focus 
to its full advantage, both camera and specimen have to be small. A programme 
of camera design has only recently been started, but the results obtained are 
encouraging. One of the cameras used so far is of the Seemann—Bohlin type 
(radius2 inches). It is clamped rigidly to the body of the tube, close to the window. 
For this camera no collimator or pinhole is necessary. A fully exposed photograph 
of the diffraction pattern of ZnO was obtained in one minute using Cu-radiation 
and a Ni-filter. In this diagram the «, and a, lines of the (102) reflection 
(i.e. Bragg angle =23-9°) are resolved. Using the same camera and a General 
Electric Ca-6 diffraction tube, an exposure of 2 hours was required to obtain an 
equivalent photograph. 
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ABSTRACT. Heterochromatic light is treated and the spectral sensitivity of the human 
eye is taken into account. In spite of certain simplifying assumptions, conclusions of a 
general character may be drawn, the most important being that the optical properties 
required for low reflectance coatings can be obtained by treating the incident light as. 
monochromatic. 


Sil IND RODU CALTON 
IE HE use of thin coatings for reducing the reflection of monochromatic light at 


glass surfaces in optical instruments has been the subject of extensive 


theoretical research by various workers (Vasicek 1947, Caballero 1947, 
Deli 1949). The reflection, considered as a function of nt (m is the refractive 
index and ¢ the geometrical thickness) of each coating has its minimum at ut =A/4, 
or odd multiples of A/4, where 4 is the wavelength. This minimum reflection 
becomes zero for the proper choice of refractive indices ; for example, for a single 
coating of refractive index n on a substrate of refractive index n,, the proper choice 
is 2 =(Ngn,)"?, where ny is the refractive index of the outer medium, usually air (for 
which m ~ 1). 

§2. PROBLEM 


The light radiated from sources used in practice is in general heterochromatic, 
with the complication that, if the optical thickness is chosen to destroy reflection at 
a certain wavelength, there still remains a certain amount of reflection at adjacent 
wavelengths (Jacobs 1943). ‘The total amount of reflection, i.e. the reflected 
intensity over all wavelengths, may still be a large fraction of the incident intensity. 
An important aim in an optical design is to make this fraction as small as possible, 
thereby reducing the loss in brightness of the object which is viewed. 

A further complication that enters into the problem of improving the efficiency 
of an optical system, and one that is often neglected, is that the human eye is not 
equally sensitive to equal radiant energy fluxes at different wavelengths. We take 
this into account when defining the ‘reflection R’, that is the fractional loss in. 
the useful incident radiant energy, by the equation 


__ Visible intensity of reflected light 
Visible intensity of incident light ° 


Mathematically, this may be written in the form 


i a(X)b(A)p(A)dA 
R= eer 
| , a)o)aA 


* Now Director of the Physical Laboratory of the University of Pretoria. 
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where A is the wavelength, ), the lower limit and A, the upper limit of the visibility 
spectrum, a(\) the energy distribution of the incident light, p(A) the reflecting 
power of the optical system, i.e. the fraction of the incident light reflected at 
wavelength A, and b(A) the spectral sensitivity of the average eye to equal radiant 
energy fluxes. 

The factor p(A) takes care of the effect of the coatings which are used to reduce 
reflection. This effect is characterized in the final functional result by the 
appearance of the refractive index m and optical thickness nt of each coating as 
variables. ‘The present problem is that of finding this functional result and its 
minimum with respect to these variables. 

In order to evaluate the reflection R by using equation (1) it is necessary to 
know the explicit forms of the functions a(A), b(A) and p(A). 

The curve of (A), known as the relative luminosity curve, is shown in Figure 1 


- 68 700, 


X (ny) 


Figure 1. Relative luminosity curve (A). Curves a, b and ¢ are respectively the graphs of b(A), 
the symmetrical part F(A) of b(A) and f(A)=(A)— F(A). Curves d and e are the graphs of 
the analytical representations of F(A) and f(A) respectively. 


and has been drawn by using the tables given by Hardy and Perrin (1932). The 
analytic expression 


B(A) = {1 2 3(4") “ a(t) F +B,{1 a (yh. ee (2) 


where Ay=555 mp, By=45, A,=616mp, 8B, =20 and B,=0-185, is seen from 
Figure 1 to give a curve which is a good substitute for the experimental curve. 

The reflecting power p(A) has been evaluated by various workers (Vasicek 1947, 
Caballero 1947) for single, double and multiple coatings. For simplicity we 
treat two cases of normal incidence. ‘These are as follows : 


Case (a). A single coating. 
In this case (Vasicek 1947, Caballero 1947) 
Ad bese A Be 


pA) = aR SU 2 ee ape (3) 
UT gS, IE oe ala 
where a rane Ueaiprere b= r ( {,say) 


and mp, 7, Ny and t have been defined above. Further, it has been assumed that 
Ny <N<Ng. 
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Case (b). Approximate treatment for a number of coatings. 

When there is more than one coating, the expressions analogous to (3) are less. __ 
simple and the corresponding calculations of R extremely complicated. An- 
approximate calculation can be carried out by neglecting all light beams which are 
reflected more than twice. Let there be k successive coatings. If the incident. 
beam is monochromatic the reflected wave is given in amplitude and phase by 


ah, el =k he? We yd Le eee 


=ytyemt+ 1... ye ert . 
where L-=(teei—1,){(tesate,) $ Opn, for r=0, 1,72 ak 
and 8, =41(myt, + Mgtgt+ .... +n, t,)/A (=a,/A, say) . 
The subscript 7 (=1, 2,....) refers to the rth coating and the subscripts 0 and 


k+1 to the outer medium and substrate respectively. 
Putting «,—«,=«,,, the reflected intensity is 


k s=k 
p(A) =1]* = male +2 es cos (a<y/A), olans Cones (4) 


where /* is the complex conjugate of /. 

The energy distribution a(A) of the incident light is in general completely 
different for different light sources. In (1), however, a(A) occurs only in the 
product a(A)b(A). This product plays an important réle and, like (A), is 
(practically) everywhere zero except in the region A,;<A<A,, where it is positive. 
In certain practical cases the curve of this product strongly resembles that of b(A) in 
departing very little from symmetry and in having a single maximum which 
is in general at a different wavelength (A,, say) from that of b(A). Useful inform- 
ation on the reflection that arises when the product has this typical form may be 
obtained by treating a(A) as a constant, in which case A,, =A. 

The reflection then becomes 


fe b(A)p(A)ad 


where A, and A, are replaced by — co and © respectively. This is permissible 
because (practically) 6(A) =0 for A<A, and for A>A,,. 


§3. EVALUATION OF THE REFLECTION 


Case (a). It is convenient to begin with the case of the approximate equation 
(4). By substituting (2) and (4) into (5) we find 


k s=k 
R,= Ly Pp 2g 2 yl ae ee (6) 
r=0 


s>rso 


here Tyg = (OS Aci) Ah 
where rs | 14+3(5) +3 (54) ca ie (=) 
Tic 3 Bo 3 Bo — oo By 
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These integrals have been evaluated by means of the theory of residues. We find 
T,,,= P{exp (— ay%s,) Sin ($6 — byx5,) + EXP ( — dy%5p) SiN (Fh + dom en)} 

+O {exp ( — a3x,,.) SiN G — byt) +exp (—a,x,,) sin G + bits) eee (7) 
where 

P=37B,/4p?" sind, Q=478,B,, p=1-5¥2, tang =2312, 

ay =sin by/FpP,, G, = Cos $y/pB, 5 (p =1,2;¥ =0), (p =3,4, v= 1), 

012 = (Bo? +p" sin $4)?-+p cos? 44, tan yp, =p! cos 44 /(xoBo +p! sin 4), 

0,2 = (xoByt — pt? sin $4)? +p cos? 44, tan yy =p"!? cos 44/(o(%-1— pl? sin 44), 


of meBy t+ 26141, tan y= (2! By4 +1) 
og =x,°By *—2"?x, By +1, tan of =(20?x,8, + —1)-4, 
and Iq = relete=o i 2M? ar {OU 2B (24? he Ure ap B, By}. 
Case (b). Equation (3) can be written in the form 
pA)=1—H+2H X(—1)(AB)’cosv8, i... (8) 
v=0 
where H =(A?+ B?—1— A?B?)/(1 — A?B?). 


On substituting (2) and (8) in (5) we find 
R=1—-H+2HI,;4 & (—1)(ABYI, 
v=() 


where J, is of the same form as J,,. except that «,, must be replaced by va (see (6) ) 
and I,,=J). On substituting from (7) and summing the series one finds 
sin $¢+ ABe~“™ sin (4¢ + dnt) 
1+ A*B%e- ent + 2A Be cos dnt 
sin $¢ + ABe~“™ sin (4¢ — dant) 
1 + A? B%e-*"t +. 2A Be cogs zat 


R=1-H+2HPIy4 


oHOL- sin a+ ABe ™ sin (fm +dsngt) _ 
23 Os 1+ A2B%e-2et +. 2. A Bee cos dant 
sin 47 + ABe¢s" sin (4a — d,nt) \ 
1 + A? B%e-*s"! + 2A Bes" cos dnt 


where c,=4ra, and d,=4nb,; r=1, 2,3 and 4. 


When nt is taken as a single independent variable, m still remains a second 
independent variable of R. That this is the case follows from the fact that H, A 
and B are, as seen from (8) and (3), functions of n. 

The coordinates on the graph of R(nt), curve aof Figure 2, have been calculated 
for an air (m)=1) outer medium, a glass (, = 1-6) substrate and a coating with a 
refractive index of n = (npn)? = 1-2649. 

In order to obtain some idea of the magnitude of the error involved in the 
derivation of the approximate formula (6) a few points on the graph of R(nt) were 
obtained by using (6). These points are indicated in Figure 2 by circles. 
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The broken curve b in Figure 2 is the curve of R(nt) when b(A), which in 
equation (2) is of the form b(A) = F(A) + SA); is taken to be b(A) = F(A) —f(). Re. | 
curve is useful for the purpose of investigating the influence of the asymmetrical — 
form of (A) on the location of the minimum of R(nt). ae ; 

The graph of R(n), when nt is fixed at 138-1 mp (Ao/4), is shown in Figure 3. 


§4. REMARKS AND CONCLUSIONS cs 
It is seen, under the conditions set out below equation (4), that the minimum re- 
flection occurs for practically the same values of m and nt, namely n= (1gt2)" 2 = 1-268 
and A,/4=138 mp respectively, as in the case of monochromatic light of wave- 
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i 2. Graph of the reflection R plotted against nt for an air (m)=1) outer medium, a glass 
me Gaia) eubsteate and a coating of refractive index n= (NoMa)! *—1-2649. Curves a and b 
are obtained from equation (9) and correspond respectively to b(A)=F(A)-+f(A) and 
b(A)=F(A)—f(A) (see Figure 1). The circles represent points obtained from the 
approximate equation (6). 
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Figure 3. Graph of the reflection R plotted against m for nt=138-1 mp (Ao/4), m=1 and m.=1°6. 


length Ay, and this is true in spite of a slight asymmetry in the curve of a(A)b(A). 
Also, when the actual value of nt in a practical case is different from A,/4 by as much 
as 40 my and v is different from (mpn.)"? by 119%, the reflection will still be small, 
less than about3%. The minimum reflection, namely 0-1%, is negligible. 

Under the above mentioned conditions the rule of considering the incident light 
as a monochromatic beam of wavelength A, is therefore satisfactory when a choice 
must be made of the proper values of m and nt for a single low reflectance coating. 
Also, these results are not inconsistent with the idea that the same rule may be 
applied to obtain good estimates in the case of multiple coatings. It must be 
pointed out that, though the actual calculations were made for glass, the rule 
applies equally to other optical materials. 

It may also be mentioned that the approximate formula may be used to attain 


high accuracy in the case of single coatings. Itis probably less accurate for multiple 
coatings. 


Low Reflectance Coatings and the Sensitivity Curve of the Eye 81 


Once the coating material is chosen and the thickness appropriate to the given 
light source determined, the problem of forming the film still remains. An 
important part of this problem is that of measuring the film thickness. "The most 
simple way of doing this is by making observations on the interference colours and 
using the relationship between the optical thickness and interference colour. A 
mathematical examination of this relationship has now been made by Kubota 
(1950), who gives a plot of his results on the I.C.I. Diagram. 
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The Interpretation of Dielectric Measurements using the 
Cole—Cole Plot 


Stark (1950) has considered the equivalence of frequency and temperature in dielectric 
measurements with particular reference to the Cole—Cole plot (Cole and Cole 1941); we 
here consider the use of this representation in the interpretation of dielectric measurements 
and compare the plots using frequency or temperature as the parameter. A plot of the 
imaginary part of the dielectric constant €” against the real part «’ for a given temperature, 
with frequency as the parameter, almost always results in a curve which is the arc of a circle. 
In the rare case when the dispersion can be represented by a single relaxation time the arc 
is a semi-circle. 

This diagram is useful in several ways. Suppose for example that only a part of the 
dispersion is covered by the frequency range available, as is often the case, then by drawing 


the best fitting arc of the Cole—Cole plot one is able to deduce, with an accuracy depending 
on the proportion of the total arc determined, both the zero-frequency dielectric constant €’9 
and the limiting high-frequency value ¢’,,, assuming that only one absorption occurs in the 
total frequency range. In addition the angle is a measure of the distribution of relaxation 
times (Kauzmann 1942). This representation enables one to see at a glance the extent of 
the distribution of relaxation times and its variation with temperature, whereas the usual 
method requires calculation and curve fitting. Further, when one has two or more domains 
of dispersion which may not be widely separated in frequency, as for instance in the case of 
solid hydrogen bromide and deuterium bromide below the transition at —184° c. (Powles 
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1950 a, b), a slight extension of the analysis using the Cole—Cole plot enables one to ) 
oerane the values of ¢’) and e’.. and of the distribution of relaxation times of each | 
dispersion as if it existed separately. The details of this analysis will be published elsewhere — 


(Powles 1950 c). . 
In the representation at fixed frequency with temperature as the parameter much more 


care is required in interpretation: in fact useful information can be deduced only in parti- 
cular cases. If the temperature range, which should, in theory, cover all temperatures, 
includes a phase transition of any kind, for instance liquid — gas, this type of plot is quite 
inappropriate; a plot with frequency as the parameter on the other hand is valid at any 
temperature. Evidently one may obtain an arc only if the frequency of measurement is 
such that the dispersion region passes through this frequency for a range of temperature 
sufficiently far removed from any phase transition, and not necessarily even in this case. 


The end points of the arcs have no meaning if only because of the variation of the dielectric 
constant due to thermal expansion or contraction. The effect mentioned above is quite easy 
to see analytically. The expression for the Cole—Cole arc is 


e*¥—€'g 1 

€9— Ee 1+ (jurt)?-*’ 
where e*=e’+je”. It is not, however, possible in general to eliminate the temperature as a 
parameter even with the aid of a relation such as log 77=ST)+C as proposed by Stark, 
since €’) and €’.. are also functions of temperature. 

The angle w7/2, which is a measure of the spread of relaxation times in the constant | 
temperature diagram, evidently only represents in the constant frequency diagram some sort 
of average of the spread of relaxation times over the whole temperature range used to 
determine the arc. In fact in nearly all cases the spread of relaxation times is a function of 
temperature. If in a particular case we find an arc of a circle the angle az/2 might be 
considered as the qualitative spread of relaxation times. Care is necessary, however, as a 
particular variation of €’9, €’~ and the spread of relaxation times might combine to 
produce an arc in the Cole—Cole plot and thus lead to misinterpretation. 

There is at present no theoretical reason why so many substances have a distribution of 
relaxation times corresponding to an arc in the Cole—Cole plot (there are in fact many 
exceptions (Kauzmann 1942)); we therefore stress the purely practical utility of this 
representation. 


Department of Theoretical Physics, J. G. PowLes 
The University, Liverpool, 7. 
29th September 1950. 


Coxe, K. S., and Corz, R. H., 1941, 7. Chem. Phys., 9, 341. 

KauzMann, W., 1942, Rev. Mod. Phys., 14, 12. 

Powtss. J. G., 1950a, Nature, Lond., 165, 686; 1950b, C. R. Acad. Sci., Paris, 230, 836; 
1950¢, F. Chim. Phys., in the press. 
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Characteristics of Selenium Photocells 


While the forward characteristics of rectifier cells are not yet clearly understood, it has. 
been shown (Landsberg 1949) how their reverse characteristics can be analysed accurately. 
The application of this method (previously outlined for a Mott barrier) to a Schottky barrier 
leads, with the previous notation, to reverse characteristics of the form * 


} “i Fe aN LES 
1 _y%expy, wherey=(—2—~) | 
Jo Vo 


3 ») 1/2 7 3 4 
ion = (=) nveW* exp ( “2 and V,;= lids 
TT 
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e is the dielectric constant of the semiconductor. The purpose of this note is to point out 


that the characteristics recently given by Preston (1950) are in excellent agreement with the 
theoretical relation, as shown by the accompanying Figure. In this case zinc oxide films 
were deposited on the selenium by sputtering for varying times t, while a gold electrode 


was deposited on the oxide film by sputtering for 7 seconds. The numerical constants. 


* The detailed theory and the interpretation of the results will be published elsewhere. 
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which can be derived from the curves assuming room temperature, Vpj=} volt, and «=6:3 
{Tammann and Béhme 1931), are given in the Table. They reproduce the main features. 


10 


e C= Ssec. 

* yaar Experimental 
10? oe eee Po Points 

© ¢ = 20sec. 

- ¢ =20sec. 


Theoretical Curve 
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obtained by analysing the curves given by Williams and Thompson (1941) for selenium 
rectifiers *, except that they show in addition that the mobility of the current carriers decreases 


(sec.) 5 10 15 20 120 
n (cm-* x 10?8) 0-65 1-46 1-46 1-46 1-46 
Av (x 10-3) 40 1-19 0:59 0-26 0-099 


t=time of sputtering, n=volume density of impurity centres, d=area of contact 
(11-3 cm? +), v=mobility of current carriers in cm?/volt sec. 


as the thickness of the zinc oxide film increases. ‘The present and previous notes show the 
advantage to be gained by plotting reverse characteristics in a form suitable for comparison 


with the curve z=x?/2e”. 


Department of Natural Philosophy, P. 'T. LANDSBERG. 


University of Aberdeen. 
7th September 1950. 


LanpssBerG, P. T., 1949, Nature, Lond., 164, 967. 

PresTON, J. S., 1950, Proc. Roy. Soc. A, 202, 449 (Figure 6). 

TAMMANN, G., and BOHME, W., 1931, Z. anorg. Chem., 197, 1. 

Wituiams, A. L., and THompson, L. E., 1941, 7. Instn. Elect. Engrs., Pt. 1, 88, 353 (Figure 5). 
* The detailed theory and the interpretation of the results will be published elsewhere. 
+ Private communication from Dr. Preston. 


The Measurement of Ultrasonic Vibration Potentials 
(Debye Effect) with Pulse Techniques{ 


Debye (1933) proposed a method involving ultrasonic waves for the determination of 
the relative masses of electrolytic ions in solution. If an electrolytic solution is irradiated 
with ultrasonic waves, periodic differences in potential should appear within the solution 
primarily as a result of the dynamical reaction of the ions to the material wave. The 
existence of this effect, however, was not established experimentally for fifteen years 


+ This research has been partially supported by the U.S. Office of Naval Research under 


Contract No. N7 onr 47002, Project No. NR 051 162. 
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following its prediction because of the small magnitude of the alternating potentials © 
(10-® volts for unit velocity amplitude) and the problem of preventing electromagnetic _ 
coupling between the ultrasonic generator and the detection apparatus. In 1948 them | 
authors (Yeager et al. 1949) detected the effect with standing wave techniques but were | 
limited to semi-quantitative measurements because of the problem of determining the 
acoustical parameters in a complex field. 

Hunter (1950) has indicated his interest in investigating the Debye effect with pulse 
techniques. Apparatus for the measurement of the Debye effect and related phenomena 
with pulse-modulated ultrasonic waves has been developed during the past two years in 
the Ultrasonic Research Laboratory at Western Reserve University (Yeager et al. 1949, 
Ultrasonics Research Laboratory 1949) and measurements are in progress at the present 
time. The use of pulse techniques affords two advantages in this application : 
(i) acoustically produced effects can be differentiated from electromagnetic coupling 
because of the finite time required for the propagation of a material wave, and 
(ii) measurements can be made under virtually free-field conditions (Yeager et al. 1949) 
because the alternating potentials resulting from the directly incident waves can be 
differentiated from those produced by reflections. 

Preliminary measurements have been made with pulse-modulated ultrasonic waves 
at 200 kc/s. with a pulse width of 400 microseconds and a repetition frequency of 400 per 
second in a large tank which permits delay times of approximately 1,000 microseconds. 
‘The alternating potentials are amplified with a conventional tuned amplifier and viewed 
with a cathode-ray oscilloscope. 'The magnitude of the effect has been measured by the ~ 
substitution method with signals of known amplitude. The acoustical velocity amplitudes 
have been determined with a barium titanate hydrophone (Ultrasonic Research Laboratory 
1949, Bugosh et al. 1949) which has been calibrated at the Naval Underwater Sound 
Reference Laboratory at Orlando, Florida (U.S.A.). 

The effect has been found to be approximately 5 microvolts for unit velocity amplitude 
in 0-005 normal potassium chloride. The amplitude of the effect has also been shown to 
be essentially independent of the concentration of the electrolyte for dilute solutions and 
dependent on the nature of the electrolyte. These characteristics are in general agreement 
with those found previously with standing wave techniques. 

The use of pulse-modulated ultrasonic waves has also been extended to a similar effect 
in colloidal systems (Rutgers 1938, 1946, Vidts 1945, Hermans 1938 a, b, Causse 1950, 
Rutgers and Vidts 1950). Measurements with colloidal suspensions of silica with low 
ionic conductivities are in progress and appear to be of the same order of magnitude as 
obtained by Rutgers and Vidts with arsenic trisulphide sols (Rutgers and Vidts 1950). 

Details concerning the experimental methods as well as the results will be published 
in the near future. 


Ultrasonic Research Laboratory, E. YEAGER. 
Department of Chemistry, J. Bucosu. 
Western Reserve University, F. Hovorka. 


Cleveland 6, Ohio, U.S.A. 
Ist October 1950. 
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A Note on Resonance Damping, at Injection, in Betatrons 
and Synchrotrons 


It has not yet been satisfactorily explained how an electron can be injected into betatrons 
and synchrotrons providing maximum particle energies of the order of 50 Mev., without 
colliding with the injector structure on subsequent revolutions. Experimental work 
(Westendorp and Elder 1949) has shown, however, that the presence of large circulating 
currents are essential for efficient injection, and a semi-quantitative theory has been 
developed by Kerst (1948) based on this fact and on the assumption that while the 
circulating current in such machines increases, it induces a corresponding change in the 
orbit linking flux, resulting in deceleration of the electrons, and a consequent shrinkage 
of the radius of their instantaneous orbits. Apart from the fact that the application of 
these ideas to explain internal injection appears unconvincing, the extra energy provided 
by the circulating currents stored in the magnetic field is reversible; it is likely, therefore, 
that an appreciable fraction of the energy initially lost by the electron will be subsequently 
regained within a very few revolutions, as the heavy injection currents dissipate themselves. 
The present note proposes a mechanism which may assist in the successful operation of 
a large number of betatrons and synchrotrons. 

Consider electron motion in an azimuthally inhomogeneous axial magnetic field, 
H{r,¢). The azimuthal error field will be assumed to be a circular function of ¢ with 
a period 1// times that of orbital motion of the electron. The ‘ free’ vertical and radial 
oscillations are then given, with some approximation, by 


er COL Y a0 eee rye ae mieaags (1) 
dx 


where a change in the independent variable has been made. The parameters a and q are 
related to the axial magnetic field index n, to the magnitude and phase of the ‘ error’ 
fields, and to the integer /. 

Equation (1) is the usual form for the generalized Mathieu equation, its solution being 
stable (periodic or non-periodic, and bounded) or, in the general case, unstable (€—> 00, y—> 00) 
depending on the relative magnitude of the parameters a and g. In particular if a=1, 
27, 37, 4°, . . . , solutions of (1) are always stable provided g=0, 1.e. if the axial magnetic 
field is azimuthally homogeneous. For 0<n<1, the most important cases of general 
instability arise for a=1, J=1, or (1—n)!”. For oscillations wholly in the median plane, 
i.e. radial oscillations, the condition a=l=1 gives n=#, while for ‘free’ vertical 
oscillations the same conditions give n=+. The case a=1, /=(1—n)'” occurs only in the 
case of vertical oscillations, giving n=+. 

For a= 1, the general solution of (1) can be written 


€=Ae’% sin(y—o)+Be"Xsin(x+o) sn anee (2) 


where A and B are arbitrary constants much less than 1, and terms of the order of g® are 
neglected. The constants w and o are related to each other and to a and q by the equations 


w= $q?—(a—1)"}12 £0 
2a0= cos{(1—a)/q]= sin- [—2n/q]. 


It is obvious therefore that ¢ is stable and bounded only if «4 is imaginary, i.e. g?<(a—1)? 5 
the boundary between the unstable regions is thus given by the equation 


g=(a—1)* or +tq=a—-1. =  — eaneee (4) 
If (=, dl/dy=n , y=y/2, when ¢=0, we have, writing A=7B, and putting a=1, 
uicl = pee cot tryeatn) ae (5a) 
0<pel 
[20/6 p+0 cecot(tythr) eens (5) 
O0<u€1 


Equation (5 a) gives the initial conditions resulting in A=0, B+0, while equation (5 b) 
gives initial conditions resulting in 440, B=0, the former specifying initial conditions 
giving exponental damping of €, and the latter initial conditions giving exponential divergence 
of ¢. The constant y represents the phase of the inhomogeneity relative to the azimuthal 
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situation (6=0) of the injector, and yw is directly proportional (a=1) to the magnitude © 
of the inhomogeneity. In fact, from (3), “=q/2, giving o= —7/4. 

For the particular case of vertical oscillations (=z, so that ¢’=2’ represents the height _ 
of the injector above the median plane and 7 is the angle of injection measured relative to 
the plane parallel to the median plane and containing the injector. Thus when a=1, 
n=, we have the case of resonance between the ‘free’ vertical oscillations and the 
fundamental of the electronic orbital motion in the accelerator. It follows from (5 a) that 
the most favourable conditions for resonance damping occur when 7=0, y= —7/2, from 


which 
no)=4+—asing. 5 « =. | sg eee (6) 


Thus, if the field index varies azimuthally as (6), damping of the ‘ free ’ vertical oscillations 
follow. If 0<7r<1, however, the oscillation amplitude will reach a minimum and 
subsequently diverge so that, for effective damping, « in (6) must be a large quantity, 
consistent with ¢<1, but must fall to zero after a period very short in comparison with the 
period of the- acceleration cycle. These are precisely the conditions found necessary by 
Jones et al. (1950) for successful injection into a non-ferromagnetic synchrotron. 

For similar resonance damping of free radial oscillations, it can be shown that the 
azimuthal variation of 2 must be of the form 


n(o)=2--0 sind, eee (7) 


when ¢ is proportional to the displacement of an electron from its instantaneous circle, 
and 7 is the angle subtended by its trajectory at injection. 

Assume that the large currents circulating at injection provide such a coupling 
coefficient «. It is almost certain that a large number of electrons do not complete even one 
revolution. The charge and current density will thus be very large in front of the injector, 
and will fall to a minimum behind the injector after completing one revolution. The 
charge and current density will, therefore, have a large first harmonic component. The 
charge density produces an electrostatic (space charge) field which has the same effect 
as an increase in the field index n. The magnitude of this effect on 7 is of the order of 
(c/v)? times that of the opposite effect on m produced by the field generated by the current 
elements, where v is the electron velocity at injection and c/v ~2. Thus a in (7) will be 
positive; the damping condition (5 a) is therefore satisfied, and the circulating space charge 
at injection (lasting for a period short relative to the period of acceleration cycle) should 
initiate resonance damping of the * free’ radial oscillations. 

Thus, for example, if «=0-05 we have 4=0:05 and g=0-01, so that the damping per 
revolution is given by : 
Ar 


(rinj—"i) 


where rjnj and 7; are the radii of the circle on which the injector is situated, and of the 
instantaneous circle respectively. For | Ying —7i | =2 cm., Ar=3 mm. 

Qualitative arguments, based on this theory, can be given to show that injection 
efficiency tends to increase with the energy of the injected electrons when the coupling 
constant is large. A further important point is that such damping is not reversible. 

As a consequence of these considerations, and the fact that a large number of betatrons 
and synchrotrons have values for the field index in the neighbourhood of 3, over a large 
part of doughnut space, it is suggested that the mechanism of resonance damping of the 
‘free’ radial oscillations, brought into operation through the creation of favourable 
inhomogeneities by the large currents circulating at injection, plays an important part in 
the successful internal and external injection of electrons. Details of this theory will shortly 
be published elsewhere. 


oe 0°05er 2 0:15! 9 Ee (8) 


Research Department, S. E. BarDEN. 
Metropolitan-Vickers Electrical Co., Ltd., 
Manchester. 


9th October 1950. 
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REVIEWS OF BOOKS 


La théorie des images optiques; Colloques internationaux. Pp.322. 1st Edition. 
(Paris: Service des Publications du Centre National de la Recherche 
Scientifique, 1949.) 2,000 fr. 


This collection of papers represents a series of contributions by various authors to the 
proceedings of the International Reunion held in the Institut d’Optique, Paris, in October 
1946. It owes much to careful editing by Professor P. Fleury, Dr. A. Maréchal, and 
Mme C. Anglade. The subjects cover practically the whole field of Technical Optics, 
and the book is interesting as giving a survey of the state of optical science at the end of the 
second world war, not a few war-time researches being brought together (though the 
majority of the topics form the subjects of fuller papers published elsewhere). 

When the armistice of November 1918 brought the first world war to an end, it was found 
that pressure of circumstances had compelled a fresh interest in Technical Optics, especially 
in England, France and America. But German workers had established a long lead in this 
field and had already begun to make progress along a good many lines which have since been 
followed up elsewhere. The works of Gleichen, Czapski, Eppenstein, von Rohr and others 
were studied in translation and otherwise, and renewed attention was also given to non- 
German writers such as Dennis Taylor. On the whole, however, the practice of lens design 
was dominated by traditional and empirical methods helped out by much ray-tracing. In 
spite of the pioneer work of Rayleigh, Strehl, and others, little was known quantitatively 
of the physical characteristics of aberration effects; and the tolerances for aberration residuals 
in calculation were based entirely on geometrical aberrational intercepts. Few designers 
were producing adequate formulae without recourse to trial lenses in which separations and 
lens forms were modified in numerous practical trials. ‘The actual use in instruments of 
aspheric surfaces (apart from hand-figured systems) was almost unknown, although fairly 
numerous attempts had been made to produce machines for their production, and A. Koenig 
with M. von Rohr had given formulae for the calculation of primary aberrations in systems 
containing them. Physiological optics had already become a fairly well-developed study 
but the conditions and limitations under which really valid and useful results could be 
obtained were only too often misunderstood, and the implication of the phenomena of vision 
and photography in regard to the design of optical instruments had received little attention. 

It is very interesting then, for those of us who know something of the optical field as long 
ago as 1918 to look through the volume under review. It opens with a series of papers on the 
classification and calculation of aberrations. Mr. 'T. Smith summarizes the investigations 
ot many intervening years in a paper on Hamilton’s characteristic function and shows that in 
principle the problem of forming the function for a system of coaxial surfaces can be solved. 
He expresses the opinion that modern calculating machines are capable of dealing effectively 
with the numerical aspect of the work, though we have yet to see a paper on the coding of an 
electronic machine for exploiting the remarkable properties of the Hamilton Function. 

It was much to be desired that Smith’s earlier work on this field (1921-27) should have 
been adapted for practical calculations, and the explicit formulae of J. Korringa and 
W. G. Stephan reported by A. S. van Heel now make this more readily possible; they have 
been adapted for calculations of fifth order aberrations, and include differential formulae 
for the effects of variations of axial separations of surfaces and refractive indices of the media. 
The paper gives some interesting curves showing the comparison between aberrations 
calculated by third order formulae, by third and fifth order combined, and the total aber- 
It will be inferred that ray-tracing will still remain a sine qua non for the design of 
systems of large aperture and field! But the work does indeed represent a stage in the 
development of optical calculations. ‘The contemporary work, published in the Proceedings 
of the Physical Society, of H. Buchdahl will be remembered. Abstracts of work on similar 
lines by M. Herzberger and by P. P. Platzeck appear in the book; and H. A. Hopkins 
has an interesting paper suggesting the calculation of the aberration coefficients of rays 
associated with a trigonometrically traced principal ray; this has been more recently 
followed up by W. Weinstein. Space does not permit any comment on other papers in 


this section. 


ration. 
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Turning now to the physical characteristics of optical images, the study indispensable for 


the better knowledge of tolerances in design, we have papers on geometrical aberration — 


patterns by M. Herzberger, A. C. S. van Heel and E. Durand; and (in a later section) 
on aberration diffraction effects by F. Zernike with B. R. A. Nijboer, also by H. H. Hopkins 
and others. The integrator designed by A. Maréchal receives only a short notice. 

A valuable section deals with the possibilities of aspheric surfaces. J. W. Perry 
discusses inter alia some war-time attempts at the correction of photographic lenses for 
off-axis aberrations by figuring of two or more surfaces; B. O’Brian and R. E. Hopkins 
describe new attempts at the ‘ dropping’ technique; and H. Rinia with P. M. van Alphen 
describe their technique for the production of Schmidt plates etc. by the drying of moulded 
gelatine films. F. Dourneau describes a machine for deep parabolic surfaces which is the 
revival of a method known to Newton. One has the impression here of increasing possi- 
bilities. Chapters on researches in the theory of diffraction and in phase contrast include 
some material which is not yet well enough known in England, e.g. the calculation and 
manufacture of amplitude filters which suppress the lateral fringes in spectral lines. ‘This 
topic is discussed by P. Jacquinot, P. Boughon, and Brigette Dossier; also by G. Lansraux. 
The application of phase contrast to observe the irregularities of a polished glass surface is 
described by B. Lyot. This will be of much interest to those who have taken an interest in 
the design of the coronagraph. F. Zernike’s discussion of the phase-contrast microscope is, 
of course, now well known from papers available in English. 

The volume concludes with a group of papers on “‘ The Intervention of the Receptor ”’. 
The researches of A. Arnulf and his collaborators on the capacity of the eye to resolve 
gratings of different fineness and contrast are now well known, and he contributes here a 
short paper discussing the general plan of the work. V. Ronchi summarizes and fore- 
shadows a series of investigations also dealing with the effect of the receptor on resolving 
power. ‘The final paper of the whole collection (contributed by J. L. Tearle) describes the 
work of the Kodak Research Laboratory, Harrow, on photographic resolving power. 

It has not been. possible in the course of a short review even to mention many other 
interesting articles. The flavour of the whole collection is (as the title conveys) theoretical, 
and there are comparatively few descriptions of actual apparatus. It is interesting to 
speculate how far the intense intellectual effort reflected in this work will affect the main 
developments or changes which may be apparent in these coming years. It is unfortunately 
true in optics that much which has been written during the last thirty years has had (so far): 
no effect whatever in technical practice, but delayed action in such a difficult field is not 
unexpected. Some of the more popular modern topics, e.g. the ‘ meniscus plate ’ systems, 
anti-reflection films, interference filters, reflecting microscopes and other things, find no 
echo here; nevertheless in many fields this book will serve, as records of few other symposia 
can, to give some impression at least of the general optical outlook at the close of the second 
world war. Its forty-seven papers and summeries (though of varying merit and value) 


accompanied by many useful references, will stimulate thought and serve to suggest future 
avenues of progress. L. C. MARTIN. 


Introduction to Physical Optics, by J. K. Ropertson. Pp.x+512. 3rd Edition. 
(New York: D. van Nostrand Co., Inc., 1948.) 22s. 


This is the eighth printing of the third edition of an already established textbook on 
elementary physical optics. It therefore needs little comment from a reviewer, and in 
particular the contents of the book will not need listing here. Because there is so little 
that needs saying about the book in general, the following comments on particular points 
should not be regarded as unnecessarily carking. 

On p. 12 occurs the sentence “The modern position, then, is that light has a dual 
nature”. ‘This is epistem logically unsound, and rather confusing to the student. Light 
has both wave and corpuscular attributes, but light has a single nature. Moreover the 
‘modern position’ harmonizes the two attributes in one theory. 

On p. 92 it is suggested that the spherical aberration of a single lens is a minimum 
when the deviation is the same at the two surfaces. 'This is an old error, and the author 
may perhaps be excused. It should, nevertheless, be recorded that the rule is not true, 
except in the trivial case of an equiconvex lens having unit magnification. 


X 
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On p. 274 the author characterizes Fresnel diffraction by saying that (a) ‘ the intensity 
at any point is the resultant of disturbances coming directly to that point from all parts of 
the exposed wavefront’, while in Fraunhofer diffraction (b) ‘a lens is placed beyond the 
aperture or obstacle, and the diffraction pattern is examined in the plane where a sharp 
image of the source would be formed in the absence of the aperture or obstacle ’. 

What distinction is made in (a) is difficult to see, while (b) really describes the conditions 
of observation. 

The real distinction is as follows: At a point near the focus of a wave, the secondary 
disturbances arrive with only small phase differences. Hence the whole of the exposed 
part of the wave is effective, and it is necessary to integrate over the whole wavefront. 
This is Fraunhofer diffraction, which, in the case of a plane wavefront, gives ‘ diffraction 
at infinity ’. 

At a point near to a geometrical shadow, but remote from a region of focus, the phases 
of the secondary disturbances arriving at the point vary rapidly with the positions of the 
points on the primary wavefront at which they originate, except for those points on the 
wave near to the diffracting edge. Hence it is only the part of the wave in the neighbour- 
hood of the edge which contributes effectively to the disturbance at a point near to the 
geometrical shadow. The diffraction fringes then tend to follow the outlines of the 
diffracting edges. This is Fresnel diffraction. 

Speaking of Faraday and the conception of tubes of forces, the author suggests that it 
was ‘ impossible for him, as for most Anglo-Saxon minds, to think of an electrified body 
repelling or attracting another at some distance from it, with nothing taking place in the 
intervening space ’. It came as a surprise to at least one reader to know that such a thing 
as the ‘ Anglo-Saxon mind’ existed. The sentence smacks of woolly racial theory— 
harmless in this case, it is true. 

It is to be hoped that these few blemishes can be removed in a next edition of the book, 
which, for the most part, makes an excellent introduction to an interesting subject. 


Die Brechzahlen einiger Halogenidkristalle, by H. Hartinc. Pp. 25. Ist 
Edition. (Berlin: Sitzungsberichte der deutschen Akademie der Wissen- 
schaften, 1948.) 2.50 DM. 


Using an interpolation formula due to Ketteler, the refractive indices of a series of 
halides have been calculated. The formula contains five constants, which were determined 
independently for the two ranges 1083-03-390-64 and 404-66-239-98 my, using widely 
spaced lines. "Tables give the calculated refractive indices for the range 1120-238 mp, 
with intervals of 10 my for 1120-400 my and 2 mu for 400-238 mu. For each crystal, 
a large number of indices was measured and the results compared with the calculated values. 
The error rarely exceeds 510-5. The crystals included are NaF, LiF, CaF,, NaCl, 
PGK Br KT; 

These tables and measurements should prove invaluable to a great number of workers 
in many fields. H. H. H. 


Introduction to Electron Optics, by V. E. Cossett. Pp. xvi+293. 2nd Edition. 
(Oxford: Clarendon Press, 1950.) 27s. 6d. 


This book needs no recommendation. In the four years which have passed since the 
first edition it has gained well-deserved popularity with students and teachers alike. ‘The 
second edition is not essentially altered. It is noted with regret that Figure Dy) has not 
been corrected, and that the plates illustrating electron microscopes and electron micro- 
graphs have been retained, though they do not illustrate any longer the most modern models 
or achievements. Nor is it surprising, in view of the rapid development in this field, 
that the sections on particle accelerators are considerably out of date. It is hoped that 
this will be remedied in the third edition, which may require also some re-balancing of the 
theoretical parts, with more emphasis on Hamiltonian. methods, which. are now represented 
only in an excellent Appendix, also on wave optics, and on the optics of ponmaccnered 
systems, which have become increasingly important since it has been realized that Workshop, 
accuracy is a major limitation of electron optical systems. But though there is room for 
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improvement, this is of little importance to the student, who will find in Cosslett’s book 
an excellent account of what he ought to know first, before venturing into the more advanced 
fields of the subject. The research worker on the other hand will find in the newly added 
Notes a short, but almost complete and up-to-date collection of the most important recent 
advances in electron optics. D. GABOR. 


Theory of Hearing, by E. G. Wever. Pp. 484. (New York: Wiley & Sons; 
London: Chapman & Hall, 1949.) $6.00. 


The appearance of this volume will be welcomed not only by physicists but also by 
physiologists, anatomists and psychologists. It is most appropriately dedicated ‘ to my 
associates in the study of the ear’. Dr. Wever has been engaged in auditory research 
throughout his career, and is a well-known leader in this field. "Twenty years ago, when 
his work commenced, an interesting situation had arisen—two rival theories of hearing, 
the ‘place’ and ‘frequency’ theories, held the field on the physical side at a time when 
the new electro-physical techniques which provided basic evidence on sensory nerve 
action were being developed. Under these circumstances Dr. Wever and Dr. C. W. Bray 
carried out their classical experiments on auditory nerve impulses, leading to the discovery 
of the electrical potentials of the ear—work which. laid the foundations of much further 
research in many laboratories. 

As one might expect, therefore, Dr. Wever’s book deals principally with his own 
researches and theories in considerable detail, but nevertheless it does not neglect to refer 
to the work of many other investigators on the complex problem of hearing. ‘The volume 
is divided into three parts containing, in all, 18 chapters. The first part deals with the 
historical side of the subject, the ‘ classical’ theorics, and is mainly concerned with the 
Helmholtz resonance theory, the ‘ place’ theories and the ‘ frequency’ (or ‘ telephone ’) 
theories. These theories, some plausible and some rather fantastic, are explained in four 
very interesting chapters. The sccond part of the book includes four chapters on (1) 
modern developments of the classical theories, (2) auditory nerve responses, (3) frequency 
representation in the auditory nerve and (4) the ‘ volley’ principle. In part three the 
“ volley’ principle is developed and further evidence on cochlear localization is discussed. 
In this part also the various auditory characteristics and their bearing on the theories of 
hearing are considered, viz. sensitivity, loudness and fatigue, pitch discrimination, auditory 
abnormalities (deafness, etc.), transient phenomena, and binaural effects. 

The chapters on auditory nerve responses are particularly interesting in showing that 
the alternating fluid pressures in the ear, which themselves faithfully represent the sound 
waves, are communicated to the hair cells, and the cochlear potentials (which are traced 
to the hair cells) are representative of the form of the stimulus—the sound wave. ‘The 
cochlear potentials are observed over a wide range of magnitudes from the smallest 
experimentally demonstrable, a tenth of a microvolt, to as much as a millivolt—a range of 
10,000-fold—the potential being proportional to sound pressure. The description of the 
remarkable experiments involved in such investigations makes extremely interesting reading. 
As a mere physicist the reviewer ‘ takes off his hat’ to the anatomical physiologist who 
can carry out successfully such delicate operations, as, for example, the demonstration 
(by display on a cathode-ray tube) that the electro-physiological response of the ear lies 
in the hair cells of the cochlea. 

The book is full of interest and leaves one with the feeling that the physicist is lagging 
behind the anatomist. ‘There is clearly plenty of scope for further research—particularly 
of the ‘ combined operations’ variety in which the physicist must work closely with the 
physiologist and anatomist. For the benefit of physicists studying the book it would be 
an advantage if, in a future edition, the principal dimensions and other physical 
characteristics of the important parts of the ear could be tabulated or collected together 
in an appendix. Such information is at present scattered throughout the pages of the book 
and is not easily available for reference. 

The book is well written and inspiring. It displays the great complexity of the subject 
and clearly shows that there is still much to be learnt. ‘The volume can be recommended 
unreservedly to all those interested in the subject of hearing in man and animals. 


A. B. WOOD. 
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‘Counting Tubes, by S. C. Curran and J. D. Craces. Pp. xi+238. (London: 
Butterworth’s Scientific Publications Ltd., 1949.) 35s. 


Books on ion counters are few and far between. There is still an urgent need for a 
comprehensive and up-to-date treatment of the subject. The present book deals almost 
exclusively with the gas-filled type of counter, with the exception of one chapter which gives 
a brief discussion of some other types such as scintillation and crystal counters. The early 
part of the book summarizes the relevant background information on electrical discharges 
and of the ionization produced by fast particles, and then goes on to analyse the behaviour 
of counting tubes in the different voltage regimes. There follow three chapters which deal 
with the efficiency of counters, their construction, and some of the auxiliary circuits to be 
used with them. The remainder of the book comprises the chapter on special types of 
counter already mentioned and some chapters dealing with representative applications of 
counter tubes in different fields of research. 

The authors have included a great deal of useful information on the different phenomena 
which occur in counter tubes and there is a good selection of references to the original 
literature. The book was, however, found to be a difficult one to read. This was due in part 
to insufficient care in the definition of terms and symbols in a number of places and, in part, 
to an over-condensed treatment of the material. In fact—rightly or wrongly—the impres- 
-sion was gained that the book had been written rather hurriedly. Nor, unfortunately, is the 
text free from error. In the section dealing with the resolving power of counters, the 
argument leading to equation (40) on page 95 is definitely fallacious and it is not correct 
for the authors to state, as they do on the following page, that it is difficult to decide whether 
this or equation (42) is the more nearly correct. Certainly confusion has existed on this 
point in the past, but it has now been clarified and the matter is adequately dealt with in the 
literature. 

‘The authors state in their preface that they hope that “‘ this book will appeal both to 
graduates about to commence experimental work and to undergraduate students who find it 
necessary to use Geiger counters as instruments. In other fields such as medicine and 
chemistry an increasing number of workers are finding it necessary to use counters in their 
researches and we hope that this book will be of assistance to them, if only as a source book ”’. 
It is felt that the potential value of the book has been diminished by trying to cast the net too 
wide. ‘Those whose interest is primarily in the use of the counter as a tool are liable to find 
their understanding of the broad principles of its behaviour obscured by a fog of detail 
-concerning little-understood processes in the discharge, while to those for whom these 
details are of importance, the chapters dealing with sundry applications of counter tubes are 
likely to be redundant. It is to be hoped that the authors will find it possible to remedy 
these defects in a second edition. J. L. M. 


The Strength of Plastics and Glass, by R. N. Hawarp. Pp. vilit+245. 
1st Edition. (London: Cleaver-Hume Press, Ltd., 1949.) 30s. 


Although the science of glasses and plastics is a comparatively new subject, a considerable 
amount of theoretical and experimental work on the strength properties of these materials 
has been reported in the literature over the past thirty years. The need for collecting together 
the results of this work into one volume has been felt for some time by specialists in these 
subjects. For this reason alone, the appearance of this monograph by Dr. Haward is to be 
welcomed. However, there is much of fundamental interest in the book which will appeal to 
a wider scientific public and particularly to the rheologist outside the field of applied science, 
for Dr. Haward has whenever possible directed attention to the contrasting behaviour of 
plastics and glass under stress. Glass may be regarded as an ideal isotropic solid whose 
deformation is almost entirely Hookean, while relationships based on the classical theory of 
-elasticity are not generally applicable to plastics. 

A major difficulty arises in comparing the results of work on strength properties from 
different sources, because the measured basic quantities frequently differ in dimensions; for 
this reason it is sometimes possible to draw misleading conclusions because there appears to 
be a correlation between two such quantities simply because some other relevant factor is not 
changing sufficiently to affect the correlation. These difficulties are explained clearly in the 
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introductory chapter. The next five chapters deal in turn with the static strength of plastic 
materials, the static strength of glass, the deformation of these materials under stress, 
hardness and related properties and the impact strength of glass and plastics. The final 
chapter which concerns the properties of some composite materials, such as those in which 
plastics and glass are combined in the form of laminates, is of particular interest to the 
technologist. 

These chapters take the form of a presentation of the results of the most important 
published work and, where necessary, reasons are advanced for the apparent contradictions 
between the results obtained by different workers. Where theoretical explanations are in 
doubt, the author confines himself to theories accompanied by the results of experiment. 

The subject-matter of each chapter covers a wide field but much of the detailed descrip-- 
tion of the behaviour of individual plastic materials will be of interest to the technologist 
rather than to the pure scientist. So far as the glass technologist is concerned, the inclusion 
of a survey of the behaviour of glass under thermal stresses would have enhanced the value 
of the book. 

The bibliography at the end of each chapter comprises nearly 400 references, mostly to 
work published within the past ten years. Particular mention should be made of the 
considerable assistance which the reader derives from the clear presentation of the line 


diagrams. The author, printer and publisher are to be congratulated on the production of 


an easily readable work. R. E. BASTICK. 


Recent Advances in Radio Receivers, by L. A. Moxon. Pp. ix +183. 1st Edition. 
(Cambridge: University Press, 1949.) 18s. 


The advances which have been made in radio receiver design during the last decade have 
arisen mainly from wartime requirements and the postwar expansion of the television service. 
The present book in the series Modern Radio Technique is primarily concerned with design 
for low noise factor and good pulse response, which are among the most important character- 
istics of a radar or television receiver. 

Chapter 1 is a brief introduction explaining the scope of the book. Chapters 2 to 5 are 
concerned with the definition, calculation and measurements of noise factor and strike a 
satisfactory balance between the rigour required in the formulation of the fundamental 
concepts and the many practical considerations which determine the performance of actual 
amplifiers. ‘The influence of feedback, partition noise and transit-time loading are discussed 
and emphasis is directed to the necessity of compromise between the conflicting require- 
ments of a wide bandwidth and a low noise factor. R.F. amplifiers, especially triodes, for 
frequencies up to 600 Mc/s. and mixers for higher frequencies form the basis of the input 
circuit considerations. 

Chaptcr 6 is concerned with the design of intermediate-frequency amplifiers and, in 
particular, the use of coupled, staggered and cascaded circuits to provide a desired bandwidth 
with sufficient uniformity of gain and time delay to give an acceptable pulse response. The 
chapter concludes with a brief discussion of overloading, unwanted feedback and the choice 
of intermediate frequency; much useful advice is given in a small space. 

The last three chapters comprise a miscellany of recent features in the design of radio. 
receivers, the selection of which was, inevitably, determined by the author’s own interests. 
The Appendix which concludes the book outlines the analysis of the frequency response of 
multi-stage I.F. amplifiers. 

The errors in the text are relatively few and unimportant, and the diagrams, with one or 
two exceptions, are clear and useful. 

The book can be recommended as a primer on receiver design. For more detailed 
information on some problems the reader would turn naturally to Volume 18 of the Radiation 
Laboratory series on Vacuum Tube Amplifiers but in the field of British books, Mr. Moxon’s. 
contribution has no rival at present. R. E. BURGESS.- 
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Le Vide et ses Applications, par L. Dunoyer. Pp. 114. 1st Edition. (Paris: 
Presses Universitaires de France, 1950.) 90 fr. 


Most physicists will have heard of Dunoyer because of his book Vacuum Practice, 
published in 1926: it was a pioneer work on the subject. His recent little book is a 
contribution to the series “‘ Que sais-je?’’ which includes over four hundred titles. A 
small, paper-covered volume with only twenty-six diagrams, it is not to be expected that 
it would cover the subject adequately, or give those practical details of apparatus and 
materials so useful to the vacuum technician. Nevertheless, it is written so well, and so 
delightfully, and the material included is selected with such great experience, that it can 
be read with profit by all interested in vacuum physics provided they have a reasonable 
knowledge of French. Especially noteworthy are the clarity and conciseness of the earlier 
parts on kinetic theory, adsorbed gas layers, and the flow of gases through orifices and 
tubes. Vacuum pumps, gauges, speed measurements, vacuum stoving, getters, and some 
industrial applications are all discussed in simple terms, with an occasional flash of humour, 
yet in a manner which convinces the reader that this is the work of a scientist with a 
wealth of experience from which to draw. The text is singularly free from error: the 
reviewer discovered only one small mistake. Some criticism might be levelled at the sections 
on electron-tube manufacture, the production of optical mirrors, vacuum metallurgy, and the 
preparation of penicillin and of dehydrated orange juice. Here the lack of illustration 
and detail is such that a rather vague impression must be left in the minds of those not 
already familiar with these topics. Again, it is a pity that nothing is said about leak-hunting. 

It is a book particularly suitable for the physics student, or for the layman with some 
mathematical knowledge and an intelligent interest in science. At the same time, even 
the physicist with years of experience in vacuum practice will frequently come across 
ideas and descriptions which throw new light on many of his problems. J. YARWOOD. 
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ABSTRACTS FOR SECTION A 


Cascade Theory including Ionization Loss, by L. JANossy and H. MEssEL. 
ABSTRACT. The Bhabha-Chakrabarty solution of the cascade equations with 


ionization loss is discussed, and a simplified method of numerical evaluation suggested. 


The Nuclear Interaction Length of the Particles in Penetrating Cosmic-Ray Showers, 
by K. H. Barker and C. C. BuTuer. 


ABSTRACT. The work of Rochester and Butler, published in 1948, and that of Butler, 
Rosser and Barker, published in 1950, has been continued. In particular, the interaction 
properties and the spectrum of the penetrating particles are considered. The interaction 


length of the shower particles with momenta greater than 10° ev/c. is found to be 18 cm. 
lead (420%). 


The Importance of Gas Scattering in Particle Accelerators, by L. RippIFORD. 


ABSTRACT. In electron and proton synchrotrons there exists a minimum useful vertical 
aperture of the vacuum chamber, and for apertures less than this the beam intensity falls. 
to zero. It is considered, on the basis of a comparison of theory and experiment, that this 
phenomenon is in the main due to free oscillations which arise from gas scattering. The 
final beam size in proton synchrotrons may be substantial for this reason also, but for 
cyclotrons the large spread which is known to exist must be attributed to some other 
phenomenon. 


An Investigation of the Secondary Electron Spectrum of 8 Au, by P. E. CAVANAGH, 
J. F. Turner, D. V. Booker and H. J. DuNSTER. 


ABSTRACT. ‘The secondary electron spectrum of 18Au has been studied using very 
strong sources, up to 0-3 curie, in a magnetic lens spectrometer. In addition to the 
well-known y-ray at 0:411 Mev., evidence of two much less intense higher energy y-rays 
at 0-67 Mey. and 1-09 Mey. has been found, both from the Compton and photoelectric 
spectra. he intensities of these two y-rays with respect to the one at 0-411 Mev. have 
been measured. ‘Two photoelectric peaks at low energy, in about the region expected 
for Ka and Kf x-radiation from elements in the neighbourhood of gold, have been 
examined with radiators of different atomic number. The radiation responsible for these 
peaks has been subjected to critical absorption measurements in the spectrometer. Its 
behaviour is compatible with the Ka and Kf x-radiation of gold, and no evidence of 
platinum x-radiation, such as would suggest K-capture, has been found. 


A Method for Determining Uranium and Thorium in Rocks by the Nuclear Photo- 


graphic Plate, by J. W. BREMNER. 
ABSTRACT. A nuclear photographic emulsion which has been in contact with material 
containing traces of uranium and thorium shows on processing «-ray stars arising from 
atoms which have undergone several successive disintegrations during the exposure. ‘The 
relative frequency of the various multiple stars, which depends on the uranium/thorium 
ratio, is here evaluated on simple assumptions to show the possibility of determining these 
elements by star counting. 


The Planar Vibrations of Ethylene and Tetra-deutero-ethylene : A Critical 
Analysis of the Potential Function, by P. 'TORKINGTON. 
ABSTRACT. A careful analysis of the planar modes of vibration of ethylene and 
ethylene-d, is carried out in order (a) to determine finally the value of the C=C stretching 
force constant in these key molecules, (b) to find what interactions, particularly those 
between the two ends of the molecules, are real, and (c) to determine the vibrational forms 
and the potential energy distributions in the normal modes, and how they are affected 
by altering the potential function. Solutions are obtained both with and without correction 
for anharmonicity, and estimates are made of the effects of other corrections which may be 
necessary. The limits found for the C=C stretching constant are 9-336, 9-388 x 10° dyne/cm., 
with a most probable value of 9-343, the carbon atoms being uncorrected for anharmonicity. 
The zero-order constant is unlikely to be greater than 9:57. There is detectable real inter- 
action between the two ends of the molecule, for both stretching and bending coordinates. 


Interactions are discussed. 
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On the Elimination of Divergencies from Classical Electrodynamics, by Suraj N. 
GUPTA. 


ABSTRACT. It is shown that the divergencies in classical electrodynamics may be 
removed by subtracting the self-energies of isolated electrons by means of a modification 
of the usual expression for the force density in an electromagnetic field. ‘The results 
obtained are equivalent to those of Dirac and Wentzel. 


On the Renormalization of Quantum Electrodynamics, by J. C. WaRD. 


ABSTRACT. ‘The subtraction procedure of Dyson is modified in order to eliminate 
a certain difficulty in the renormalization programme, namely the so-called ‘ b’ divergencies. 
The conclusions of Dyson concerning the finiteness of the renormalized S matrix are 
confirmed. 


Equations of Motion in General Relativity, by A. PAPAPETROU. 


ABSTRACT. This paper contains a new derivation of the equations of motion of bodies 
moving slowly in their (weak) gravitational field, up to terms of the order (v/c)?.__ In contrast 
to the method developed previously by Einstein and co-workers, in which only the field 
outside the bodies has been considered, in the present paper a new method is used in which 
the field outside as well as inside the bodies, and the internal structure of the bodies, is taken 
into consideration. The equations of motion derived by this new method are identical with 
those obtained by Einstein and co-workers, but the new method is essentially simpler and 
the necessary calculations much less lengthy. 


Resonant Nuclear Scattering of Gamma-Rays : Theory and Preliminary Experi- 
ments, by P. B. Moon. 


ABSTRACT. Since the lower excited states of nuclei have very small widths (<1 ev.), 
resonant scattering of gamma-rays requires precise matching of the energy available from 
the gamma-ray with the energy necessary to excite the scattering nucleus. 

Resonant scattering should be observable if (1) the emitting and scattering nuclei are 
of identical type, (2) the gamma-transition goes to the ground state and (3) the source and 
scatterer are given such a relative velocity that Doppler effect restores the energy lost by the 
gamma-ray to nuclear recoils. Thermal velocities of the emitting and scattering nuclei 
broaden and correspondingly weaken the resonant scattering peak, and the cross section at 
the optimum speed of 32E/A cm/sec. is 3:6 x 10°F /E*)(A/T)1/? cm?, where E and FI are 
the energy and intrinsic width of the excited state in electron volts, J the isotopic abundance 
of the resonantly scattering isotope, A its atomic weight and T the absolute temperature. 

Preliminary experiments have been made with the 0:411 Mev. radiation from the nucleus 
18FIg, the source being carried by a high-speed rotor up to a speed of about 7 x 10% cm/sec. 
and the scatterer being liquid mercury (10% 1°8Hg). A small but apparently significant 
increase of scattering was found, corresponding to a width F of the order of 10~° ev. 

No such increase was observed with !*"T'a gamma-rays scattered from tantalum carbide. 

The negative result for 1*"T'a and the positive result for 1*%Hg are consistent with the 
latest information about the life-times of the excited states concerned, viz. 1:1 x 10 sec. 
for 18"T'a and less than 2 x 10~1° sec. for 1°8Hg. 


Detection of -Mesons and other Fast Charged Particles in Cosmic Radiation, by 
the Cerenkov Effect in Distilled Water, by J. V. JELLEY. 


Me 

ABSTRACT. 'The Cerenkov effect has been established for the passage of cosmic-ray 
#-mesons through distilled water. The light pulses were observed using a photomultiplier 
and fast amplifier; the experiments were carried out using the detector alone, and in a 
coincidence-telescope arrangement. ‘The ratio of the light emitted in the downward and 
upward directions was measured as a function of discriminator bias. In addition measure- 
ments were carried out at various depths of water and with various orientations of the 
detector with respect to the vertical. Estimates of the absolute efficiency were made and 
other characteristics of the detector are discussed. 
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Figure 2. Slit images showing the variation Figure 10. Shadowgraph of two 104 
of spot size and quality with grid platinum wires at right angles, 
aperture—anode distance d. Exposure supported by a platinum grid. 
times: 5 min. and 25 min. Slit width Direct magnification about 20 times. 


7 w, magnification 10 times. 


Figure 8. Photograph of tube. 
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FUSED QUARTZ 
X-RAY 
SPECIMEN TUBES 


We are pleased to announce that we can 
now supply transparent VITREOSIL (pure 
fused quartz) specimen tubes for use in 
X-ray analysis. 


Fused quartz is extremely transparent to 
X-rays, and such tubes are ‘particularly 
suitable for use in high-temperature X-ray 
cameras or in other instruments in which 
high transparency to X-rays is essential. 


Length Bore Wall Thickness 
2in. 0-25-0-3 mm. 0-035-0-05 mm. 
3 in. 0-25-0:3 mm. 0:035-0-05 mm. 


In addition to the above tubes, we can also 
supply, to customers’ specification, Thermal 
alumina ware supports for high-temper- 
ature X-ray cameras. 
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The Proceedings are divided into two 
parts, AandB. ‘The charge for insertion 
is £18 for a full page in either Section A 
or Section B, £30 for a full page for 
insertion of the same advertisement in 
both Sections. The corresponding 
charges for part pages are: 


Rpage. £9 -5°0., £1500 

tpage £415 0 £8 0 0 
t page £210 0 £4 5 0 
Discount is 20% for a series of six 


similar insertions and 10% for a series 
of three. 


The printed area of the page is 
84” x 54”, and the screen number is 120. 


Copy should be received at the Offices 
of the Physical Society six weeks before 
the date of publication of the Proceedings. 


THE HANDBOOK OF THE 
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OF SCIENTIFIC INSTRUMENTS 
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Orders, with remittances, should be sent to 
THE PHYSICAL SOCIETY 


BINDING CASES 
for the 
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Binding cases for Sections A and B (separate) for 
Volume 63 (1950) may be obtained for 7s. each, 
post free, from the Offices of the Society. 


The Proceedings may be bound in the Society’s green 
cloth for 18s. 6d. each. Journals for binding should be 
sent direct to Messrs. Taylor and Francis 

Red Lion Court, Fleet Street, 

Remittances should be sent to the Physical Soci 
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